Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 578,210 T→C 91.4% V35A (GTC→GCC)  rzpD → DLP12 prophage; putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,2100TC91.4% 887.5 / 66.2 541V35A (GTC→GCC) rzpDDLP12 prophage; putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (20/27);  new base C (279/215);  total (299/242)
Fisher's exact test for biased strand distribution p-value = 9.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 541 total aligned reads displayed.

GTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGA  >  NC_000913/578067‑578357
                                                                                                                                               |                                                                                                                                                   
         ttGAGCGTGAAGTCTGTTTGTGCCGCCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTCATTATCTGCATCATCGTCAGCCTGTCATGGGCGGGCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCACCGCGACaaa                                                                                                                                        >  5:4801/1‑148 (MQ=37)
                               ggggCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                                                                 <  4:255292/149‑1 (MQ=255)
                                    atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                                                                 >  4:315251/6‑128 (MQ=11)
                                      gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc                                                                                                                           <  1:217212/129‑1 (MQ=11)
                                         tGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATTGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                                                                                                        <  6:109236/148‑1 (MQ=11)
                                          gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                                                        <  7:210068/131‑1 (MQ=11)
                                           aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                                                              <  6:6182/124‑1 (MQ=11)
                                              agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                                                                 >  5:91822/1‑134 (MQ=255)
                                                      cgcgATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg                                                                                           >  8:169899/1‑148 (MQ=255)
                                                                         gtgaTCTTCATCATCGTCTGCCTGTCATGGGCGGGTAATCATTACCGTGATAACGCCATTACCTAAAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                                                                       <  4:188766/146‑1 (MQ=11)
                                                                                                       gCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                                           >  2:314392/1‑147 (MQ=255)
                                                                                                           tCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                                     <  7:352744/149‑1 (MQ=255)
                                                                                                           tCAATCATTACCGTGATAACGCAATCGACTACAAAGTCCAGCGCGAAAAAAATGCCAGAGAACTGAAGATAGCGAACGCGGCAATTACTGACATGAAGATGCGTCAGCTTGATGTTGCTGCTCTCGATGCAAAATACACGaa                                            >  8:217084/1‑142 (MQ=14)
                                                                                                                  ttACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGAAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATg                              >  4:359584/1‑149 (MQ=255)
                                                                                                                              cGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                                     <  1:303938/130‑1 (MQ=255)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  3:41785/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  3:264660/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  6:340813/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  5:359862/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  7:100067/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  4:124533/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  2:55661/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 <  1:241022/142‑1 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  8:53608/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  1:39524/8‑149 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGa                    >  5:348553/8‑146 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGa                    <  7:43814/139‑1 (MQ=11)
                                                                                                                               gccattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGATGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  7:374777/8‑149 (MQ=11)
                                                                                                                                 cattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGACCTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAat                >  1:166547/6‑148 (MQ=11)
                                                                                                                                 cattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCTATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatg               <  7:176261/144‑1 (MQ=11)
                                                                                                                                 cattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatg               <  6:31938/144‑1 (MQ=255)
                                                                                                                                 cattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAat                >  5:331280/6‑148 (MQ=255)
                                                                                                                                 cattaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGAGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAat                >  8:193900/6‑148 (MQ=11)
                                                                                                                                  attaccTACAACGCCCATCGCGACAAAACTGCCAGAGACCTGCAGCTGGCGAACGCGGCAATTACTGACATGCAGAGGCGGCAGCGGTATGTTGCTGCGCTCGATGCAAAATCCACGAAGGAGTTAGCTGATGCTAAAGCTGAAActga              <  6:25959/145‑1 (MQ=11)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              <  8:144266/145‑1 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              <  1:353360/145‑1 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              >  1:325986/5‑149 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              >  8:331530/5‑149 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              <  1:198498/145‑1 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              <  1:102958/145‑1 (MQ=255)
                                                                                                                                  attaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTCCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatga              >  1:95871/5‑149 (MQ=11)
                                                                                                                                   ttaccTACAAATCCCTGCTCGACAAAAATGCCAGACATCTGAAGCTTGCGAACGCGGCAATTACTGACTTGCAGATGCGTCAGCGTGATGTTGCTTCGGTCGATTCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgat             <  5:119225/146‑1 (MQ=11)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                                                          <  2:356326/101‑1 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgat             <  2:117504/146‑1 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgat             >  5:285025/4‑149 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgat             >  3:184011/4‑149 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgat             <  4:328385/146‑1 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAat                <  1:135013/143‑1 (MQ=255)
                                                                                                                                   ttaccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAat                >  2:135013/4‑146 (MQ=255)
                                                                                                                                    taccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctg                                                                  <  4:144018/93‑1 (MQ=255)
                                                                                                                                     accTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATacac                                               <  2:156923/112‑1 (MQ=12)
                                                                                                                                     accTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgatg            <  3:220804/147‑1 (MQ=12)
                                                                                                                                     accTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgatg            <  2:206194/147‑1 (MQ=12)
                                                                                                                                     accTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgatg            <  6:179961/147‑1 (MQ=12)
                                                                                                                                      ccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCTATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGc           <  5:207914/148‑1 (MQ=255)
                                                                                                                                      ccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGct          >  3:355673/1‑149 (MQ=12)
                                                                                                                                      ccTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGct          >  7:165398/1‑149 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAc                                                                                              <  8:12259/64‑1 (MQ=255)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGc                                   >  8:345217/1‑123 (MQ=17)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 <  8:384094/141‑1 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGaaaa                 >  5:195036/1‑141 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctc         >  7:72250/1‑149 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctc         >  3:309441/1‑149 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctc         >  2:328158/1‑149 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCCGCGCTCGATGCAAAATACACGAAGGAGTTAg                                    >  4:226166/1‑122 (MQ=12)
                                                                                                                                       cTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCCGCGCTCGATGCAAAATACACGAAGGAGTTAg                                    <  3:226166/122‑1 (MQ=12)
                                                                                                                                        tACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                                                                              <  4:5375/79‑1 (MQ=255)
                                                                                                                                        tACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTaaa                         <  6:196569/132‑1 (MQ=12)
                                                                                                                                        tACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctct        >  3:54982/1‑149 (MQ=12)
                                                                                                                                        tACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctct        >  5:307452/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGCACGCGCAAATTACTGACATGCAGATGCGTCAGCGTGATGTTTCTGCGCTCGATTCAAAATACAAGAAGGAGTTAGCTGATGCTAAAGCTGAAAAAGATGATCTg       >  4:289819/1‑149 (MQ=11)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                                                                   >  7:135925/1‑89 (MQ=255)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaaaa                                                    >  6:356214/1‑104 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCa                                                       <  2:390390/101‑1 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                                           >  2:22562/1‑113 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGc                                   <  7:285920/121‑1 (MQ=17)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTaaa                         <  6:218390/131‑1 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTg       >  8:303010/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTg       >  5:65440/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTg       >  3:103227/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTg       >  6:280063/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTg       >  7:330216/1‑149 (MQ=12)
                                                                                                                                         aCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAAGTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCa                                                       >  5:247138/1‑101 (MQ=255)
                                                                                                                                          cAAAGCCCCGCGCGACAAAAATGCCAGAGAAGTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGGTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      <  5:263512/149‑1 (MQ=11)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAatgatg            >  8:178938/1‑143 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctct        <  3:129275/147‑1 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctct        <  5:168445/147‑1 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      >  3:165507/1‑149 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      <  7:199085/149‑1 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      >  3:77407/1‑149 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      <  4:27667/149‑1 (MQ=12)
                                                                                                                                          cAAAGCCCAGCGCGACAAAAAATACATAGAACGGAAGCTTGCGAACGCGGCAATTAATTACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCTATGCAAAATACACTAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGctct        <  2:287561/147‑1 (MQ=11)
                                                                                                                                           aaaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      >  6:145761/1‑148 (MQ=12)
                                                                                                                                           aaaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc      <  5:145761/148‑1 (MQ=12)
                                                                                                                                           aaaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCg     >  7:247149/1‑149 (MQ=12)
                                                                                                                                           aaaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCg     >  8:237944/1‑149 (MQ=12)
                                                                                                                                           aaaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCg     <  3:320826/149‑1 (MQ=12)
                                                                                                                                            aaGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGt    >  2:184185/1‑149 (MQ=9)
                                                                                                                                             aGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgatg            >  2:264420/1‑140 (MQ=38)
                                                                                                                                             aGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgatg            <  1:264420/140‑1 (MQ=38)
                                                                                                                                               |                                                                                                                                                   
GTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGA  >  NC_000913/578067‑578357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: