Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,718,490 T→G 18.7% D56A (GAC→GCC)  CJ019_01614 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,718,4900TG18.7% 119.0 / 3.0 64D56A (GAC→GCC) CJ019_01614hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/40);  new base G (12/2);  total (22/42)
Fisher's exact test for biased strand distribution p-value = 1.21e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAGCA  >  ACB122/1718358‑1718622
                                                                                                                                    |                                                                                                                                    
aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGGGCGGCAGGGGCGTcc                                                                                                                                    >  1:12390/1‑135 (MQ=255)
      ctGCGCGCACTACCCGCATTGCGGTGGCCAGGCTTGCCCTGCGCGGGGTCATAGGGGAGTTGGAACGGCAGCGGCACGTTGCCTTGGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGc                                                                                                                              <  2:65935/134‑1 (MQ=255)
      aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGGGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                                              <  7:296059/135‑1 (MQ=255)
      aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                                              <  3:373363/135‑1 (MQ=255)
      aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                                              <  7:202106/135‑1 (MQ=255)
         cgcgcACTTACCGCATAGCGGTGGGCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                                                           <  7:24075/135‑1 (MQ=255)
         cgcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                                                           <  7:79312/135‑1 (MQ=255)
         cgcgcACTTACCGCATAGCGGTGGCCAGGCGTGACCGGCGCGGGGGCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                                                           <  5:86/135‑1 (MQ=255)
         cgcgaAATTACCGCATAGGGGGGGCCCGGCTTGACCGGGGGGGGGTCATCGGGGCGAGGCACCGGCCGCGGCACGTGGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGac                                                                                                                           <  4:23378/135‑1 (MQ=255)
          gcgcACTTACCTCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTTGCGCGCCAGGGTCACCGCCGGGGCGTGCAGCcgcgcc                                                                                                                           >  6:267052/1‑132 (MQ=255)
          gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGCCCGGGTCCCGGCAGGGGCGGCCAACAGCgccc                                                                                                                          >  6:22373/1‑132 (MQ=255)
          gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCCGGGGCGGCCCTCAGCGcca                                                                                                                          >  4:153257/1‑135 (MQ=255)
          gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCCCGGCCGGGGCGGCCCTCCGCGaca                                                                                                                          >  5:66975/1‑135 (MQ=255)
           cgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                         <  8:57311/135‑1 (MQ=255)
              aCTTACCGATTAGCGGTGGCCAGGCTGGACGTGCGCGGGGTCATTGGTGACAGGGAACGGCAGCGGCAGGTTGCGTTTGACGGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                                                      <  4:327797/135‑1 (MQ=255)
               cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGGACGGCCGGGGCGGGCCGCAGCGACCCCTcc                                                                                                                     >  8:16813/1‑135 (MQ=255)
                    cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCCGGGGCGGCCCTCCGCGACCCCGCCTGCAg                                                                                                                >  8:202106/1‑135 (MQ=255)
                           cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                                                          <  4:332367/134‑1 (MQ=255)
                           cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCCCGGCCGGGGCGGCCCTCCCCGCCACATCCGTCCGTTccccc                                                                                                         >  2:346922/1‑132 (MQ=255)
                           cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACg                                                                                                         <  2:281646/135‑1 (MQ=255)
                                   aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGGCCCGGCGGGGGCGGCCCTCAGCGCCCCATCCTGCGGGGCCGCGCGCAGGGt                                                                                                 >  3:325256/1‑135 (MQ=255)
                                       ttGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGCCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                                             <  7:130837/135‑1 (MQ=255)
                                       gtgcccTGCGGGGGGGGATGGGTGAGCTGCACCGGACGCGCCACGTGGCGGTTGCCCGGGCGGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGGCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                                             <  7:86637/131‑1 (MQ=255)
                                       gtGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                                             <  4:59614/134‑1 (MQ=255)
                                        tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCCGGGGCGTCCCTCAGCGCCACCTCCTGCAGTTCCCCGCGCCGGGGggcgg                                                                                            >  6:77998/1‑135 (MQ=255)
                                        tGACCTGCGCGGGGTCATAGGTGAGATGCAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcgg                                                                                            <  5:267052/135‑1 (MQ=255)
                                            gtGCGCGGGGTCTTAGGTGACATGTAACGCCAGCGGCAGGTTGCCTTTGACCGGGCCGTCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                                                        <  3:183853/134‑1 (MQ=255)
                                            cTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                                                        <  8:265502/135‑1 (MQ=255)
                                            cTGCGCGGGGGCATAGGTGAGATGGACCGGCCGCGGCACGTTGCCTGTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                                                        <  1:44959/135‑1 (MQ=255)
                                             tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGCGCGGCCGGGGCGGCCCTCAGCGCCCCCGCCTGCAGGTCCCCCCTCCGGGGggcgggcggc                                                                                       >  5:187422/1‑130 (MQ=255)
                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCAGGCGCGCCACGGGCACCGCCGGGGCGTCCATCCGCGACACCTCCAGCAGTTACACGCTCACGGGggcggccggcggccg                                                                                  >  3:107889/1‑133 (MQ=255)
                                                      tCATAGGTTAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCa                                                                              <  3:129386/135‑1 (MQ=255)
                                                       cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGGGGGCGCGCCCGGGGCCCGGCCCGGGGGGGGCGCAGCGACACATCCTGCCGGTACACGCGGAGGGTGGCGGCCGGCGGCCGGGcca                                                                             >  8:244836/1‑133 (MQ=255)
                                                        aTAGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAAt                                                                            <  3:184402/135‑1 (MQ=255)
                                                            gggagaGGGAACGGCAGCGCCACGTTGCCTTTGCCCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                                                        <  2:187350/133‑1 (MQ=255)
                                                            gTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGGACCGCAGGGGCGGCCCTCAGCGCCACAACCAGCCCTTCCCCGCGCAGGGTGGCGGCCCGCGGCAGGGCCATgcgc                                                                        >  3:262796/1‑135 (MQ=255)
                                                            gTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                                                        <  7:244836/135‑1 (MQ=255)
                                                                  tGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCAGGGGCGGCCATCAGCGCCACAACCTGCACGTCCCCGCGCATGGTGGCGGCCGGCGGCCGGGCCATGCGCGGCAgg                                                                  >  6:249962/1‑135 (MQ=255)
                                                                         gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGGCCGGGGCCCGGCCGGGGGGGCCATCCGCGCCACAGCCGGCAGATCCCCGCGCCGGGGGGCGGGCGGCGGCCGGGCCATGCGCGGCAgggggcag                                                            >  8:59111/1‑128 (MQ=255)
                                                                           cggggCAGGTTGCCTTTGCCGGGGCGGGCGGGGGGCGCCTGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                                         <  2:180218/132‑1 (MQ=255)
                                                                           cggggCACGTGGCCTTTGACGGGGCCGGCGTGGTGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGTCAATGCGCTGCAGGTAGAAGACt                                                         <  1:158365/132‑1 (MQ=255)
                                                                           cgCGGCAGGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                                         <  5:344722/134‑1 (MQ=255)
                                                                           aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                                         <  2:31027/135‑1 (MQ=255)
                                                                           aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                                         <  5:189187/135‑1 (MQ=255)
                                                                           aGCGGCACGTTGCCGTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                                         <  4:223204/135‑1 (MQ=255)
                                                                             cGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCGCCCCAACCAGCAGTTCCACGCGCAGGGTGGCGGCCGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTc                                                       >  6:186783/1‑135 (MQ=255)
                                                                                         tttGCCCGGGCCGGCTGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCt                                           <  5:372116/135‑1 (MQ=255)
                                                                                            gACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTg                                          <  3:243589/133‑1 (MQ=255)
                                                                                             aCCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGcc                                       <  7:19258/135‑1 (MQ=255)
                                                                                              ccGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGGAGAAGACTTCACCGTCCAGGCTGGcc                                       <  7:19263/134‑1 (MQ=255)
                                                                                                  gccggccgGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt                                  <  5:34358/135‑1 (MQ=255)
                                                                                                  gccggccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt                                  <  8:92066/135‑1 (MQ=255)
                                                                                                    cggccGGGCGCGCCAGGGTCACGGCCGGGGCGTCCATAAGCGACCCAACCTGGAGTTCCACGCGCCGGGTGGCGGCCGGCGGCCGGGCAATGCGCGGCCGGGTGAACACTTCCACGTCCAGGGGGGCCTGGGtgg                                >  6:207347/1‑135 (MQ=255)
                                                                                                            cgcgcCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCGGGTTGGGTGGGGtg                        <  4:152606/135‑2 (MQ=255)
                                                                                                              cgcCAGGGTCACGGCAGGGGCGTCCATCAGCGACACAACCTGCAGTTCCCCGCTCACGGTGGCGGCAGGCGGCAGGGCAAAGAGCTGCAGGTAGAAGACTTAACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTg                      >  8:8769/1‑135 (MQ=255)
                                                                                                                  tggggCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGGGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTg                  <  1:192161/134‑1 (MQ=255)
                                                                                                                   gggTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCt                   >  3:132180/1‑133 (MQ=255)
                                                                                                                   gggTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTgg                 <  3:178268/135‑1 (MQ=255)
                                                                                                                    ggTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGa                <  6:283245/135‑1 (MQ=255)
                                                                                                                        aCGGCAGGGGCGGCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAg            <  5:43713/135‑1 (MQ=255)
                                                                                                                          ggCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGc          <  6:66975/135‑1 (MQ=255)
                                                                                                                             aGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGc       <  5:77998/135‑1 (MQ=255)
                                                                                                                                  cGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAGCa  >  5:228837/1‑135 (MQ=255)
                                                                                                                                  cGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAGCa  >  6:104355/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                    
AGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAGCA  >  ACB122/1718358‑1718622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: