Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 560,919 G→T 7.8% 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,9190GT 7.8% 1116.5 / 42.5 301V172V (GTC→GTAintDDLP12 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base G (140/137);  new base T (13/11);  total (153/148)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

AGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAAC  >  CP009273/560770‑561068
                                                                                                                                                     |                                                                                                                                                     
aGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTtgatg                                                                                                                                                       >  1:692900/1‑150 (MQ=255)
  aTTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGAc                                                                                                                                                     >  2:1602735/1‑150 (MQ=255)
  aTTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGAc                                                                                                                                                     <  2:1388683/150‑1 (MQ=255)
      atgatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgc                                                                                                                                                 >  1:2169721/1‑150 (MQ=255)
       tgatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgcg                                                                                                                                                <  2:3701726/150‑1 (MQ=255)
        gatgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGgcgc                                                                                                                                               >  1:1433543/1‑150 (MQ=255)
         atgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCt                                                                                                                                              <  1:1519310/150‑1 (MQ=255)
          tgTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCtt                                                                                                                                             >  1:1381359/1‑150 (MQ=255)
           gTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCttt                                                                                                                                            >  2:4024537/1‑150 (MQ=255)
           gTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCttt                                                                                                                                            <  1:4179330/150‑1 (MQ=255)
              cGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTc                                                                                                                                         >  1:3892177/1‑150 (MQ=255)
               gACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTcc                                                                                                                                        >  2:1072585/1‑150 (MQ=255)
               gACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTcc                                                                                                                                        >  1:2467980/1‑150 (MQ=255)
                 cTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAg                                                                                                                                      >  1:1082199/1‑150 (MQ=255)
                 cTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCGTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAg                                                                                                                                      <  2:544998/150‑1 (MQ=255)
                 cTTTCTCAGACCAGGTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAg                                                                                                                                      <  2:3823939/150‑1 (MQ=255)
                    tctcAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCa                                                                                                                                   >  2:1939762/1‑150 (MQ=255)
                    tctcAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCCCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCa                                                                                                                                   <  1:2940317/150‑1 (MQ=255)
                       cAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACtt                                                                                                                                <  2:38943/150‑1 (MQ=255)
                       cAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACtt                                                                                                                                >  1:335919/1‑150 (MQ=255)
                         gACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTcc                                                                                                                              <  1:1863280/150‑1 (MQ=255)
                            cAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGt                                                                                                                           >  2:2306871/1‑150 (MQ=255)
                              gTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTcg                                                                                                                         >  2:1833026/1‑150 (MQ=255)
                               ttGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTcgc                                                                                                                        >  1:51526/1‑150 (MQ=255)
                                tGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTcgcg                                                                                                                       >  1:4428475/1‑150 (MQ=255)
                                  ccAGCGCAAACTTGACGACATATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTATGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGtt                                                                                                                     >  1:867413/1‑150 (MQ=255)
                                   cAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTc                                                                                                                    >  2:4128492/1‑150 (MQ=255)
                                   cAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTc                                                                                                                    >  1:4147417/1‑150 (MQ=255)
                                         aaaCTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTgcagc                                                                                                              >  1:960341/1‑150 (MQ=255)
                                           aCTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCAc                                                                                                            <  1:398231/150‑1 (MQ=255)
                                            cTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACg                                                                                                           >  1:4322681/1‑150 (MQ=255)
                                               gacgacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgag                                                                                                        >  2:3101536/1‑150 (MQ=255)
                                               gacgacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgag                                                                                                        <  1:1063815/150‑1 (MQ=255)
                                               gacgacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgag                                                                                                        >  2:3101439/1‑150 (MQ=255)
                                                  gacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGgagact                                                                                                     >  2:4363271/1‑148 (MQ=255)
                                                  gacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAAt                                                                                                     <  2:2835337/150‑1 (MQ=255)
                                                     aGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGc                                                                                                  <  1:193121/150‑1 (MQ=255)
                                                        tttCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCtt                                                                                               <  1:4501953/150‑1 (MQ=255)
                                                         ttCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCttt                                                                                              >  2:2700377/1‑150 (MQ=255)
                                                            aGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATc                                                                                           <  2:641304/150‑1 (MQ=255)
                                                            aGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATc                                                                                           >  1:256516/1‑150 (MQ=255)
                                                            aGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATc                                                                                           >  2:1712889/1‑150 (MQ=255)
                                                              tGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCAt                                                                                         <  1:2885466/150‑1 (MQ=255)
                                                                    ggggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAg                                                                                   >  2:767355/1‑150 (MQ=255)
                                                                     gggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGa                                                                                  >  1:161132/1‑150 (MQ=255)
                                                                      ggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGAt                                                                                 <  1:1293613/150‑1 (MQ=255)
                                                                      ggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGAt                                                                                 >  1:2831070/1‑150 (MQ=255)
                                                                      ggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGAt                                                                                 <  1:1293760/150‑1 (MQ=255)
                                                                         aCTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCt                                                                              <  2:2839273/150‑1 (MQ=255)
                                                                          cTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCtt                                                                             <  2:3187239/150‑1 (MQ=255)
                                                                           tcatcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCttt                                                                            <  2:692900/150‑1 (MQ=255)
                                                                           tcatcaATAAGGCGTTGTGCGTCCTCCTGTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCttt                                                                            <  2:4194064/150‑1 (MQ=255)
                                                                            catcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTg                                                                           >  2:1222685/1‑150 (MQ=255)
                                                                             atcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGc                                                                          >  1:3895321/1‑150 (MQ=255)
                                                                               caATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCt                                                                        <  2:886477/150‑1 (MQ=255)
                                                                                 aTAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCttt                                                                      <  1:3550948/150‑1 (MQ=255)
                                                                                  tAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTg                                                                     >  1:4161459/1‑150 (MQ=255)
                                                                                  tAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTg                                                                     >  1:723662/1‑150 (MQ=255)
                                                                                   aaGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGt                                                                    >  1:64017/1‑150 (MQ=255)
                                                                                    aGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTc                                                                   >  2:3941966/1‑150 (MQ=255)
                                                                                       cGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGa                                                                <  1:2416683/150‑1 (MQ=255)
                                                                                        gTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAg                                                               <  2:2256123/150‑1 (MQ=255)
                                                                                         ttttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGt                                                              <  1:2926144/150‑1 (MQ=255)
                                                                                          tttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGtt                                                             <  1:3445422/150‑1 (MQ=255)
                                                                                           ttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTg                                                            <  1:3570986/150‑1 (MQ=255)
                                                                                           ttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTg                                                            >  1:2396975/1‑150 (MQ=255)
                                                                                                  ctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTgg                                                     <  2:2615651/150‑1 (MQ=255)
                                                                                                    cctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCt                                                   >  2:2508727/1‑150 (MQ=255)
                                                                                                    cctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCt                                                   <  1:2074182/150‑1 (MQ=255)
                                                                                                     ctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCtt                                                  >  2:4423313/1‑150 (MQ=255)
                                                                                                     ctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCGTTGTCTGAGTTGATACTGGCtt                                                  >  2:2183815/1‑150 (MQ=255)
                                                                                                     ctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCtt                                                  <  2:3739591/150‑1 (MQ=255)
                                                                                                           cAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTc                                            >  1:3246534/1‑150 (MQ=255)
                                                                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCATTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                           >  2:2377042/1‑150 (MQ=255)
                                                                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTATCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                           >  2:715160/1‑150 (MQ=255)
                                                                                                              ccATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATa                                         <  1:2253697/150‑1 (MQ=255)
                                                                                                                 tCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAAc                                      <  1:1478384/150‑1 (MQ=255)
                                                                                                                  cTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACa                                     >  2:179391/1‑150 (MQ=255)
                                                                                                                   tGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAg                                    <  2:1554276/150‑1 (MQ=255)
                                                                                                                     aCTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGc                                  >  2:1544383/1‑150 (MQ=255)
                                                                                                                     aCTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGc                                  >  2:3377806/1‑150 (MQ=255)
                                                                                                                       tCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTATGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAg                                <  1:4353962/150‑1 (MQ=255)
                                                                                                                        cGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGt                               >  1:2845967/1‑150 (MQ=255)
                                                                                                                         gCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGtt                              <  2:2010663/150‑1 (MQ=255)
                                                                                                                         gCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGtt                              >  1:4412346/1‑150 (MQ=255)
                                                                                                                             gTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCttt                          <  2:3006959/150‑1 (MQ=255)
                                                                                                                               ttCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTAc                        >  1:1563720/1‑150 (MQ=255)
                                                                                                                               ttCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTAc                        >  2:1138497/1‑150 (MQ=255)
                                                                                                                                  tGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCtt                     <  2:3831405/150‑1 (MQ=255)
                                                                                                                                  tGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCtt                     <  2:2784162/150‑1 (MQ=255)
                                                                                                                                   gACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCATCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCttt                    <  2:3469912/150‑1 (MQ=255)
                                                                                                                                      cGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCt                 <  1:4342118/150‑1 (MQ=255)
                                                                                                                                      cGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCt                 >  2:2527136/1‑150 (MQ=255)
                                                                                                                                         tGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt              >  1:2634477/1‑150 (MQ=255)
                                                                                                                                          ggTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCTGTTCTTTTCCTTTCCTGATg             <  2:3667210/150‑1 (MQ=255)
                                                                                                                                            tATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggc           >  2:2889940/1‑150 (MQ=255)
                                                                                                                                             aTCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCTTGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg          <  2:3805246/150‑1 (MQ=255)
                                                                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCCGTTCTTTACCTTTCCTGATggcg          >  2:505410/1‑150 (MQ=255)
                                                                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg          <  1:4266645/150‑1 (MQ=255)
                                                                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg          <  1:2182317/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGATTTATGATGTTCGACTTTCTCAGACCAGTTGCCAGCGCAAACTTGACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAAC  >  CP009273/560770‑561068

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: