Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,514,531 T→C 6.0% L234P (CTG→CCG)  ydcU → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,514,5310TC 6.0% 201.6 / 3.2 66L234P (CTG→CCG) ydcUputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base C (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CGCTTTCTGGTGTTTCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAGGGGGCGATTGG  >  NC_000913/1514389‑1514672
                                                                                                                                              |                                                                                                                                             
cgcTTTCTGGTGTTTCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCt                                                                                                                                         >  6:356890/1‑149 (MQ=255)
          tgtTTCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCg                                                                                                                               <  6:307074/149‑1 (MQ=255)
           gtTTCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCgg                                                                                                                              >  5:161538/1‑149 (MQ=255)
            tttCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCggg                                                                                                                             <  3:64913/149‑1 (MQ=255)
               ctctATATCTGGTTGCCGTTCATGATCCTGCGCGTTCAGGCGGCGCTTGAGCGTTTGCCGGCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCGCGGGTa                                                                                                                           <  8:278440/148‑1 (MQ=255)
                tctATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTAtt                                                                                                                         >  4:4374/1‑149 (MQ=255)
                  tataTCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGc                                                                                                                       >  3:118528/1‑149 (MQ=255)
                     atCTTGTTGCCGTTCATGATGATGCCCGTTCAGGCGGCGCTTGAGCGGTTGCCGCCGTCATTGTTGCAGGCGTCTGCTGATCTCGGCGCACGTCGACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCgg                                                                                                                              <  2:280753/139‑1 (MQ=255)
                     atCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCgg                                                                                                                              >  1:280753/1‑139 (MQ=255)
                           ttGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGctc                                                                                                              >  6:324567/1‑149 (MQ=255)
                               cGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATc                                                                                                          >  7:347260/1‑149 (MQ=255)
                                gTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTAACt                                                                                                         >  4:182391/1‑149 (MQ=255)
                                     tGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTa                                                                                                            <  1:342180/141‑1 (MQ=255)
                                     tGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTa                                                                                                            >  2:342180/1‑141 (MQ=255)
                                      gATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCt                                                                                                   >  3:272199/1‑149 (MQ=255)
                                        tCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTtc                                                                                                 >  7:78805/1‑149 (MQ=255)
                                        tCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTtc                                                                                                 <  1:282646/149‑1 (MQ=255)
                                          cTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTtctc                                                                                               <  4:118528/149‑1 (MQ=255)
                                             ctcGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTctcact                                                                                            >  4:62288/3‑149 (MQ=255)
                                                    aGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACt                                                                                      <  6:161538/148‑1 (MQ=255)
                                                    aGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTg                                                                                     >  2:59593/1‑149 (MQ=255)
                                                    aGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTg                                                                                     >  2:240615/1‑149 (MQ=255)
                                                     ggcggcGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTgg                                                                                    <  5:324567/149‑1 (MQ=255)
                                                       cggcgCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCTCTCACACTGGGc                                                                                  <  3:182391/149‑1 (MQ=255)
                                                        ggcgCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGc                                                                                  >  2:212983/1‑148 (MQ=255)
                                                                      tGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTccc                                                                    <  3:4374/148‑1 (MQ=255)
                                                                       gccgccGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTgg                                                                                    >  2:309924/1‑131 (MQ=255)
                                                                       gccgccGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTgg                                                                                    <  1:309924/131‑1 (MQ=255)
                                                                          gccgTCATTGTTGCAGGCGTCGGCTGCTCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGc                                                               >  5:304491/1‑149 (MQ=255)
                                                                             gTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtgg                                                            >  3:15637/1‑149 (MQ=255)
                                                                             gTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtgg                                                            <  1:212983/149‑1 (MQ=255)
                                                                              tCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtggt                                                           >  5:331882/1‑149 (MQ=255)
                                                                                aTTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtggt                                                           <  3:199404/147‑1 (MQ=255)
                                                                                aTTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtggt                                                           >  4:199404/1‑147 (MQ=255)
                                                                                       cAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGtcctc                                                   <  1:59593/148‑1 (MQ=255)
                                                                                                gCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATAttt                                         >  7:330329/1‑149 (MQ=255)
                                                                                                    aTCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTAt                                      <  4:272199/148‑1 (MQ=255)
                                                                                                     tctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCg                                    >  7:98894/1‑149 (MQ=255)
                                                                                                     tctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCCGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGatat                                           <  4:225686/142‑1 (MQ=255)
                                                                                                     tctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCCGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGatat                                           >  1:279024/1‑142 (MQ=255)
                                                                                                     tctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCCGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGatat                                           <  2:279024/142‑1 (MQ=255)
                                                                                                     tctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCCGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGatat                                           >  3:225686/1‑142 (MQ=255)
                                                                                                      ctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtggt                                                           >  4:69411/1‑125 (MQ=255)
                                                                                                      ctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCtggt                                                           <  3:69411/125‑1 (MQ=255)
                                                                                                      ctcGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCgg                                   >  5:197510/1‑149 (MQ=255)
                                                                                                         ggCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATc                                     <  8:78805/144‑1 (MQ=255)
                                                                                                         ggCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCaa                                >  3:277277/1‑149 (MQ=255)
                                                                                                          gcgcACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAat                               <  5:356890/149‑1 (MQ=255)
                                                                                                               cGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGt                          <  4:15637/149‑1 (MQ=255)
                                                                                                                gTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGtt                         <  6:197510/149‑1 (MQ=255)
                                                                                                                   cACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTAt                      <  6:304491/149‑1 (MQ=255)
                                                                                                                     cGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTc                    <  7:316809/149‑1 (MQ=255)
                                                                                                                      gACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTcc                   >  1:163034/1‑149 (MQ=255)
                                                                                                                      gACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTcc                   <  8:347260/149‑1 (MQ=255)
                                                                                                                       aCAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTccc                  <  4:277277/149‑1 (MQ=255)
                                                                                                                         aaaCCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATAtt                                          <  5:94806/123‑1 (MQ=255)
                                                                                                                         aaaCCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATAtt                                          >  6:94806/1‑123 (MQ=255)
                                                                                                                           aCCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCcag                >  1:222340/1‑147 (MQ=255)
                                                                                                                                tCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAggg           <  3:62288/147‑1 (MQ=255)
                                                                                                                                 cGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAggggg         <  2:222340/148‑1 (MQ=255)
                                                                                                                                   cTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGt                          >  4:5129/1‑129 (MQ=255)
                                                                                                                                    tATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAGGGGGCGAt     >  6:358947/1‑149 (MQ=255)
                                                                                                                                      tgtgGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAGGGGGCGATTg   <  2:163034/149‑1 (MQ=255)
                                                                                                                                        tggtgCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAGGGGGCGATTgg  <  6:331882/148‑1 (MQ=255)
                                                                                                                                         ggtgCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGtt                         <  8:16816/124‑1 (MQ=255)
                                                                                                                                         ggtgCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGtt                         >  7:16816/1‑124 (MQ=255)
                                                                                                                                              |                                                                                                                                             
CGCTTTCTGGTGTTTCTCTATATCTGGTTGCCGTTCATGATCCTGCCCGTTCAGGCGGCGCTTGAGCGTTTGCCGCCGTCATTGTTGCAGGCGTCGGCTGATCTCGGCGCACGTCCACGACAAACCTTTCGCTATGTGGTGCTGCCGCTGGCAATCCCGGGTATTGCCGCTGGCTCTATCTTTACCTTCTCACTCACACTGGGCGATTTTATCGTCCCGCAGCTGGTTGGTCCTCCAGGATATTTTATCGGCAATATGGTTTATTCCCAGCAGGGGGCGATTGG  >  NC_000913/1514389‑1514672

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: ‑