Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,610,619 C→A 12.6% G389G (GGC→GGA bepA → periplasmic metalloprotease and chaperone for OM protein maintenance and assembly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,610,6190CA12.6% 970.9 / 77.1 291G389G (GGC→GGAbepAperiplasmic metalloprotease and chaperone for OM protein maintenance and assembly
Reads supporting (aligned to +/- strand):  ref base C (141/113);  new base A (20/17);  total (161/130)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 291 total aligned reads displayed.

CTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGC  >  CP009273/2610471‑2610765
                                                                                                                                                    |                                                                                                                                                  
cTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCACTCCGGTGTTGCAGCTCAACCTGGCGAACGCTGATTTGCAAGGCg                                                                                                                                                   >  2:4329028/1‑150 (MQ=255)
cTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCg                                                                                                                                                   <  1:412062/150‑1 (MQ=255)
  ggCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGt                                                                                                                                                 <  1:1423338/150‑1 (MQ=255)
   gCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTc                                                                                                                                                <  1:3996455/150‑1 (MQ=255)
    cAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCa                                                                                                                                               <  1:2635136/150‑1 (MQ=255)
    cAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCa                                                                                                                                               >  1:1920852/1‑150 (MQ=255)
      aCGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAAc                                                                                                                                             <  1:1202141/150‑1 (MQ=255)
      aCGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAAc                                                                                                                                             <  2:2622997/150‑1 (MQ=255)
       cGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTACAAGGCGGTCAAcc                                                                                                                                            <  2:1563539/150‑1 (MQ=255)
          aTGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACcaca                                                                                                                                         >  2:1273192/1‑150 (MQ=255)
            ggTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACaag                                                                                                                                       <  1:2526228/150‑1 (MQ=255)
             gTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACaaga                                                                                                                                      <  1:3975322/150‑1 (MQ=255)
              tATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACaagaa                                                                                                                                     >  2:21767/1‑150 (MQ=255)
                 ctcGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCg                                                                                                                                  <  2:540119/150‑1 (MQ=255)
                   cGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGc                                                                                                                                <  1:2070504/150‑1 (MQ=255)
                     aTCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCa                                                                                                                              <  1:2804529/150‑1 (MQ=255)
                          gCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATAtt                                                                                                                         <  1:4070643/150‑1 (MQ=255)
                             aCTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTg                                                                                                                      >  1:707763/1‑150 (MQ=255)
                                  tatCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCg                                                                                                                 <  2:3167764/150‑1 (MQ=255)
                                     cGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCCCCCCCGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTa                                                                                                              <  1:1064877/150‑1 (MQ=255)
                                     cGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTa                                                                                                              >  2:3578250/1‑150 (MQ=255)
                                      gATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTac                                                                                                             >  2:3323749/1‑150 (MQ=255)
                                       aTCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTaca                                                                                                            >  1:1328267/1‑150 (MQ=255)
                                       aTCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTaca                                                                                                            >  1:38457/1‑150 (MQ=255)
                                        tCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTacac                                                                                                           >  1:4254224/1‑150 (MQ=255)
                                          ttGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCt                                                                                                         <  2:3533490/150‑1 (MQ=255)
                                          ttGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCt                                                                                                         <  2:1703548/150‑1 (MQ=255)
                                          ttGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCt                                                                                                         <  1:2462061/150‑1 (MQ=255)
                                            gggCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCttt                                                                                                       <  1:1456985/150‑1 (MQ=255)
                                                aaaaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaata                                                                                                   >  2:3502000/1‑150 (MQ=255)
                                                aaaaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaata                                                                                                   >  1:1299384/1‑150 (MQ=255)
                                                 aaaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataa                                                                                                  >  2:4534366/1‑150 (MQ=255)
                                                  aaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataat                                                                                                 >  1:2972290/1‑150 (MQ=255)
                                                  aaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataat                                                                                                 >  2:771560/1‑150 (MQ=255)
                                                  aaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataat                                                                                                 >  2:771667/1‑150 (MQ=255)
                                                  aaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataat                                                                                                 <  1:1149983/150‑1 (MQ=255)
                                                   aCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataata                                                                                                >  1:2327260/1‑150 (MQ=255)
                                                   aCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataata                                                                                                <  2:1330039/150‑1 (MQ=255)
                                                     aaaGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATaaa                                                                                              <  1:3492409/150‑1 (MQ=255)
                                                      aaGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAg                                                                                             >  1:1420912/1‑150 (MQ=255)
                                                         ccAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATg                                                                                          >  1:1146010/1‑150 (MQ=255)
                                                            aTGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACa                                                                                       >  2:2710281/1‑150 (MQ=255)
                                                            aTGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACa                                                                                       <  1:871000/150‑1 (MQ=255)
                                                             tGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAg                                                                                      >  2:3803289/1‑150 (MQ=255)
                                                              gAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGc                                                                                     >  1:3937686/1‑150 (MQ=255)
                                                                ggCGATCTATCGTATGAAAAATGCCCTCTATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACTCCTTTAATCATAAAGATGACAGCaa                                                                                   <  1:1401780/150‑1 (MQ=255)
                                                                ggCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCaa                                                                                   <  1:883111/150‑1 (MQ=255)
                                                                ggCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCaa                                                                                   <  1:950603/150‑1 (MQ=255)
                                                                   gATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACgg                                                                                >  1:1214553/1‑150 (MQ=255)
                                                                   gATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACgg                                                                                <  1:1518198/150‑1 (MQ=255)
                                                                    aTCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGc                                                                               >  1:372787/1‑150 (MQ=255)
                                                                    aTCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGc                                                                               <  1:3088706/150‑1 (MQ=255)
                                                                      cAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTg                                                                             >  1:3859111/1‑150 (MQ=255)
                                                                         tCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATCTTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGa                                                                          <  2:791655/150‑1 (MQ=255)
                                                                         tCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGa                                                                          <  2:791610/150‑1 (MQ=255)
                                                                              tGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGc                                                                     >  2:1097310/1‑150 (MQ=255)
                                                                                aaaaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTg                                                                   >  2:2941151/1‑150 (MQ=255)
                                                                                 aaaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTgg                                                                  <  1:2819707/150‑1 (MQ=255)
                                                                                  aaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGc                                                                 >  2:2249282/1‑150 (MQ=255)
                                                                                  aaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGc                                                                 >  1:1991093/1‑150 (MQ=255)
                                                                                     tGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGcaca                                                              <  1:2733001/150‑1 (MQ=255)
                                                                                      gCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAg                                                             <  1:3557192/150‑1 (MQ=255)
                                                                                        ccGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCCATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGc                                                           >  1:3628374/1‑150 (MQ=255)
                                                                                        ccGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGc                                                           >  1:3628379/1‑150 (MQ=255)
                                                                                         cgcgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCg                                                          >  2:2057539/1‑150 (MQ=255)
                                                                                         cgcgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCg                                                          >  2:4421293/1‑150 (MQ=255)
                                                                                          gcgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCgg                                                         >  2:4508361/1‑150 (MQ=255)
                                                                                           cgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGa                                                        >  1:987824/1‑150 (MQ=255)
                                                                                            gATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGaa                                                       >  2:1977015/1‑150 (MQ=255)
                                                                                                tACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcg                                                   <  2:2589983/150‑1 (MQ=255)
                                                                                                 aCGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgc                                                  >  2:814028/1‑150 (MQ=255)
                                                                                                  cGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgcg                                                 <  1:3323749/150‑1 (MQ=255)
                                                                                                  cGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgcg                                                 <  2:1838455/150‑1 (MQ=255)
                                                                                                  cGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgcg                                                 >  2:1867228/1‑150 (MQ=255)
                                                                                                    cACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCt                                               <  1:251033/150‑1 (MQ=255)
                                                                                                     aCCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTa                                              >  1:2056695/1‑150 (MQ=255)
                                                                                                        aaTCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAAt                                           >  2:2463525/1‑150 (MQ=255)
                                                                                                         aTCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATa                                          >  2:3834148/1‑150 (MQ=255)
                                                                                                         aTCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATa                                          >  2:3834395/1‑150 (MQ=255)
                                                                                                             ggTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcg                                      <  2:585099/150‑1 (MQ=255)
                                                                                                              gtgtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcgc                                     >  1:774310/1‑150 (MQ=255)
                                                                                                                gtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGa                                   >  1:3165583/1‑150 (MQ=255)
                                                                                                                 ttGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGAt                                  <  1:748000/150‑1 (MQ=255)
                                                                                                                  tgaaGCTCAACCTGGCGAACGCTTATTTGCCAGGCGGTCAACCACAAGAAGCGGCCTATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATc                                 <  2:3024068/147‑1 (MQ=255)
                                                                                                                  tGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATc                                 >  1:834993/1‑150 (MQ=255)
                                                                                                                   gCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCa                                >  1:1187763/1‑150 (MQ=255)
                                                                                                                       cTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAg                            <  2:3475462/150‑1 (MQ=255)
                                                                                                                        tCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgc                           <  1:4532548/150‑1 (MQ=255)
                                                                                                                          aaCCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctg                         >  2:3593976/1‑150 (MQ=255)
                                                                                                                           aCCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctgg                        <  2:4170245/150‑1 (MQ=255)
                                                                                                                             cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                      >  2:4086006/1‑150 (MQ=255)
                                                                                                                             cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                      >  2:195839/1‑150 (MQ=255)
                                                                                                                               ggCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTgc                    >  2:1118707/1‑150 (MQ=255)
                                                                                                                                gCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTgcg                   >  2:3737027/1‑150 (MQ=255)
                                                                                                                                  gAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGcgcg                 >  1:3472210/1‑150 (MQ=255)
                                                                                                                                      gCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCa             >  2:2416209/1‑150 (MQ=255)
                                                                                                                                       cTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAg            >  1:1053606/1‑150 (MQ=255)
                                                                                                                                                aaGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATg   >  1:587860/1‑150 (MQ=255)
                                                                                                                                                 aGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATgc  >  2:600748/1‑150 (MQ=255)
                                                                                                                                                 aGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAAGAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATgc  >  2:242459/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGC  >  CP009273/2610471‑2610765

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: