Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,655,934 A→C 5.3% intergenic (‑444/+8) iscR ← / ← trmJ isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,655,9340AC 5.3% 604.8 / 5.1 224intergenic (‑444/+8)iscR/trmJisc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (131/75);  new base C (4/14);  total (135/89)
Fisher's exact test for biased strand distribution p-value = 8.24e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

CGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGTTTTCACGG  >  CP009273/2655785‑2656083
                                                                                                                                                     |                                                                                                                                                      
   ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                                                                    >  1:1214487/1‑150 (MQ=35)
   ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                                                                    >  2:2408743/1‑150 (MQ=35)
   ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                                                                    >  2:4313415/1‑150 (MQ=35)
       cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                                                                                >  1:483002/1‑150 (MQ=32)
       cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                                                                                >  1:4360116/1‑150 (MQ=32)
       cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCGACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                                                                                >  1:483060/1‑150 (MQ=17)
          ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                                                                             >  2:1515477/1‑150 (MQ=33)
           tCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAt                                                                                                                                            >  1:2427035/1‑150 (MQ=34)
             tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCCACGGCGGTTATTc                                                                                                                                          >  2:3075856/1‑150 (MQ=25)
             tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                                                                                                                                          >  2:730975/1‑150 (MQ=35)
             tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                                                                                                                                          >  1:1442958/1‑150 (MQ=35)
              gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCg                                                                                                                                         >  2:1420427/1‑150 (MQ=37)
               ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                                                                                                                        >  1:4262962/1‑150 (MQ=37)
                  aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCt                                                                                                                                     >  2:1330879/1‑150 (MQ=38)
                       aGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTa                                                                                                                                >  1:2967620/1‑150 (MQ=255)
                       aGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTa                                                                                                                                >  1:3622927/1‑150 (MQ=255)
                           cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                                                                                                                            >  2:263702/1‑150 (MQ=255)
                            tGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt                                                                                                                           >  2:3511928/1‑150 (MQ=255)
                            tGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt                                                                                                                           >  2:1888184/1‑150 (MQ=255)
                              aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                                                                                         >  1:1033055/1‑150 (MQ=255)
                                 aGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                                                                                                                      >  1:388085/1‑150 (MQ=255)
                                     gCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCTCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                                                                                                  >  1:2393436/1‑150 (MQ=255)
                                     gCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                                                                                                  >  1:4063948/1‑150 (MQ=255)
                                      cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                                                                                                 >  1:792181/1‑150 (MQ=255)
                                       tttccgcgGCGCATCTGGCACCGGCTGTCGGTTGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                                                                                                                <  2:1059046/146‑1 (MQ=17)
                                       aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                                                                                                                >  2:4252632/1‑150 (MQ=255)
                                       aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                                                                                                                >  1:3072879/1‑150 (MQ=255)
                                        ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                                                                                                               >  1:578904/1‑150 (MQ=255)
                                             cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTCCCTGTATTCTGCTGCTCACTAg                                                                                                          >  2:4279451/1‑150 (MQ=255)
                                             cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCGTATCCGACCGACGGTTCTGTTCACTGTAGGCCTGATCAGACGCCTTACGCGTCGCATCTGGCACCGGCTGTTATTCGGCCTTGTGTCCTGTATTCTGCTGCTccctag                                                                                                          >  2:4278115/1‑145 (MQ=17)
                                                 gcctCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGc                                                                                                      >  2:757070/4‑150 (MQ=255)
                                                 gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCACTAGAAGc                                                                                                      >  2:2952109/1‑150 (MQ=255)
                                                  cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGcc                                                                                                     >  2:3251540/1‑150 (MQ=255)
                                                      aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                                                                                 >  2:4004340/1‑150 (MQ=255)
                                                      aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                                                                                 >  1:1827348/1‑150 (MQ=255)
                                                          aaggggtgtgGGTTGGGGGGTGACCGCATTAACCGACCTCCGGTTCTGTTCACTGTCGGCCTTATACGACGCTTTCCGGGGCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATccc                                                                                             <  1:1641937/144‑1 (MQ=255)
                                                                 gTCGGAGGCGGCGTGACCGCCTTACCCGACCTCGGTTTATGTACACTGTCGGCCTTATAAGACGCATTACGCGTCGCATCAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAgg                                                                                      <  2:797012/150‑1 (MQ=25)
                                                                    ggATGCGGGTTGACCGCCTTATCCGCCCTAGGTTTCTGTAAACTGTAGGCCTGAAAAGACGCATTCCGCGTCGAATCAGGAAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGata                                                                                   <  2:543081/150‑1 (MQ=255)
                                                                    ggATGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTCGGCCTGAAAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGata                                                                                   <  2:2319593/150‑1 (MQ=255)
                                                                    ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGata                                                                                   >  2:1356125/1‑150 (MQ=255)
                                                                    ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGata                                                                                   >  1:1052164/1‑150 (MQ=255)
                                                                       tGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc                                                                                >  2:4289644/1‑150 (MQ=255)
                                                                         gggcgTGACCGCATAATCCGACCTACGTTTCTGTACACTGTAGGCCTGATAAGCCGAATTCCGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                                                                              <  1:2108604/149‑1 (MQ=37)
                                                                          ggcgTGAACGCCTAATCCGACCTCGGGTTCTGTTCACTGTCGGCGTGATAAGACGCATTACGCGTCGAATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                                                                             <  2:3410103/150‑1 (MQ=25)
                                                                           gcgTGACCGCCTTAGCCGACCTTGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCCGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                                                                            <  2:3752093/150‑1 (MQ=25)
                                                                            cttGAACGCCTAACCCGACCTCCTTTTCTTTTCTCTTTAGCCCTTAAACCACCCATTCGGCGACGCATCAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                                                                           <  1:3413211/148‑1 (MQ=255)
                                                                            cgTTAACCCATTATCCGACCTCGGTTTCTGTTCAGTGTTGGCCTGAAACGACGATTTACGCGTCGCACCAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                                                                           <  2:1171009/150‑1 (MQ=255)
                                                                              tGACCGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTTATAAGACGCATTCCGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt                                                                         <  2:1909976/150‑1 (MQ=37)
                                                                               taaCGCCTTACCCGACCTCGGGTTCTGTTCAGTGTAGGCCTTATAAGCCGCATCCGGCGTCGAATCCGGCCACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTg                                                                        <  1:3727187/149‑1 (MQ=255)
                                                                               gaccGCATTAGCCGACCTTGGGTTCTTTTCTCTGTAGGCCTGATATGCCGCTTTCCGCGTCGCTTCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTg                                                                        <  2:1250630/147‑1 (MQ=255)
                                                                                    ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt                                                                   >  1:2467001/1‑150 (MQ=255)
                                                                                       tATCCGACATAGGGTTCGGTTCACTGTAGGCCTTATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg                                                                <  2:2461976/150‑1 (MQ=37)
                                                                                       caTCCGACATAGGGTTCTTTTCACTGTTGGCCTGAAAAGACGAATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg                                                                <  2:287889/149‑1 (MQ=37)
                                                                                         tCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg                                                              >  2:881176/1‑150 (MQ=255)
                                                                                          ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgc                                                             >  2:3953076/1‑150 (MQ=255)
                                                                                            gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgc                                                           >  1:1105924/1‑150 (MQ=255)
                                                                                             ccctcgggTTCTTTACCCTGTCGGCCTGATAACCCGCATAACGCGCCGCATACGCACACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg                                                          <  1:1060923/144‑1 (MQ=255)
                                                                                              ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcgc                                                         >  1:2544129/1‑150 (MQ=255)
                                                                                                tccgtttCTGTTCACTGTAGGCCTTAAAAGCCGCATTCCGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCAc                                                       <  2:1610802/145‑1 (MQ=255)
                                                                                                 aggtttCTTTTCAGTGTCGGCCTTATACGACGCATTCCGCGTCGCAACAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACg                                                      <  1:631920/146‑1 (MQ=255)
                                                                                                      ttttttCAGTGTTGGCCTTATATCCCGCTTTACGCGTCCAATCAGGCCACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTaa                                                 <  1:2665420/146‑1 (MQ=255)
                                                                                                         gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACa                                              >  2:500922/1‑150 (MQ=255)
                                                                                                          ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAg                                             >  1:2754886/1‑150 (MQ=255)
                                                                                                            cagtGTCGGCCTTATATGCCGCATTACGGGTCGAATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGAc                                           <  1:2787587/147‑1 (MQ=255)
                                                                                                               tGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgc                                        >  1:172158/1‑150 (MQ=255)
                                                                                                                gcccgCCTCAAAAGCCCCATCACGCGCCCCAACAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCCGCTGCTCAATCGACGCCAGAATCCCGCGCAGGATATACAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgcg                                       <  1:4348180/146‑1 (MQ=255)
                                                                                                                gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgcg                                       >  1:1526447/1‑150 (MQ=255)
                                                                                                                 tAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgcgc                                      >  2:3508693/1‑150 (MQ=255)
                                                                                                                 tAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgcgc                                      >  2:3550789/1‑150 (MQ=255)
                                                                                                                 tAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGgcgc                                      >  2:2983776/1‑150 (MQ=255)
                                                                                                                   ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCaa                                    >  1:2213988/1‑150 (MQ=255)
                                                                                                                    gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAAt                                   >  1:880597/1‑150 (MQ=255)
                                                                                                                      cttaTAAGACGAATTCCGCGTCGCATCTGGCCACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAAttt                                 <  1:3803053/147‑1 (MQ=255)
                                                                                                                      cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAAttt                                 >  2:3595966/1‑150 (MQ=255)
                                                                                                                       ttaTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTa                                <  2:2819670/148‑1 (MQ=255)
                                                                                                                        taTAAGACTCATTACGCGTCGAATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTAt                               <  2:1499042/149‑1 (MQ=255)
                                                                                                                        taTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTAt                               <  2:1500210/149‑1 (MQ=255)
                                                                                                                         aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTAtt                              >  1:3821577/1‑150 (MQ=255)
                                                                                                                           aatactcATTACGCGTCGCATCCGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATtca                            <  2:1405059/144‑1 (MQ=255)
                                                                                                                           aaGACGCATTACGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATtca                            <  2:916384/150‑1 (MQ=255)
                                                                                                                           aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATtca                            <  2:4257090/150‑1 (MQ=255)
                                                                                                                            cgccGCAGTCCGCGGCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATtcat                           <  2:2205808/147‑1 (MQ=255)
                                                                                                                            aGACGAATTACGCGTCGAATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATtcat                           <  2:2687306/150‑1 (MQ=255)
                                                                                                                                 cATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCt                      >  2:3161354/1‑150 (MQ=255)
                                                                                                                                   ttACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGc                    >  1:2362801/1‑150 (MQ=255)
                                                                                                                                      cgcgTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGcccc                 <  2:2074576/150‑1 (MQ=255)
                                                                                                                                      cgcgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGcccc                 <  2:2151148/150‑1 (MQ=255)
                                                                                                                                       gcgTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCg                <  2:3692138/150‑1 (MQ=255)
                                                                                                                                       gcgTCGCATCAGGCCACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCg                <  2:1981411/150‑1 (MQ=255)
                                                                                                                                        cgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCgg               >  1:346534/1‑150 (MQ=255)
                                                                                                                                          tCGCATCAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGAt             <  1:3030288/150‑1 (MQ=255)
                                                                                                                                           cGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATg            >  2:270752/1‑150 (MQ=255)
                                                                                                                                               tctggCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGttt        <  2:4090296/147‑1 (MQ=255)
                                                                                                                                               tCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGttt        <  2:3709297/150‑1 (MQ=255)
                                                                                                                                               tCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGttt        >  1:706490/1‑150 (MQ=255)
                                                                                                                                               gctggCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGttt        <  2:2802798/147‑1 (MQ=255)
                                                                                                                                                 aGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGTTTTc      >  1:327057/1‑150 (MQ=255)
                                                                                                                                                  ggCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGTTTTCa     >  2:2458687/1‑150 (MQ=255)
                                                                                                                                                  ggCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGTTTTCa     >  2:3807384/1‑150 (MQ=255)
                                                                                                                                                    cAACGGCTGTTATTCGGCCTTGTTACCGTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAAGTTATTCATCACCTGCCCCGGATGGTTTTCACg   >  2:3416743/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGACGGCGCAATTTATTCATCACCTGCCCCGGATGGTTTTCACGG  >  CP009273/2655785‑2656083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: