Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,470,873 G→T 84.1% R928S (CGT→AGT)  valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,470,8730GT84.1% 643.0 / 147.1 290R928S (CGT→AGT) valSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (32/14);  new base T (127/117);  total (159/131)
Fisher's exact test for biased strand distribution p-value = 3.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 290 total aligned reads displayed.

TTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCA  >  CP009273/4470729‑4471019
                                                                                                                                                |                                                                                                                                                  
ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCtt                                                                                                                                               >  2:3455083/1‑150 (MQ=255)
  aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTc                                                                                                                                             <  1:4007263/150‑1 (MQ=255)
   ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCGTCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTcg                                                                                                                                            >  2:698912/1‑150 (MQ=255)
   ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTcg                                                                                                                                            >  2:112840/1‑150 (MQ=255)
     aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTcgcg                                                                                                                                          <  1:4386534/150‑1 (MQ=255)
       ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGAt                                                                                                                                        >  1:1316038/1‑150 (MQ=255)
          tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGAc                                                                                                                                     <  2:208509/150‑1 (MQ=255)
          tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCGTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGAc                                                                                                                                     >  1:2644361/1‑150 (MQ=255)
           gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGAcc                                                                                                                                    <  2:4304172/150‑1 (MQ=255)
           gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGAcc                                                                                                                                    >  2:3561747/1‑150 (MQ=255)
           gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGAcc                                                                                                                                    >  2:270354/1‑150 (MQ=255)
             aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGc                                                                                                                                  >  1:180814/1‑150 (MQ=255)
               gataTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCtt                                                                                                                                <  1:813370/147‑1 (MQ=255)
                agaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACGGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTc                                                                                                                               >  2:812612/1‑150 (MQ=255)
                 gaTAAAAATGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACAGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCACTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTAGCGATGACCGCTTcc                                                                                                                              >  2:664413/1‑150 (MQ=255)
                 gaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTcc                                                                                                                              >  2:3823328/1‑150 (MQ=255)
                  aTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCCCTTTCGCGATGACCGCTTCCg                                                                                                                             <  2:1249043/150‑1 (MQ=255)
                  aTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGGTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCg                                                                                                                             >  2:3100998/1‑150 (MQ=255)
                   tAAAAAAGGCCGGAGCATGCTCCGGCCGTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCgg                                                                                                                            >  2:3133988/1‑150 (MQ=255)
                    aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGt                                                                                                                           <  1:2874039/150‑1 (MQ=255)
                      aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTgc                                                                                                                         >  1:599711/1‑150 (MQ=255)
                      aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTgc                                                                                                                         <  2:262980/150‑1 (MQ=255)
                      aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTgc                                                                                                                         >  2:2998219/1‑150 (MQ=255)
                      aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTgc                                                                                                                         <  1:3266387/150‑1 (MQ=255)
                         aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGcgcg                                                                                                                      <  1:2123783/150‑1 (MQ=255)
                            ccGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGc                                                                                                                   >  2:1837444/1‑150 (MQ=255)
                             cGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCg                                                                                                                  <  1:2466671/150‑1 (MQ=255)
                              ggAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGa                                                                                                                 >  2:4072354/1‑150 (MQ=255)
                              ggAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGa                                                                                                                 >  1:3364626/1‑150 (MQ=255)
                                 gCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACaa                                                                                                              >  1:3423086/1‑150 (MQ=255)
                                   aTGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAg                                                                                                            >  2:2415469/1‑150 (MQ=255)
                                       tCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCtt                                                                                                        <  1:759204/150‑1 (MQ=255)
                                           gCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCttcgtt                                                                                                    <  2:1259551/150‑1 (MQ=255)
                                           gCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCttcgtt                                                                                                    >  2:29248/1‑150 (MQ=255)
                                            ccTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCttcgttc                                                                                                   <  2:443838/150‑1 (MQ=255)
                                                cGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCa                                                                                               <  1:2263330/150‑1 (MQ=255)
                                                cGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCa                                                                                               >  1:3646908/1‑150 (MQ=255)
                                                 gTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAg                                                                                              >  2:3783987/1‑150 (MQ=255)
                                                       atcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGtt                                                                                        >  1:4082246/1‑150 (MQ=255)
                                                        tcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtc                                                                                       >  2:2314940/1‑150 (MQ=255)
                                                        tcaCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtc                                                                                       <  2:3929381/150‑1 (MQ=255)
                                                         caCTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTtct                                                                                      <  1:125475/150‑1 (MQ=255)
                                                           cTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCg                                                                                    >  2:1979690/1‑150 (MQ=255)
                                                            tgtgTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTAGa                                                                                   >  2:3055365/1‑150 (MQ=255)
                                                             gtgtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGAt                                                                                  <  2:4494731/150‑1 (MQ=255)
                                                             gtgtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGAt                                                                                  <  1:436732/150‑1 (MQ=255)
                                                               gtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATAc                                                                                >  2:581667/1‑150 (MQ=255)
                                                                   tGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCt                                                                            <  1:1009354/150‑1 (MQ=255)
                                                                      ttACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGAt                                                                         >  2:3302811/1‑150 (MQ=255)
                                                                        aCAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGAttt                                                                       <  2:3655983/150‑1 (MQ=255)
                                                                               ggcgATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTc                                                                >  1:1924116/1‑150 (MQ=255)
                                                                                gcgATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCa                                                               <  2:1683501/150‑1 (MQ=255)
                                                                                 cgATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCaa                                                              >  2:4416698/1‑150 (MQ=255)
                                                                                  gATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAAt                                                             <  1:418948/150‑1 (MQ=255)
                                                                                    tAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCt                                                           >  1:2937860/1‑150 (MQ=255)
                                                                                     aaCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCtt                                                          >  1:1683948/1‑150 (MQ=255)
                                                                                     aaCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCtt                                                          >  1:734997/1‑150 (MQ=255)
                                                                                     aaCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCtt                                                          >  1:1597540/1‑150 (MQ=255)
                                                                                     aaCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCtt                                                          <  2:2520318/150‑1 (MQ=255)
                                                                                     aaCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCtt                                                          <  2:2806517/150‑1 (MQ=255)
                                                                                      aCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTc                                                         <  2:396002/150‑1 (MQ=255)
                                                                                       cAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCg                                                        >  2:4556473/1‑150 (MQ=255)
                                                                                          ccTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCa                                                     <  2:660738/150‑1 (MQ=255)
                                                                                          ccTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCa                                                     >  1:4063087/1‑150 (MQ=255)
                                                                                            tgctgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCAct                                                   <  1:2602757/150‑1 (MQ=255)
                                                                                               tgTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActtct                                                >  2:3520011/1‑150 (MQ=255)
                                                                                                gTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCActtctt                                               >  1:3881234/1‑150 (MQ=255)
                                                                                                 ttCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCttt                                              >  1:186206/1‑150 (MQ=255)
                                                                                                  tCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTc                                             <  1:3133988/150‑1 (MQ=255)
                                                                                                  tCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAACTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTc                                             <  1:2340638/150‑1 (MQ=255)
                                                                                                    aaTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGc                                           >  1:221774/1‑150 (MQ=255)
                                                                                                      tCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCa                                         >  1:3621438/1‑150 (MQ=255)
                                                                                                         gTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAc                                      <  2:2356651/150‑1 (MQ=255)
                                                                                                         gTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAc                                      >  1:2800044/1‑150 (MQ=255)
                                                                                                          tttcgctttcgcTTCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAcg                                     >  2:4586356/1‑150 (MQ=255)
                                                                                                            tcgctttcgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGAcgcg                                   <  2:750387/150‑1 (MQ=255)
                                                                                                             cgctttcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACgcgc                                  >  2:3681081/1‑150 (MQ=255)
                                                                                                               ctttcgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCa                                <  2:2191957/150‑1 (MQ=255)
                                                                                                               ctttcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCa                                >  2:742158/1‑150 (MQ=255)
                                                                                                                tttcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAg                               >  2:50632/1‑150 (MQ=255)
                                                                                                                   cgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTc                            >  1:150801/1‑150 (MQ=255)
                                                                                                                   cgcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGAGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTc                            <  1:1401248/150‑1 (MQ=255)
                                                                                                                   cgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTc                            <  2:3401278/150‑1 (MQ=255)
                                                                                                                    gcttCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCa                           >  2:4115356/1‑150 (MQ=255)
                                                                                                                     cttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAt                          >  1:983925/1‑150 (MQ=255)
                                                                                                                       tCCGCATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCAtct                        >  2:1931195/1‑150 (MQ=255)
                                                                                                                          gaaTAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCTCCACTTCTTTCGCCAGACGCGCCAGCTCAtcttct                     >  1:566337/3‑150 (MQ=255)
                                                                                                                           cATAGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATcttctt                    <  1:3100998/150‑1 (MQ=255)
                                                                                                                            aTAGCCTTCCAGCTTCTCACTCTCTTTCGCTATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCttt                   <  1:4277173/150‑1 (MQ=255)
                                                                                                                              aGCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTttgt                 <  2:1319065/150‑1 (MQ=255)
                                                                                                                               gCCTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTttgtt                <  2:1832441/150‑1 (MQ=255)
                                                                                                                                ccTTCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTttgttg               <  1:3033412/150‑1 (MQ=255)
                                                                                                                                  ttCCAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTtgat             <  2:1914899/150‑1 (MQ=255)
                                                                                                                                    ccAGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTtgatga           <  1:4622610/150‑1 (MQ=255)
                                                                                                                                      aGCTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAgg         >  2:1887327/1‑150 (MQ=255)
                                                                                                                                        cTTCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgcc       <  1:698872/150‑1 (MQ=255)
                                                                                                                                         ttCTCACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgcca      <  2:2005763/150‑1 (MQ=255)
                                                                                                                                             cACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgccagcca  >  2:4150058/1‑150 (MQ=255)
                                                                                                                                             cACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgccagcca  >  2:4042730/1‑150 (MQ=255)
                                                                                                                                             cACTCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGgccagcca  >  2:1730618/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCA  >  CP009273/4470729‑4471019

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: