Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 291,348 G→A 93.1% intergenic (‑58/+59) yagN ← / ← intF CP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273291,3480GA93.1% 686.5 / 33.5 216intergenic (‑58/+59)yagN/intFCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base G (10/5);  new base A (69/132);  total (79/137)
Fisher's exact test for biased strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 224 total aligned reads displayed.

ATTGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAAC  >  CP009273/291203‑291491
                                                                                                                                                 |                                                                                                                                               
aTTGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATaaa                                                                                                                                               >  2:116280/1‑148 (MQ=255)
aTTGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATaaa                                                                                                                                               <  1:1238794/148‑1 (MQ=255)
aTTGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATaaa                                                                                                                                               <  1:2328258/148‑1 (MQ=255)
 ttGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATaaaa                                                                                                                                              >  2:1316746/1‑148 (MQ=255)
   gAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAat                                                                                                                                            <  2:839217/148‑1 (MQ=255)
   gAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAat                                                                                                                                            >  1:1889580/1‑148 (MQ=255)
     gttgttTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAatct                                                                                                                                          <  2:322077/148‑1 (MQ=255)
          tttACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCtt                                                                                                                                     <  1:1116099/148‑1 (MQ=255)
          tttACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCtt                                                                                                                                     <  1:1439358/148‑1 (MQ=255)
           ttACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTg                                                                                                                                    <  2:1860936/148‑1 (MQ=255)
            tACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGa                                                                                                                                   <  2:1078086/148‑1 (MQ=255)
             aCTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGAt                                                                                                                                  >  2:425206/1‑148 (MQ=255)
               ttGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCa                                                                                                                                <  1:518747/148‑1 (MQ=255)
                tGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCaa                                                                                                                               >  2:1989415/1‑148 (MQ=255)
                   aTTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAg                                                                                                                            <  1:741026/148‑1 (MQ=255)
                   aTTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAg                                                                                                                            <  2:561673/148‑1 (MQ=255)
                      aaTACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtg                                                                                                                         <  1:3194234/148‑1 (MQ=255)
                       aTACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtga                                                                                                                        >  2:2632476/1‑148 (MQ=255)
                       aTACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtga                                                                                                                        >  1:938204/1‑148 (MQ=255)
                        tACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtgat                                                                                                                       >  1:531122/1‑148 (MQ=255)
                        tACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtgat                                                                                                                       <  1:3040109/148‑1 (MQ=255)
                          cacCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtgatga                                                                                                                     >  1:866753/1‑148 (MQ=255)
                          cacCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTtgatga                                                                                                                     >  2:1542125/1‑148 (MQ=255)
                            ccAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGGTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGata                                                                                                                   <  2:638014/148‑1 (MQ=255)
                             cAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGatat                                                                                                                  <  1:3142509/148‑1 (MQ=255)
                             cAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGatat                                                                                                                  <  2:1443712/148‑1 (MQ=255)
                              aGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATAtt                                                                                                                 <  1:407862/148‑1 (MQ=255)
                              aGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATAtt                                                                                                                 <  1:836896/148‑1 (MQ=255)
                               gCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCGTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATAttt                                                                                                                >  1:2340906/1‑148 (MQ=255)
                                  aGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGt                                                                                                             <  2:903258/148‑1 (MQ=255)
                                   gCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTa                                                                                                            <  2:1240754/148‑1 (MQ=255)
                                    ctctGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTaa                                                                                                           >  1:213727/1‑148 (MQ=255)
                                      ctGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAata                                                                                                         <  1:1143291/148‑1 (MQ=255)
                                       tGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTTGTTGCTGTATCATTGACTAAACTGTGTATTAAGATACATAAAAATCTATCTTGATCAATAGTTGTTGATATTTTGTAatat                                                                                                        >  1:1523223/1‑148 (MQ=255)
                                        gATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATg                                                                                                       <  1:2188953/148‑1 (MQ=255)
                                        gATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATg                                                                                                       <  2:3476182/148‑1 (MQ=255)
                                         aTAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGt                                                                                                      <  1:325288/148‑1 (MQ=255)
                                           agatGCCAGAGTGTGTTCTATTGAAACAGATGCTGGGGTTGCCAATTGCTTATCTCCTTTGAAGTTGCTGTATCATTGACTAAACTGTGTATTAAGATACAAAAAAAACTATCAAGATCAATATATTAATTTATTTTTTAATATGTTg                                                                                                    >  2:1504525/4‑148 (MQ=255)
                                            gCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGa                                                                                                   >  2:1736347/1‑148 (MQ=255)
                                             cTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGAt                                                                                                  <  2:1435836/148‑1 (MQ=255)
                                                ccAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGAtttt                                                                                               <  2:531122/148‑1 (MQ=255)
                                                  agagTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTa                                                                                             <  2:2634794/148‑1 (MQ=255)
                                                  agagTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTa                                                                                             <  1:1113382/148‑1 (MQ=255)
                                                   gagTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTaa                                                                                            >  2:2317661/1‑148 (MQ=255)
                                                      tGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATgg                                                                                         >  2:2720741/1‑148 (MQ=255)
                                                        ggTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGaa                                                                                       <  2:1121915/148‑1 (MQ=255)
                                                        ggTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGaa                                                                                       <  2:2895715/148‑1 (MQ=255)
                                                           tCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCa                                                                                    <  1:2946000/148‑1 (MQ=255)
                                                              aTGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACg                                                                                 <  2:1178618/148‑1 (MQ=255)
                                                                ggATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCa                                                                               <  1:186905/148‑1 (MQ=255)
                                                                ggATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCa                                                                               <  1:187015/148‑1 (MQ=255)
                                                                      aGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCaa                                                                         <  1:2225683/148‑1 (MQ=255)
                                                                      aGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCaa                                                                         <  1:2317661/148‑1 (MQ=255)
                                                                        ttGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACg                                                                       >  1:2594115/1‑148 (MQ=255)
                                                                           cTGGGGTTGCCATGGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACgttg                                                                    <  1:1989534/148‑1 (MQ=255)
                                                                            tGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATAATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACgttgt                                                                   <  2:2109920/148‑1 (MQ=255)
                                                                              gggTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGttt                                                                 <  2:686039/148‑1 (MQ=255)
                                                                              gggTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGttt                                                                 <  2:2292630/148‑1 (MQ=255)
                                                                              gggTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGttt                                                                 >  2:1681664/1‑148 (MQ=255)
                                                                                 ttGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGtttttt                                                              <  1:1766857/148‑1 (MQ=255)
                                                                                  tGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTa                                                             <  2:213727/148‑1 (MQ=255)
                                                                                  tGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTa                                                             <  2:2010832/148‑1 (MQ=255)
                                                                                    ccATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATg                                                           <  2:2385359/148‑1 (MQ=255)
                                                                                    ccATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATg                                                           >  2:296966/1‑148 (MQ=255)
                                                                                     cATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATgg                                                          <  1:1060323/148‑1 (MQ=255)
                                                                                      aTTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGt                                                         >  1:3213257/1‑148 (MQ=255)
                                                                                       ttGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTa                                                        <  2:1889580/148‑1 (MQ=255)
                                                                                            ttATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATAcc                                                   <  1:3112478/148‑1 (MQ=255)
                                                                                             tATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGAGTATTACGATACATAAAAATCTATCTTGATAAATAGTTGATGATAATATGTAATATGTTGATTTTTAATGGAATCAACGCATAGCAACGTTGTTTTTTATGGTAATACCt                                                  >  2:357495/1‑148 (MQ=255)
                                                                                              aTCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTg                                                 <  1:3446961/148‑1 (MQ=255)
                                                                                              aTCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTg                                                 <  2:3557263/148‑1 (MQ=255)
                                                                                              aTCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATAGTTTGTAATATGTTGATTTGTAATGGAATCAACGCATATCAACGTGGTTTTTTATGGTAATACCTg                                                 >  2:32217/1‑148 (MQ=255)
                                                                                               tctcCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGc                                                >  2:1943127/1‑148 (MQ=255)
                                                                                                ctcCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCt                                               >  2:244241/1‑148 (MQ=255)
                                                                                                    tgtgTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCAt                                           <  1:605866/148‑1 (MQ=255)
                                                                                                       gtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                        <  2:3338450/148‑1 (MQ=255)
                                                                                                        tAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTaa                                       <  2:3274741/148‑1 (MQ=255)
                                                                                                        tAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTaa                                       <  1:1316746/148‑1 (MQ=255)
                                                                                                         aGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAg                                      <  1:207227/148‑1 (MQ=255)
                                                                                                          gTTGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGa                                     <  2:2658361/148‑1 (MQ=255)
                                                                                                            tGCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATc                                   >  2:1027629/1‑148 (MQ=255)
                                                                                                             gCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATcc                                  <  2:3043994/148‑1 (MQ=255)
                                                                                                             gCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATcc                                  <  1:2569469/148‑1 (MQ=255)
                                                                                                                 tATCATTGACTAAACTGTTTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGt                              <  2:2541141/148‑1 (MQ=255)
                                                                                                                 tATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGt                              <  2:1405350/148‑1 (MQ=255)
                                                                                                                 tATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGt                              <  1:2722949/148‑1 (MQ=255)
                                                                                                                 tATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGt                              >  2:1754969/1‑148 (MQ=255)
                                                                                                                       tGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAAt                        >  1:182678/1‑148 (MQ=255)
                                                                                                                           tAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtc                    <  2:515862/148‑1 (MQ=255)
                                                                                                                            aaaCTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                   >  2:1966674/1‑148 (MQ=255)
                                                                                                                            aaaCTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                   <  1:138032/148‑1 (MQ=255)
                                                                                                                             aaCTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtctc                  >  2:2160800/1‑148 (MQ=255)
                                                                                                                              aCTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCg                 >  2:346897/1‑148 (MQ=255)
                                                                                                                               cTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATTAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGt                >  1:1688314/1‑148 (MQ=255)
                                                                                                                               cTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATTAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGt                >  1:1688302/1‑148 (MQ=255)
                                                                                                                                    tATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCAt           >  2:3189984/1‑148 (MQ=255)
                                                                                                                                           aTACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTa    <  1:2632589/148‑1 (MQ=255)
                                                                                                                                            tACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTaa   >  1:3411800/1‑148 (MQ=255)
                                                                                                                                            tACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTaa   >  2:2069448/1‑148 (MQ=255)
                                                                                                                                            tACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTaa   <  1:2977225/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATTGAGTTGTTTTACTTGAATTAATACACCAGCGAGCTCTGATAGCTGCCAGAGTGGGTTCTATGGATACAGTTGCTGGGGTTGCCATTGGCTTATCTCCTGTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAAC  >  CP009273/291203‑291491

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: