Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 218,786 T→C 6.2% intergenic (‑11/+101) proS ← / ← tsaA prolyl‑tRNA synthetase/tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913218,7860TC 6.2% 186.3 / 2.9 64intergenic (‑11/+101)proS/tsaAprolyl‑tRNA synthetase/tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base C (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACG  >  NC_000913/218640‑218932
                                                                                                                                                  |                                                                                                                                                  
cGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTc                                                                                                                                                  <  3:192850/149‑1 (MQ=255)
  gTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCa                                                                                                                                                <  6:158915/149‑1 (MQ=255)
                atataACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGaa                                                                                                                                  >  6:10606/1‑149 (MQ=255)
                 tataACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGCTCCAGTTGGAAGGTAGaa                                                                                                                                  >  6:127427/1‑148 (MQ=255)
                 tataACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGCTCCAGTTGGAAGGTAGaa                                                                                                                                  >  6:130681/1‑148 (MQ=255)
                       ccGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCt                                                                                                                            >  5:187917/1‑148 (MQ=255)
                       ccGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCt                                                                                                                            >  4:70042/1‑148 (MQ=255)
                                     gCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTg                                                                                                               <  3:63260/147‑1 (MQ=255)
                                       ggATCAGCCCGGCGCGCAGCATCAGCTGATGGCGGATCACCTCGGCGTCGGCAGGTGGCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAc                                                                                                           <  6:170669/149‑1 (MQ=255)
                                        gATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAc                                                                                                           >  5:60605/1‑148 (MQ=255)
                                         atcatcCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGTCTCAAGGCGAGCCTGGGACa                                                                                                          >  6:90704/1‑148 (MQ=255)
                                               ccGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaaaa                                                                                                   <  4:39234/149‑1 (MQ=255)
                                                cGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAg                                                                                                  <  6:187917/149‑1 (MQ=255)
                                                 ggCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGt                                                                                                 <  3:70042/149‑1 (MQ=255)
                                                  gcgcgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGt                                                                                                 >  1:116753/1‑148 (MQ=255)
                                                    gcgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGAt                                                                                              <  5:153693/149‑1 (MQ=255)
                                                            tCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTa                                                                                      >  7:126876/1‑149 (MQ=255)
                                                            tCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTa                                                                                      >  5:27776/1‑149 (MQ=255)
                                                                   aTGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTCATTTAGTTTACCAGTGc                                                                               <  2:116753/149‑1 (MQ=255)
                                                                     ggCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTATGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCa                                                                              >  2:11245/1‑148 (MQ=255)
                                                                        tGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCaa                                                                             >  3:215834/1‑146 (MQ=255)
                                                                        tGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCaa                                                                             >  1:115627/1‑146 (MQ=255)
                                                                          aTCACCTCGTCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTAATGGCTAGTAAGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGTAAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGa                                                                         <  1:46056/148‑1 (MQ=255)
                                                                           tCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGtt                                                                                        >  8:103493/1‑132 (MQ=255)
                                                                           tCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGtt                                                                                        <  7:103493/132‑1 (MQ=255)
                                                                             aCCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATg                                                                     <  7:239046/149‑1 (MQ=255)
                                                                               cTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGt                                                                    >  4:135634/1‑148 (MQ=255)
                                                                                  ggCGTCGGCAGGTGGCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCaaa                                                                <  7:58062/149‑1 (MQ=255)
                                                                                   gCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCaa                                                                 >  4:232983/1‑147 (MQ=255)
                                                                                        ggCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAagag                                                            >  6:174463/1‑147 (MQ=255)
                                                                                          cAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAgg                                                        <  8:126876/149‑1 (MQ=255)
                                                                                          cAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAAGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAgg                                                        <  5:10606/149‑1 (MQ=255)
                                                                                              tgtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGt                                                    <  4:215834/149‑1 (MQ=255)
                                                                                               gtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAg                                                                                          >  1:233982/1‑110 (MQ=255)
                                                                                               gtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAg                                                                                          <  2:233982/110‑1 (MQ=255)
                                                                                                tctcCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTg                                                   >  8:216067/1‑148 (MQ=255)
                                                                                                  tcCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAAt                                                >  7:183348/1‑149 (MQ=255)
                                                                                                    cTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGTTTCCATTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAAttt                                              <  5:90704/149‑1 (MQ=255)
                                                                                                           gTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCgg                                       <  4:196114/149‑1 (MQ=255)
                                                                                                              gagaGCAGGTAATGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTCGGACAAAAAAAAGTGATTTAGTTTACCAGGGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTc                                    >  4:163223/1‑149 (MQ=255)
                                                                                                                  gCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAg                                 <  7:216067/148‑1 (MQ=255)
                                                                                                                    aGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGt                                                    <  4:231795/127‑1 (MQ=255)
                                                                                                                    aGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGt                                                    >  3:231795/1‑127 (MQ=255)
                                                                                                                    aGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAgcg                               >  3:91153/1‑148 (MQ=255)
                                                                                                                      gTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCa                             >  7:10672/1‑148 (MQ=255)
                                                                                                                       tATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaaa                           <  3:85335/149‑1 (MQ=255)
                                                                                                                       tATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaaa                           <  8:10672/149‑1 (MQ=255)
                                                                                                                         ttGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGa                         >  1:73513/1‑149 (MQ=255)
                                                                                                                           ggCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGt                                                    <  2:69364/120‑1 (MQ=255)
                                                                                                                           ggCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGt                                                    >  1:69364/1‑120 (MQ=255)
                                                                                                                             cTAGTACGCATGTTGTTACGGCTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTc                     <  5:130681/149‑1 (MQ=255)
                                                                                                                             cTAGTACGCATGTTGTTACGGCTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTc                     <  5:127427/149‑1 (MQ=255)
                                                                                                                                gTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAAt                                                                      >  8:11851/1‑97 (MQ=255)
                                                                                                                                gTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAAt                                                                      <  7:11851/97‑1 (MQ=255)
                                                                                                                                gTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGCGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCaa                   >  4:159511/1‑148 (MQ=255)
                                                                                                                                 tACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCaaaa                 >  8:237207/1‑149 (MQ=255)
                                                                                                                                    gCATGTTGTTATGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCt              <  1:11245/149‑1 (MQ=255)
                                                                                                                                    gCATGTTGTTACGGTTCCTGTTGGTAGGTAGAACAGGCTCATGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCTAAAGAAATGTCAAAAGAGAAGGGCGTGTATTTAACGCGGTTCCAgcgc                              <  2:105035/133‑1 (MQ=255)
                                                                                                                                    gCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCt              >  6:15183/1‑149 (MQ=255)
                                                                                                                                    gCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGATTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAgcgc                              >  1:105035/1‑133 (MQ=255)
                                                                                                                                      aTGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCGCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGc            <  2:52999/149‑1 (MQ=255)
                                                                                                                                            ttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGt      >  1:56390/1‑149 (MQ=255)
                                                                                                                                              aCGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGa    >  1:138360/1‑149 (MQ=255)
                                                                                                                                                ggTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGAcg  >  3:106360/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCCGGCGCGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACG  >  NC_000913/218640‑218932

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: