Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 560,919 G→T 10.9% V172V (GTC→GTA intD ← DLP12 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273560,9190GT10.9% 977.1 / 68.0 282V172V (GTC→GTAintDDLP12 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base G (128/123);  new base T (18/13);  total (146/136)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 299 total aligned reads displayed.

AGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAA  >  CP009273/560770‑561067
                                                                                                                                                     |                                                                                                                                                    
aGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTtgatg                                                                                                                                                      <  1:4015494/150‑1 (MQ=255)
aGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTtgatg                                                                                                                                                      <  2:4022269/150‑1 (MQ=255)
aGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTtgatg                                                                                                                                                      >  1:998538/1‑150 (MQ=255)
 gATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTtgatga                                                                                                                                                     >  1:2088876/1‑150 (MQ=255)
   tttATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACa                                                                                                                                                   <  1:2943568/150‑1 (MQ=255)
   tttATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACa                                                                                                                                                   <  1:1326045/150‑1 (MQ=255)
   tttATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACa                                                                                                                                                   <  2:311277/150‑1 (MQ=255)
     tATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAgg                                                                                                                                                 <  1:3366909/150‑1 (MQ=255)
       tgatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgcg                                                                                                                                               >  2:28637/1‑150 (MQ=255)
        gatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgcgc                                                                                                                                              >  2:2597952/1‑150 (MQ=255)
        gatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgcgc                                                                                                                                              >  1:2312644/1‑150 (MQ=255)
         atgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCt                                                                                                                                             >  1:3761014/1‑150 (MQ=255)
         atgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCt                                                                                                                                             <  2:3919373/150‑1 (MQ=255)
             tCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCttttt                                                                                                                                         <  1:2131381/150‑1 (MQ=255)
             tCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCttttt                                                                                                                                         >  2:348595/1‑150 (MQ=255)
             tCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCttttt                                                                                                                                         <  1:911237/150‑1 (MQ=255)
              cGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTc                                                                                                                                        >  2:1914675/1‑150 (MQ=255)
                 cTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAg                                                                                                                                     >  2:3608913/1‑150 (MQ=255)
                   ttCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGcc                                                                                                                                   >  2:4660129/1‑150 (MQ=255)
                    tctcAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCa                                                                                                                                  <  1:2015171/150‑1 (MQ=255)
                      tcAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACt                                                                                                                                <  2:4597763/150‑1 (MQ=255)
                       cAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACtt                                                                                                                               >  1:1400776/1‑150 (MQ=255)
                        aGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTc                                                                                                                              <  2:818746/150‑1 (MQ=255)
                         gACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTcc                                                                                                                             >  2:1817870/1‑150 (MQ=255)
                          aCCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCa                                                                                                                            <  1:4319216/150‑1 (MQ=255)
                             aGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTc                                                                                                                         <  2:4304810/150‑1 (MQ=255)
                                  ccAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGtt                                                                                                                    <  1:1139347/150‑1 (MQ=255)
                                   cAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTc                                                                                                                   >  1:1952432/1‑150 (MQ=255)
                                          aaCTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTgcagca                                                                                                            >  2:355201/1‑150 (MQ=255)
                                           aCTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCAc                                                                                                           <  1:2974434/150‑1 (MQ=255)
                                              tGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGga                                                                                                        <  1:518378/150‑1 (MQ=255)
                                               gacgacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgag                                                                                                       >  1:1079035/1‑150 (MQ=255)
                                               gacgacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgag                                                                                                       >  1:2994782/1‑150 (MQ=255)
                                                  gacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAAt                                                                                                    <  1:2703001/150‑1 (MQ=255)
                                                   acAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATg                                                                                                   <  2:737316/150‑1 (MQ=255)
                                                    cAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATgg                                                                                                  <  2:638962/150‑1 (MQ=255)
                                                    cAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATgg                                                                                                  <  1:2800202/150‑1 (MQ=255)
                                                       aTTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCt                                                                                               <  2:2314474/150‑1 (MQ=255)
                                                       aTTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGCTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCt                                                                                               <  2:726536/150‑1 (MQ=255)
                                                            aGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATc                                                                                          <  2:155122/150‑1 (MQ=255)
                                                              tGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCAt                                                                                        <  2:391520/150‑1 (MQ=255)
                                                              tGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCAt                                                                                        <  2:270807/150‑1 (MQ=255)
                                                               ggTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATg                                                                                       <  1:347066/150‑1 (MQ=255)
                                                               ggTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATg                                                                                       >  2:2039881/1‑150 (MQ=255)
                                                                gTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATgg                                                                                      <  2:3738258/150‑1 (MQ=255)
                                                                   cGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCaa                                                                                   <  2:1802587/150‑1 (MQ=255)
                                                                    ggggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAg                                                                                  >  2:2068677/1‑150 (MQ=255)
                                                                       gCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATg                                                                               <  2:4627568/150‑1 (MQ=255)
                                                                             atcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTACAGCCACTTCCAGTCGCGTTCTTCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGc                                                                         >  1:2064327/1‑150 (MQ=255)
                                                                                aaTAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCtt                                                                      >  1:142409/1‑150 (MQ=255)
                                                                                aaTAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCtt                                                                      <  1:193187/150‑1 (MQ=255)
                                                                                  tAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTACAGCCACATACAGTCCCGTTCTGACGAACGGAGAATGGCCGTTATCATTGCAAGATGCGTTTCCATTg                                                                    >  1:1233342/1‑150 (MQ=255)
                                                                                     ggCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCt                                                                 <  1:3703139/150‑1 (MQ=255)
                                                                                        gTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATTGCAAGATGCTTTGCCTTTGTCTGAg                                                              >  1:575760/1‑150 (MQ=255)
                                                                                        gTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAg                                                              <  1:3671513/150‑1 (MQ=255)
                                                                                           ttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTg                                                           >  2:1095589/1‑150 (MQ=255)
                                                                                           ttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTg                                                           >  1:1734236/1‑150 (MQ=255)
                                                                                              cTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATAATGGCAAGATGCTTTTCCTTTGTCTGAGTTGATa                                                        >  1:1663153/1‑150 (MQ=255)
                                                                                                  ctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTgg                                                    >  1:218146/1‑150 (MQ=255)
                                                                                                    cctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCt                                                  <  1:1838837/150‑1 (MQ=255)
                                                                                                        ttCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATAATGGCAAGATGCTTTTCCTTTGTCTGAGTTGATAATGGCTTTgt                                              >  2:3268127/1‑149 (MQ=255)
                                                                                                        ttCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTgg                                              >  2:3938476/1‑150 (MQ=255)
                                                                                                          ccAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGtt                                            <  1:129702/150‑1 (MQ=255)
                                                                                                          ccAGCCATCTGACTCGCGTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGtt                                            <  2:4466836/150‑1 (MQ=255)
                                                                                                           cAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTc                                           <  1:687384/150‑1 (MQ=255)
                                                                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                          >  1:1985279/1‑150 (MQ=255)
                                                                                                               cATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATaa                                       >  2:2446569/1‑150 (MQ=255)
                                                                                                               cATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATAATGGCAAGATTCTTTTCCTTTGTCTGAGTTGATACTTGCTTTGGTTCATaa                                       >  1:2493566/1‑150 (MQ=255)
                                                                                                                aTCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATAATGGCTTTGGTTCATaaa                                      >  1:1023389/1‑150 (MQ=255)
                                                                                                                    gACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAgg                                  >  1:3833566/1‑150 (MQ=255)
                                                                                                                     aCTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGc                                 >  1:2170702/1‑150 (MQ=255)
                                                                                                                       tCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAg                               <  2:854139/150‑1 (MQ=255)
                                                                                                                        cGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGt                              <  2:3442286/150‑1 (MQ=255)
                                                                                                                         gCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGtt                             >  2:516428/1‑150 (MQ=255)
                                                                                                                          cTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTc                            <  2:3033913/150‑1 (MQ=255)
                                                                                                                          cTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTc                            >  1:630459/1‑150 (MQ=255)
                                                                                                                            tGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCtt                          >  1:144747/1‑150 (MQ=255)
                                                                                                                            tGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGACTGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCtt                          >  2:1739013/1‑150 (MQ=255)
                                                                                                                            tGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCtt                          >  2:3595111/1‑150 (MQ=255)
                                                                                                                              tttCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTa                        <  2:4245815/150‑1 (MQ=255)
                                                                                                                               gtCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTAc                       <  2:4314233/149‑1 (MQ=255)
                                                                                                                                tCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTAcc                      <  1:3938476/150‑1 (MQ=255)
                                                                                                                                 cTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCGTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCt                     >  1:611088/1‑150 (MQ=255)
                                                                                                                                  tGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCtt                    <  2:2171246/150‑1 (MQ=255)
                                                                                                                                  tGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCtt                    >  2:453021/1‑150 (MQ=255)
                                                                                                                                    aCCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTc                  >  1:4689173/1‑150 (MQ=255)
                                                                                                                                      cGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCt                <  1:2274419/150‑1 (MQ=255)
                                                                                                                                       gCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTg               <  1:854191/150‑1 (MQ=255)
                                                                                                                                        cTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGa              >  1:4519433/1‑150 (MQ=255)
                                                                                                                                         tGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt             <  2:949887/150‑1 (MQ=255)
                                                                                                                                         tGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt             <  2:3344494/150‑1 (MQ=255)
                                                                                                                                         tGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt             <  2:1311510/150‑1 (MQ=255)
                                                                                                                                         tGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt             >  2:1261559/1‑150 (MQ=255)
                                                                                                                                         tGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt             <  1:435905/150‑1 (MQ=255)
                                                                                                                                          ggTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATg            <  2:4642187/150‑1 (MQ=255)
                                                                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg         >  1:1580304/1‑150 (MQ=255)
                                                                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg         <  1:255589/150‑1 (MQ=255)
                                                                                                                                                ttGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGcc      >  2:4192225/1‑150 (MQ=255)
                                                                                                                                                 tgatgaCAGGCGCTCTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCt     >  2:2713808/1‑150 (MQ=255)
                                                                                                                                                  gattaCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGGTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTTCCTTTGTCTGAGTTGATACTGGCTGTGGTTCATAAAAAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTg    >  2:4018475/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAA  >  CP009273/560770‑561067

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: