Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,011,492 T→A 6.0% T139T (ACT→ACA yigM → putative inner membrane EamA‑like transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,4920TA 6.0% 215.9 / 4.1 67T139T (ACT→ACAyigMputative inner membrane EamA‑like transporter
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (2/2);  total (32/35)
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.11e-01

CCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGG  >  NC_000913/4011347‑4011625
                                                                                                                                                 |                                                                                                                                     
ccGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTgg                                                                                                                                     <  8:331846/148‑1 (MQ=255)
       gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTtgctg                                                                                                                              >  6:274137/1‑148 (MQ=255)
         cgccgcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGACTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTgctgct                                                                                                                            <  1:216202/148‑1 (MQ=255)
         cgccgcTCTACATCACGCTGATGTATGTCATCATGAGTAAGCGCCGTCTGAGCTGGGACTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGGGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTTGCTtgctgcg                                                                                                                            <  1:222161/148‑2 (MQ=255)
          gccgcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGctgctg                                                                                                                           <  6:49757/148‑1 (MQ=255)
            cgcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGt                                                                                                                         <  7:370597/148‑1 (MQ=255)
                 tACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCa                                                                                                                      <  6:346787/146‑1 (MQ=255)
                  aCATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACAGGCTTGCTGCTGGTGCAAct                                                                                                                   >  7:217906/1‑148 (MQ=255)
                  aCATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACAGGCTTGCTGCTGGTGCAAct                                                                                                                   >  7:223894/1‑148 (MQ=255)
                         gCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAat                                                                                                            <  3:16269/148‑1 (MQ=255)
                             aTTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATc                                                                                                         <  8:135144/147‑1 (MQ=255)
                                     catcatGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGc                                                                                                 >  5:104798/1‑147 (MQ=255)
                                        catGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCAtt                                                                                             >  7:319539/1‑148 (MQ=255)
                                           gAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTgg                                                                                           >  2:320742/1‑147 (MQ=255)
                                           gAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTgg                                                                                           >  1:277187/1‑147 (MQ=255)
                                            aGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGc                                                                                          <  5:274137/147‑1 (MQ=255)
                                              tAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCAtg                                                                                       <  2:277187/148‑1 (MQ=255)
                                              tAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCAtg                                                                                       <  8:319539/148‑1 (MQ=255)
                                               aaGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCAtg                                                                                       <  7:257533/147‑1 (MQ=255)
                                               aaGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCAtg                                                                                       >  8:257533/1‑147 (MQ=255)
                                                        cTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTAc                                                                              >  6:20766/1‑147 (MQ=255)
                                                          gcgcTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGc                                                                                                 <  5:337978/126‑1 (MQ=255)
                                                          gcgcTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGc                                                                                                 >  6:337978/1‑126 (MQ=255)
                                                          gcgcTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACaaa                                                                           >  4:325893/1‑148 (MQ=255)
                                                           cgcTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACAGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACaa                                                                            <  8:217906/146‑1 (MQ=255)
                                                           cgcTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACAGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACaa                                                                            <  8:223894/146‑1 (MQ=255)
                                                                gggCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTg                                                                     <  6:104798/148‑1 (MQ=255)
                                                                gggCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTg                                                                     <  5:20766/148‑1 (MQ=255)
                                                                   cTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATg                                                                  >  6:175064/1‑148 (MQ=255)
                                                                       gCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaaa                                                              >  5:186513/1‑148 (MQ=255)
                                                                       gCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaaa                                                              >  5:168241/1‑148 (MQ=255)
                                                                       gCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaa                                                               <  5:240798/147‑1 (MQ=255)
                                                                       gCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaa                                                               >  6:240798/1‑147 (MQ=255)
                                                                           ttAGCGCCTTGCTGTCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGa                                                                    >  2:296622/1‑138 (MQ=255)
                                                                           ttAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCa                                                                                               <  7:209933/111‑1 (MQ=255)
                                                                           ttAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCa                                                                                               >  8:209933/1‑111 (MQ=255)
                                                                           ttAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGa                                                                    <  1:296622/138‑1 (MQ=255)
                                                                            tAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACg                                                                         <  8:400853/132‑1 (MQ=255)
                                                                            tAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACg                                                                         >  7:400853/1‑132 (MQ=255)
                                                                              gcgcCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGcc                                                        <  5:175064/147‑1 (MQ=255)
                                                                              gcgcCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGcc                                                        <  3:325893/147‑1 (MQ=255)
                                                                                   ttGCTGGCGGAGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGc                                                  >  8:23332/1‑148 (MQ=255)
                                                                                         gCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAg                                             <  6:186513/147‑1 (MQ=255)
                                                                                            gTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATa                                         <  1:320742/148‑1 (MQ=255)
                                                                                                ttGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaa                                                               >  3:109725/1‑122 (MQ=255)
                                                                                                ttGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGaa                                                               <  4:109725/122‑1 (MQ=255)
                                                                                                  ggTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAg                                             <  4:73893/138‑1 (MQ=255)
                                                                                                  ggTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAg                                             >  3:73893/1‑138 (MQ=255)
                                                                                                        gggATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGt                             <  6:168241/148‑1 (MQ=255)
                                                                                                         ggATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTg                            >  5:305372/1‑148 (MQ=255)
                                                                                                           tttattCGCGATGATCAGGTGACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGt                          <  2:218949/147‑1 (MQ=255)
                                                                                                             tattCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATa                                         >  1:212128/1‑131 (MQ=255)
                                                                                                             tattCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATa                                         <  2:212128/131‑1 (MQ=255)
                                                                                                              attCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCt                                   >  3:113283/1‑136 (MQ=255)
                                                                                                              attCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCt                                   <  4:113283/136‑1 (MQ=255)
                                                                                                               ttCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATa                                         >  1:389263/1‑129 (MQ=255)
                                                                                                               ttCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATa                                         <  2:389263/129‑1 (MQ=255)
                                                                                                                tCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACaa                                                                            >  5:351116/1‑93 (MQ=255)
                                                                                                                tCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACaa                                                                            <  6:351116/93‑1 (MQ=255)
                                                                                                                              tCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCt       >  8:352077/1‑148 (MQ=255)
                                                                                                                                aCCGACCATTTCTGGACTGGCTTGCTTCTGGTGGAACTCTCCAATATCACTTTTGCCATTGGGATGGTGGGTTACAAACGCCTGTTGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGt                          <  2:3309/127‑1 (MQ=255)
                                                                                                                                 cagACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAAAATCACTATTGCCATTGGCATGGAGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACACCAAAACGCCATAGCGTGGtt                         >  1:3309/3‑127 (MQ=255)
                                                                                                                                   gACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCtggtgg  >  4:246242/1‑148 (MQ=255)
                                                                                                                                   gACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTCGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCtggtgg  >  7:232228/1‑148 (MQ=255)
                                                                                                                                    aCCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCtggtgg  >  4:175988/1‑147 (MQ=255)
                                                                                                                                           cTGGACTGGCTTGCTGCTTTTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAACCTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTATATCTTGGCGCGTTTCtg      <  7:163302/136‑1 (MQ=255)
                                                                                                                                            tGGAATGGCTTGCTGCTGGTGAAAATCTCCAATATCAATTTTTCCATTGGCATGGTCGGTTACAAACGACTGATGGAAACTAGCCCGATGCCACAGCATAACGACTTAGCGTGGTTCTATCTTGGCGCGTTTCtgg     >  8:163302/1‑136 (MQ=255)
                                                                                                                                                 |                                                                                                                                     
CCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGGACTGGCTTGCTGCTGGTGCAACTCTCCAATATCACTTTTGCCATTGGCATGGTGGGTTACAAACGCCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGG  >  NC_000913/4011347‑4011625

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