Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,610,619 C→A 5.2% 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,6190CA 5.2% 1498.5 / 31.6 374G389G (GGC→GGAbepAperiplasmic metalloprotease and chaperone for OM protein maintenance and assembly
Reads supporting (aligned to +/- strand):  ref base C (190/163);  new base A (9/11);  total (199/175)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

CTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCT  >  CP009273/2610471‑2610768
                                                                                                                                                    |                                                                                                                                                     
  ggCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGt                                                                                                                                                    >  1:2622233/1‑150 (MQ=255)
  ggCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGt                                                                                                                                                    <  2:1362479/150‑1 (MQ=255)
    cAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCa                                                                                                                                                  >  2:2602249/1‑150 (MQ=255)
     aaCGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCaa                                                                                                                                                 <  1:4391446/150‑1 (MQ=255)
     aaCGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCaa                                                                                                                                                 <  1:1722454/150‑1 (MQ=255)
     aaCGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCaa                                                                                                                                                 <  1:4283906/150‑1 (MQ=255)
       cGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAAcc                                                                                                                                               <  1:2005406/150‑1 (MQ=255)
         cATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACcac                                                                                                                                             >  2:1865690/1‑150 (MQ=255)
           tGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACaa                                                                                                                                           >  2:4340820/1‑150 (MQ=255)
               aTCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACaagaag                                                                                                                                       <  1:4693419/150‑1 (MQ=255)
                 ctcGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCg                                                                                                                                     <  2:143239/150‑1 (MQ=255)
                   cGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGc                                                                                                                                   <  2:3860785/150‑1 (MQ=255)
                       cTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAat                                                                                                                               <  2:5088771/150‑1 (MQ=255)
                                 atatCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATc                                                                                                                     >  1:5039333/1‑150 (MQ=255)
                                  tatCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCg                                                                                                                    <  1:56003/150‑1 (MQ=255)
                                  tatCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCg                                                                                                                    <  1:1364074/150‑1 (MQ=255)
                                      gATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATAtt                                                                                                                            <  2:2313285/138‑1 (MQ=255)
                                           tGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCtt                                                                                                           >  1:862970/1‑150 (MQ=255)
                                            gggCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCttt                                                                                                          >  2:2500716/1‑150 (MQ=255)
                                             ggCAAAACAAAGCCAATGAGGCGATCAATCGTCTGTAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCATCCTGGCTATCGCTTATTTGCAAGGCGGTAAACCACAAGAATCGGCCAATATTATGAATCTCTACAACTTta                                                                                                         <  2:4748917/150‑1 (MQ=255)
                                               cAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaat                                                                                                       <  1:4012629/150‑1 (MQ=255)
                                               cAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaat                                                                                                       <  1:2579102/150‑1 (MQ=255)
                                               cAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaat                                                                                                       >  2:315028/1‑150 (MQ=255)
                                                aaaaCAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaata                                                                                                      >  2:3206477/1‑150 (MQ=255)
                                                    cAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTtaataataa                                                                                                  >  1:629550/1‑150 (MQ=255)
                                                     aaaGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATaaa                                                                                                 <  2:1620926/150‑1 (MQ=255)
                                                     aaaGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTATAATAATaaa                                                                                                 >  2:3698957/1‑150 (MQ=255)
                                                       aGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGa                                                                                               >  1:1158526/1‑150 (MQ=255)
                                                       aGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGa                                                                                               <  1:2800459/150‑1 (MQ=255)
                                                         ccAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATg                                                                                             <  1:1351673/150‑1 (MQ=255)
                                                          cAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGa                                                                                            >  2:833062/1‑150 (MQ=255)
                                                          cAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGa                                                                                            >  1:4366184/1‑150 (MQ=255)
                                                            aTGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACa                                                                                          <  2:5114101/150‑1 (MQ=255)
                                                              gAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGc                                                                                        <  2:1746399/150‑1 (MQ=255)
                                                               aGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCa                                                                                       <  2:4300636/150‑1 (MQ=255)
                                                                ggCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCaa                                                                                      <  2:3626951/150‑1 (MQ=255)
                                                                 gCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAAc                                                                                     <  1:477603/150‑1 (MQ=255)
                                                                 gCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAAc                                                                                     <  1:4469163/150‑1 (MQ=255)
                                                                  cGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACg                                                                                    >  2:628650/1‑150 (MQ=255)
                                                                   gATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGACAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTGTAAGCATAAAGATGACAGCCACgg                                                                                   >  2:4832393/1‑150 (MQ=255)
                                                                   gATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACgg                                                                                   <  1:2134466/150‑1 (MQ=255)
                                                                    aTCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGc                                                                                  >  1:2543974/1‑150 (MQ=255)
                                                                       aaTCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTgg                                                                               >  1:2258447/1‑150 (MQ=255)
                                                                          cGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGTGGTAAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGAt                                                                            <  1:3477105/150‑1 (MQ=255)
                                                                           gTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGAtt                                                                           >  1:4558839/1‑150 (MQ=255)
                                                                           gTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGAtt                                                                           <  1:945621/150‑1 (MQ=255)
                                                                            tCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGAttt                                                                          <  2:2689420/150‑1 (MQ=255)
                                                                              tGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGc                                                                        <  1:2472443/150‑1 (MQ=255)
                                                                                aaaaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTg                                                                      >  2:3760073/1‑150 (MQ=255)
                                                                                 aaaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTgg                                                                     >  1:785812/1‑150 (MQ=255)
                                                                                    aTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGcac                                                                  <  2:5050777/150‑1 (MQ=255)
                                                                                    aTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGcac                                                                  >  2:3552526/1‑150 (MQ=255)
                                                                                      gCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAg                                                                <  1:3683194/150‑1 (MQ=255)
                                                                                        ccGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGc                                                              <  2:1923777/150‑1 (MQ=255)
                                                                                         cgcgATTTACGCACCAAGCCGGTGTTGCAGCTCAACCTGGCGAACGCTTCTTTGCACGGCGGTCAACCACAAGAAGCGGCCAAGATTCTGAATCGCTACACCGTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCg                                                             <  2:2741874/150‑1 (MQ=255)
                                                                                          gcgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCgg                                                            >  2:4053673/1‑150 (MQ=255)
                                                                                            gATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGaa                                                          >  1:1019245/1‑150 (MQ=255)
                                                                                            gATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGaa                                                          >  1:4159646/1‑150 (MQ=255)
                                                                                             aTTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGGAg                                                         >  2:2347790/1‑150 (MQ=255)
                                                                                             aTTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAg                                                         <  1:5092139/150‑1 (MQ=255)
                                                                                             aTTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAg                                                         >  2:2868867/1‑150 (MQ=255)
                                                                                             aTTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAg                                                         >  2:4145846/1‑150 (MQ=255)
                                                                                               ttACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGcc                                                       >  1:4335971/1‑150 (MQ=255)
                                                                                               ttACGCACCAATACGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGTATATGCTGGCACAGGCGGAAGcc                                                       >  1:5005383/1‑150 (MQ=255)
                                                                                                  cGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgcg                                                    <  1:1612679/150‑1 (MQ=255)
                                                                                                      ccAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTaa                                                <  1:4073088/150‑1 (MQ=255)
                                                                                                      ccAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTaa                                                >  2:442251/1‑150 (MQ=255)
                                                                                                      ccAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTaa                                                >  2:4940344/1‑150 (MQ=255)
                                                                                                        aaTCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAAt                                              <  2:326002/150‑1 (MQ=255)
                                                                                                            cGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAAcc                                          <  2:1163487/150‑1 (MQ=255)
                                                                                                            cGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGAGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACCGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAAcc                                          >  2:2206610/1‑150 (MQ=255)
                                                                                                             ggTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcg                                         >  2:3946420/1‑150 (MQ=255)
                                                                                                              gtgtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcgc                                        <  1:1457712/150‑1 (MQ=255)
                                                                                                              gtgtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcgc                                        >  2:100649/1‑150 (MQ=255)
                                                                                                               tgtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcgcg                                       >  1:1094232/1‑150 (MQ=255)
                                                                                                               tgtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcgcg                                       >  1:1538581/1‑150 (MQ=255)
                                                                                                                gtTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGa                                      >  2:1773903/1‑150 (MQ=255)
                                                                                                                  tGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATc                                    <  2:4609235/150‑1 (MQ=255)
                                                                                                                  tGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATc                                    >  2:1661502/1‑150 (MQ=255)
                                                                                                                     aGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAgg                                 >  2:768184/1‑150 (MQ=255)
                                                                                                                       cTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAg                               <  1:3633243/150‑1 (MQ=255)
                                                                                                                       cTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAg                               <  1:3326326/150‑1 (MQ=255)
                                                                                                                        tCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgc                              >  1:1424208/1‑150 (MQ=255)
                                                                                                                         cAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgct                             >  1:2954625/1‑150 (MQ=255)
                                                                                                                           aCCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctgg                           <  2:3040193/150‑1 (MQ=255)
                                                                                                                            ccTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggc                          <  1:2288368/150‑1 (MQ=255)
                                                                                                                             cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                         <  1:3777220/150‑1 (MQ=255)
                                                                                                                             cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                         <  1:1205378/150‑1 (MQ=255)
                                                                                                                             cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                         >  2:3297230/1‑150 (MQ=255)
                                                                                                                                  gAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGcgcg                    <  1:3083025/150‑1 (MQ=255)
                                                                                                                                      gCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCa                <  1:2878438/150‑1 (MQ=255)
                                                                                                                                       cTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAg               >  1:2245212/1‑150 (MQ=255)
                                                                                                                                       cTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAg               >  1:4844862/1‑150 (MQ=255)
                                                                                                                                         tATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGaa             >  2:2099216/1‑150 (MQ=255)
                                                                                                                                          aTTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAg            <  1:2358044/150‑1 (MQ=255)
                                                                                                                                          aTTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAg            <  1:2748277/150‑1 (MQ=255)
                                                                                                                                               cAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTAt       >  1:3244860/1‑150 (MQ=255)
                                                                                                                                               cAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTAt       >  2:4519138/1‑150 (MQ=255)
                                                                                                                                               cAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTAt       >  1:3953250/1‑150 (MQ=255)
                                                                                                                                                aaGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATg      <  2:3424943/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCAACGCATGGTATCTCGATCTGGCTACTGATATCGATCTTGGGCAAAACAAAGCCAATGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCT  >  CP009273/2610471‑2610768

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: