Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,009,272 G→A 6.2% intergenic (+175/‑147) mocA → / → ygfK CTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,009,2720GA 6.2% 1166.9 / 47.2 320intergenic (+175/‑147)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (173/126);  new base A (15/5);  total (189/131)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 320 total aligned reads displayed.

AATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAGGGTTAT  >  CP009273/3009124‑3009420
                                                                                                                                                    |                                                                                                                                                    
aaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGc                                                                                                                                                     >  2:300517/1‑150 (MQ=255)
   atTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGTCCt                                                                                                                                                  >  1:3088337/1‑150 (MQ=255)
     ttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTa                                                                                                                                                <  1:1399133/150‑1 (MQ=255)
      tAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTaa                                                                                                                                               >  2:3960967/1‑150 (MQ=255)
          ccgtCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                                                                                                           <  2:287771/147‑1 (MQ=255)
           cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                                                                                                          <  2:1027592/150‑1 (MQ=255)
            ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCg                                                                                                                                         <  2:967942/150‑1 (MQ=255)
            ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCg                                                                                                                                         <  2:3813667/150‑1 (MQ=255)
              ccGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAg                                                                                                                                       <  2:1654545/150‑1 (MQ=255)
                tgctcTCATTATCATAATGTTAAGGTTTCAATGATTTTGAGAAAAGATAAAAATTGATAATTGGGCTGGCAATAATATTTTTCAGTAGTGATTTTAATTCTCATTTAGCACCTGTGATTTCATTCACAATAAGCCCTTAAGTTTCGAgcc                                                                                                                                     <  1:3385296/146‑1 (MQ=255)
                     tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                                                                                                                <  1:2072057/150‑1 (MQ=255)
                          aTCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                                                                                                           <  1:2292561/150‑1 (MQ=255)
                           tCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttt                                                                                                                          <  1:394617/150‑1 (MQ=255)
                              tAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttttt                                                                                                                       <  1:853620/150‑1 (MQ=255)
                                 tGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                                                                                                                    <  2:3236476/150‑1 (MQ=255)
                                  gATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGATCCGCCTCACATTTTTTtata                                                                                                                   >  2:4582392/1‑150 (MQ=255)
                                   aTAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTAAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAATCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTatat                                                                                                                  >  2:1548438/1‑150 (MQ=255)
                                    tAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAtt                                                                                                                 <  2:3685437/150‑1 (MQ=255)
                                     aaGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                                                                                                                >  1:400233/1‑150 (MQ=255)
                                     aaGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                                                                                                                <  2:505883/150‑1 (MQ=255)
                                        gAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTccc                                                                                                             >  1:4476614/1‑150 (MQ=255)
                                        gAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTccc                                                                                                             <  1:3591771/150‑1 (MQ=255)
                                            cACTGATTTTGAGAAAACAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCa                                                                                                         <  1:4600804/150‑1 (MQ=255)
                                              cTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa                                                                                                       >  2:3175113/1‑150 (MQ=255)
                                              cTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa                                                                                                       >  2:2565041/1‑150 (MQ=255)
                                              cTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTGTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa                                                                                                       >  2:4442619/1‑150 (MQ=255)
                                                 aTTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCt                                                                                                    >  1:1182864/1‑150 (MQ=255)
                                                  ttttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTg                                                                                                   >  1:4671422/1‑150 (MQ=255)
                                                   tttGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTgg                                                                                                  >  1:1744549/1‑150 (MQ=255)
                                                    ttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGc                                                                                                 >  2:2892482/1‑150 (MQ=255)
                                                    ttGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGc                                                                                                 >  1:1641266/1‑150 (MQ=255)
                                                      gagaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCac                                                                                               >  2:1421052/1‑149 (MQ=255)
                                                         aaaagataaaaAATGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGATCATTCACAAAAAGCCCTTAAGTTTCGATCCGCCTCACATTATTTTATATTTCCCGCCAAACCTGGCAagag                                                                                            >  2:3549911/1‑150 (MQ=255)
                                                              ataaaaaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtggtg                                                                                       >  1:131673/3‑150 (MQ=255)
                                                               gattatttGTTAATTGGGCTGGCAATATTAATTTTCAGGTGTGTTTTTATTTCTCATTTTGCTCCTGTGAGTTCATTCACTAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTGCCGCCAAACCTGGCAAGTGTGGTGc                                                                                      <  2:4760827/144‑1 (MQ=255)
                                                                aaaaaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCg                                                                                     <  2:121181/150‑1 (MQ=255)
                                                                aaaaaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCg                                                                                     >  1:4535268/1‑150 (MQ=255)
                                                                aaaaaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCg                                                                                     >  2:3398459/1‑150 (MQ=255)
                                                                   aaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGAtt                                                                                  >  2:2644624/1‑150 (MQ=255)
                                                                   aaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGAtt                                                                                  >  2:4627173/1‑150 (MQ=255)
                                                                   aaTTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGAGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGAGGTGCGAtt                                                                                  >  2:4627150/1‑150 (MQ=255)
                                                                       gATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGAttgttg                                                                              >  1:972663/1‑150 (MQ=255)
                                                                        aTAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAACCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGc                                                                             >  2:870130/1‑150 (MQ=255)
                                                                         tAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGct                                                                            >  2:3482944/1‑150 (MQ=255)
                                                                          aaaTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGctc                                                                           >  1:3683913/1‑150 (MQ=255)
                                                                           aaTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGctct                                                                          >  2:608958/1‑150 (MQ=255)
                                                                           aaTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGctct                                                                          >  1:1983392/1‑150 (MQ=255)
                                                                           aaTGGGCTGGCAAGAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGctct                                                                          >  2:1082398/1‑150 (MQ=255)
                                                                            aTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTa                                                                         >  1:642551/1‑150 (MQ=255)
                                                                            aTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTa                                                                         >  1:4469583/1‑150 (MQ=255)
                                                                               ggCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATcc                                                                      >  1:3975132/1‑150 (MQ=255)
                                                                                   ggCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCt                                                                  <  1:3533430/150‑1 (MQ=255)
                                                                                    gCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTa                                                                 >  1:2134170/1‑150 (MQ=255)
                                                                                    gCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTa                                                                 >  2:4152442/1‑150 (MQ=255)
                                                                                    gCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTa                                                                 >  1:1865189/1‑150 (MQ=255)
                                                                                     cAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTaa                                                                >  1:5123850/1‑150 (MQ=255)
                                                                                     cAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTaa                                                                <  2:1575947/150‑1 (MQ=255)
                                                                                      aataataatTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTaaa                                                               <  1:1805537/150‑1 (MQ=255)
                                                                                       ataataatTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAac                                                              >  1:3072181/1‑150 (MQ=255)
                                                                                          ataatTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAaccac                                                           <  1:1508203/150‑1 (MQ=255)
                                                                                          ataatTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAaccac                                                           >  2:81649/1‑150 (MQ=255)
                                                                                           taatTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAaccacc                                                          >  2:3099589/1‑150 (MQ=255)
                                                                                               tttCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTTCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGAt                                                      >  1:1002448/1‑150 (MQ=255)
                                                                                               tttCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGAt                                                      <  1:4665421/150‑1 (MQ=255)
                                                                                                 tCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGAttt                                                    <  1:1591327/150‑1 (MQ=255)
                                                                                                   aGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTtct                                                  <  2:4476284/150‑1 (MQ=255)
                                                                                                    ggAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTtctc                                                 >  2:769434/1‑150 (MQ=255)
                                                                                                     gAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCa                                                <  1:928373/150‑1 (MQ=255)
                                                                                                      aGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCTCATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCaa                                               <  2:4886124/150‑1 (MQ=255)
                                                                                                        tGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGATCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAaca                                             >  2:5064322/1‑150 (MQ=255)
                                                                                                        tGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAaca                                             >  1:3277473/1‑150 (MQ=255)
                                                                                                        tGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAaca                                             <  2:2329554/150‑1 (MQ=255)
                                                                                                         gATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAacac                                            >  1:2806588/1‑150 (MQ=255)
                                                                                                          aTTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACAcc                                           >  2:2216689/1‑150 (MQ=255)
                                                                                                          aTTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACAcc                                           >  2:4824442/1‑150 (MQ=255)
                                                                                                           ttttAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCg                                          >  1:753818/1‑150 (MQ=255)
                                                                                                              tAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAACCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTc                                       >  1:345701/1‑150 (MQ=255)
                                                                                                                  tctcATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTc                                   <  2:3653145/150‑1 (MQ=255)
                                                                                                                  tctcATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTc                                   >  1:2370392/1‑150 (MQ=255)
                                                                                                                      aTTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATg                               >  2:1994970/1‑150 (MQ=255)
                                                                                                                          aGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGATCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGatat                           >  2:2410647/1‑150 (MQ=255)
                                                                                                                           gCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATc                          <  1:2565041/150‑1 (MQ=255)
                                                                                                                           gCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATc                          >  1:2926619/1‑150 (MQ=255)
                                                                                                                             aCCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTg                        <  2:4243362/150‑1 (MQ=255)
                                                                                                                              ccTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGt                       <  1:4893253/150‑1 (MQ=255)
                                                                                                                              ccTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGt                       <  1:41615/150‑1 (MQ=255)
                                                                                                                                 gtgAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGtata                    <  1:1427173/150‑1 (MQ=255)
                                                                                                                                 gtgAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGtata                    <  1:4823579/150‑1 (MQ=255)
                                                                                                                                   gAGTTCCTTCACAACAAGCCCGTAAGGTGCGAGCCGCCTCACATTTTTTTATATGTCCCGCCAACCCTGGCAAGAGTGGTGCGATTGTTGCTCTAGCCCCCTACACCCCCGGAGTTCTCCACACCGGTCACTCACTGAGATCTGtctacg                  <  2:3882267/150‑3 (MQ=255)
                                                                                                                                    aGTTCATTCACCAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGc                 >  1:5045735/1‑150 (MQ=255)
                                                                                                                                    aGTTCATTCACAAAAAACCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGc                 <  2:27614/150‑1 (MQ=255)
                                                                                                                                     gTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCt                <  1:2670878/150‑1 (MQ=255)
                                                                                                                                     gTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCt                >  2:2170935/1‑150 (MQ=255)
                                                                                                                                       tcattcaCAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTaa              <  2:1285295/150‑1 (MQ=255)
                                                                                                                                        cattcaCAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAg             >  2:1749415/1‑150 (MQ=255)
                                                                                                                                           tcaCAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGgag          >  2:3266228/1‑150 (MQ=255)
                                                                                                                                              caAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAgg       >  1:1760946/1‑150 (MQ=255)
                                                                                                                                               aaaaaGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAggg      >  2:236443/1‑150 (MQ=255)
                                                                                                                                                aaaaGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAGGGt     >  1:3546770/1‑150 (MQ=255)
                                                                                                                                                 aaaGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAGGGtt    >  1:3233167/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGCGATTGTTGCTCTATCCCCCTAAACCACCGGATTTCTCAACACCGGTCACTCAATGATATCTGTATAAGCTAAGGAGAGGGTTAT  >  CP009273/3009124‑3009420

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