Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 560,966 G→T 10.0% L157I (CTC→ATC)  intD ← DLP12 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273560,9660GT10.0% 990.2 / 52.6 267L157I (CTC→ATC) intDDLP12 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base G (100/139);  new base T (15/13);  total (115/152)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

ACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTGCTT  >  CP009273/560818‑561114
                                                                                                                                                    |                                                                                                                                                    
  gacAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAAt                                                                                                                                                   >  1:1256256/1‑150 (MQ=255)
   acAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATg                                                                                                                                                  >  1:582735/1‑150 (MQ=255)
    cAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATgg                                                                                                                                                 >  1:1161155/1‑150 (MQ=255)
      gATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGcc                                                                                                                                               <  2:1627720/150‑1 (MQ=255)
          tCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTa                                                                                                                                           >  1:2228983/1‑150 (MQ=255)
           cAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTAt                                                                                                                                          <  1:3960060/150‑1 (MQ=255)
              tGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCAt                                                                                                                                       >  1:2153502/1‑150 (MQ=255)
              tGGGTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCAt                                                                                                                                       <  1:1564450/150‑1 (MQ=255)
               ggTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCATTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATg                                                                                                                                      <  1:632341/150‑1 (MQ=255)
                    ggggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAg                                                                                                                                 >  1:3197123/1‑150 (MQ=255)
                    ggggCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAg                                                                                                                                 <  2:1678433/150‑1 (MQ=255)
                         aCTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCt                                                                                                                            <  2:1660149/150‑1 (MQ=255)
                          cTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCtt                                                                                                                           <  1:3748811/150‑1 (MQ=255)
                           tcatcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCttt                                                                                                                          <  1:1001697/150‑1 (MQ=255)
                            catcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTg                                                                                                                         >  2:2059729/1‑150 (MQ=255)
                            catcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTg                                                                                                                         <  1:3419941/150‑1 (MQ=255)
                            catcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTg                                                                                                                         >  2:862208/1‑150 (MQ=255)
                            catcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTg                                                                                                                         <  1:4046328/150‑1 (MQ=255)
                             atcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGc                                                                                                                        >  2:3533973/1‑150 (MQ=255)
                             atcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGc                                                                                                                        >  2:748851/1‑150 (MQ=255)
                             atcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGcc                                                                                                                       <  2:2996446/150‑1 (MQ=255)
                              tcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGcc                                                                                                                       >  1:1775246/1‑150 (MQ=255)
                              tcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGcc                                                                                                                       <  2:1224012/150‑1 (MQ=255)
                              tcaATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGcc                                                                                                                       >  2:602928/1‑150 (MQ=255)
                               caATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCt                                                                                                                      <  1:1718951/150‑1 (MQ=255)
                               caATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCt                                                                                                                      >  1:2898937/1‑150 (MQ=255)
                                aaTAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCtt                                                                                                                     >  1:102996/1‑150 (MQ=255)
                                     ggCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCt                                                                                                                <  2:1041137/150‑1 (MQ=255)
                                      gCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTg                                                                                                               <  2:900578/150‑1 (MQ=255)
                                       cGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGa                                                                                                              <  1:2657412/150‑1 (MQ=255)
                                        gTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAg                                                                                                             >  2:1459718/1‑150 (MQ=255)
                                          tttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGtt                                                                                                           <  1:2137976/150‑1 (MQ=255)
                                           ttGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTg                                                                                                          >  1:3354887/1‑150 (MQ=255)
                                            tGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGa                                                                                                         <  2:1845239/150‑1 (MQ=255)
                                            tGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGa                                                                                                         <  1:1677691/150‑1 (MQ=255)
                                              cTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATa                                                                                                       <  1:3824134/150‑1 (MQ=255)
                                                tcctccATTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACt                                                                                                     <  1:2716490/150‑1 (MQ=255)
                                                 cctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTg                                                                                                    <  2:2331715/150‑1 (MQ=255)
                                                  ctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTgg                                                                                                   <  2:2098294/150‑1 (MQ=255)
                                                  ctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTgg                                                                                                   <  2:3324288/150‑1 (MQ=255)
                                                  ctcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTGTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTgg                                                                                                   >  2:2854832/1‑150 (MQ=255)
                                                   tcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGc                                                                                                  >  1:3541083/1‑150 (MQ=255)
                                                   tcctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGc                                                                                                  >  1:681024/1‑150 (MQ=255)
                                                    cctTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCt                                                                                                 <  2:732617/150‑1 (MQ=255)
                                                          ccAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGtt                                                                                           <  1:126474/150‑1 (MQ=255)
                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                                                                         >  2:4010646/1‑150 (MQ=255)
                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                                                                         >  1:2036043/1‑150 (MQ=255)
                                                            aGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCa                                                                                         >  2:1344124/1‑150 (MQ=255)
                                                              ccATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATa                                                                                       <  2:2741059/150‑1 (MQ=255)
                                                               cATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATaa                                                                                      <  1:2994521/150‑1 (MQ=255)
                                                                aTCTGTCTCGCTTGTTTCTGACCGCTGGTATCTTGATGAGAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCAGGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTTGCTTTGGTTCATaca                                                                                     <  2:1471288/150‑3 (MQ=255)
                                                                aTCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATaaa                                                                                     <  2:489548/150‑1 (MQ=255)
                                                                   tGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAg                                                                                  <  2:1762337/150‑1 (MQ=255)
                                                                     aCTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGc                                                                                >  2:3014894/1‑150 (MQ=255)
                                                                      cTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCa                                                                               >  2:2165194/1‑150 (MQ=255)
                                                                      cTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCa                                                                               <  1:1284686/150‑1 (MQ=255)
                                                                        cGCTTGTTTCTGACCGCTGGTATCTTGATTACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGt                                                                             <  2:3037544/150‑1 (MQ=255)
                                                                        cGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGt                                                                             >  2:1131236/1‑150 (MQ=255)
                                                                         gCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGtg                                                                            <  1:3504653/150‑2 (MQ=255)
                                                                           ttGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCt                                                                          <  1:3700810/150‑1 (MQ=255)
                                                                             gTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCttt                                                                        <  1:3465709/150‑1 (MQ=255)
                                                                                     ccGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTcc                                                                <  2:699701/150‑1 (MQ=255)
                                                                                         tGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGAt                                                            <  2:830770/150‑1 (MQ=255)
                                                                                            tATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggc                                                         >  2:3696618/1‑150 (MQ=255)
                                                                                             aTCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATggcg                                                        <  2:3222600/150‑1 (MQ=255)
                                                                                                    tgaCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGaa                                                 <  1:3407033/150‑1 (MQ=255)
                                                                                                    tgaCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGaa                                                 >  1:499595/1‑150 (MQ=255)
                                                                                                     gaCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAAc                                                <  1:315666/150‑1 (MQ=255)
                                                                                                      aCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACt                                               >  2:1547926/1‑150 (MQ=255)
                                                                                                      aCAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACt                                               <  2:2501533/150‑1 (MQ=255)
                                                                                                        aGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACttt                                             >  1:1737114/1‑150 (MQ=255)
                                                                                                                tttCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTcc                                     <  1:865187/150‑1 (MQ=255)
                                                                                                                 ttCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCa                                    <  1:2819931/150‑1 (MQ=255)
                                                                                                                    cAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCatat                                 >  2:2423311/1‑150 (MQ=255)
                                                                                                                    cAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCatat                                 <  1:2493690/150‑1 (MQ=255)
                                                                                                                    cAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCatat                                 >  2:3357108/1‑150 (MQ=255)
                                                                                                                     aGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATAtt                                <  2:2869083/150‑1 (MQ=255)
                                                                                                                      gCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATAttt                               <  1:415739/150‑1 (MQ=255)
                                                                                                                        cACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATAtttct                             <  2:2814914/150‑1 (MQ=255)
                                                                                                                           ttCCAGTCGCGTTCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATAtttctttc                          <  1:1506449/150‑1 (MQ=255)
                                                                                                                           ttCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATAtttctttc                          <  1:747989/150‑1 (MQ=255)
                                                                                                                            tCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCg                         >  1:1605277/1‑150 (MQ=255)
                                                                                                                            tCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCg                         >  1:1117939/1‑150 (MQ=255)
                                                                                                                             ccAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGt                        <  1:1074222/150‑1 (MQ=255)
                                                                                                                              cAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTc                       >  2:3655275/1‑150 (MQ=255)
                                                                                                                              cAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTc                       <  1:3326112/150‑1 (MQ=255)
                                                                                                                                gTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCtt                     <  2:757199/150‑1 (MQ=255)
                                                                                                                                 tCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCttt                    <  1:3020706/150‑1 (MQ=255)
                                                                                                                                 tCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCttt                    >  1:2366227/1‑150 (MQ=255)
                                                                                                                                  cgcgTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTc                   >  1:1232039/1‑150 (MQ=255)
                                                                                                                                    cgTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTg                 <  1:191067/150‑1 (MQ=255)
                                                                                                                                       tCTGCAGCACGGATAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTa              <  1:769390/150‑1 (MQ=255)
                                                                                                                                        cTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTAt             <  2:3081818/150‑1 (MQ=255)
                                                                                                                                            agcaCGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCAt         <  2:3887018/150‑1 (MQ=255)
                                                                                                                                               aCGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCt      >  2:1117816/1‑150 (MQ=255)
                                                                                                                                               aCGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCt      >  2:2126635/1‑150 (MQ=255)
                                                                                                                                                cGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTg     >  2:2475281/1‑150 (MQ=255)
                                                                                                                                                cGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTg     <  1:3658986/150‑1 (MQ=255)
                                                                                                                                                cGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTg     >  1:1860801/1‑150 (MQ=255)
                                                                                                                                                cGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTg     >  2:3483120/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGACAGATTTCAGTGGTTCGGGGCACTCATCAATAAGGCGTTTTGCTTCCTCCTTTTCCAGCCATCTGACTCGCTTGTTTCTGACCGCTGGTATCTTGATGACAGGCGCTTTTTCCAGCCACTTCCAGTCGCGTTCTGCAGCACGGAGAATGGCCTTTATCATGGCAAGATGCTTTGCCTTTGTCTGAGTTGATACTGGCTTTGGTTCATAAACAGGCAGTTCTTTACCTTTCCTGATGGCGGCCTGAACTTTCTGTTTCCATATTTCTTTCGTCTTTCTGTTATGCATTCTGCTT  >  CP009273/560818‑561114

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: