Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,716,846 T→A 5.5% N126I (AAC→ATC)  tyrS ← tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,716,8460TA 5.5% 238.8 / 3.5 73N126I (AAC→ATC) tyrStyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TAGGAAAACTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATCAGAC  >  NC_000913/1716699‑1716993
                                                                                                                                                   |                                                                                                                                                   
tAGGAAAACTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAgtt                                                                                                                                                    >  2:248437/1‑149 (MQ=255)
tAGGAAAACTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAgt                                                                                                                                                     <  7:68785/148‑1 (MQ=255)
   gAAAACTCAGTGAATGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAgttgtt                                                                                                                                                 <  3:246250/149‑1 (MQ=255)
   gAAAACTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAgttgtt                                                                                                                                                 >  5:141932/1‑149 (MQ=255)
     aaaCTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTcg                                                                                                                                               <  1:60285/149‑1 (MQ=255)
         tCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGATGTTcgccgc                                                                                                                                           >  7:402381/1‑149 (MQ=255)
         tCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGATGTTcgccgc                                                                                                                                           >  7:396996/1‑149 (MQ=255)
             tGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATa                                                                                                                                       <  2:306197/149‑1 (MQ=255)
                aCGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAg                                                                                                                                      >  5:149940/1‑147 (MQ=255)
                  gAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCaga                                                                                                                                  <  7:204409/149‑1 (MQ=255)
                        cccTGATCTTCACGGTTGAGCCGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTtc                                                                                                                            <  8:294527/149‑1 (MQ=255)
                        cccTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGtttt                                                                                                                             <  7:284937/148‑1 (MQ=255)
                          cTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGACAATATCGAGCAGGAAGGTCAGCACATTCATATTGCCGAACAAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTtctc                                                                                                                          >  6:319414/1‑149 (MQ=255)
                           tGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGCAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTtctc                                                                                                                          >  3:335648/1‑148 (MQ=255)
                           cgATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTcc                                                                                                                         >  6:83893/2‑149 (MQ=255)
                             aTCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCca                                                                                                                        <  2:197416/148‑1 (MQ=255)
                                ttCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACACTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGt                                                                                                                    >  7:122402/1‑149 (MQ=255)
                                  cACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTc                                                                                                                   >  5:114674/1‑148 (MQ=255)
                                       ttGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGATGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATc                                                                                                             <  8:396996/149‑1 (MQ=255)
                                       ttGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGATGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATc                                                                                                             <  8:402381/149‑1 (MQ=255)
                                                     aaCCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCaa                                                                                               >  2:268269/1‑149 (MQ=255)
                                                     aaCCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCaa                                                                                               >  8:333402/1‑149 (MQ=255)
                                                      aCCGCTTCTTTGTTGATCATATTGTTAAAGGAGAAGTGTTTGCCAATATCGCTCAGGAAGGTCAGCACATTCATATTGCAGAACAAGTGATAGTTTTTCGCCGCGATAGCAGAGTTTTCTCAACAGTAGAAATCGATTAACGGGGAAAc                                                                                              <  1:201039/149‑1 (MQ=255)
                                                       ccGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAAc                                                                                              >  1:130730/1‑148 (MQ=255)
                                                       ccGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAAc                                                                                              <  2:130730/148‑1 (MQ=255)
                                                        cGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCt                                                                                            >  7:236238/1‑149 (MQ=255)
                                                        cGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCt                                                                                            >  4:96875/1‑149 (MQ=255)
                                                         gCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTg                                                                                           <  7:333402/149‑1 (MQ=255)
                                                           ttcttTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGtttt                                                                                                                             >  1:11361/1‑113 (MQ=255)
                                                           ttcttTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGtttt                                                                                                                             <  2:11361/113‑1 (MQ=255)
                                                             cttTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTa                                                                                       >  4:149781/1‑149 (MQ=255)
                                                             cttTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGCAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCtt                                                                                        <  4:335648/148‑1 (MQ=255)
                                                                tgttgATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACg                                                                                     >  3:230927/1‑148 (MQ=255)
                                                                  ttgATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGa                                                                                   <  6:114674/148‑1 (MQ=255)
                                                                  ttgATCATCTGGTTAACCGAGAAGTGTTTGGCAATATCGCGCAGGAAGGTCAGCTCATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCCCAGTCGTAATAGAGGAACGGGGCAACCTGCTGACGGAt                                                                                  >  8:48977/1‑149 (MQ=255)
                                                                    gATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGAttt                                                                                >  2:257129/1‑149 (MQ=255)
                                                                     atcatcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGAtttt                                                                               >  3:28321/1‑149 (MQ=255)
                                                                     atcatcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGAtttt                                                                               <  1:268269/149‑1 (MQ=255)
                                                                     atcatcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATAGAGCCGAACCAGTCATAGTTGTTCGACTCGATAGCAGAGTTTTCTCCACATTCGAAATCGAGGAACGGGGCAACCTGCTTACGGAtttt                                                                               <  5:319414/149‑1 (MQ=255)
                                                                        atcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGt                                                                             <  3:255353/148‑1 (MQ=255)
                                                                        atcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTc                                                                            <  1:211463/149‑1 (MQ=255)
                                                                        atcTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTc                                                                            >  2:211463/1‑149 (MQ=255)
                                                                              ttAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTccaccca                                                                      >  6:54334/1‑149 (MQ=255)
                                                                                  cGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACAcc                                                                  <  5:83893/149‑1 (MQ=255)
                                                                                       aaGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAAc                                                             <  7:30259/149‑1 (MQ=255)
                                                                                        aGTGTTTGCCAATATCGCGGAGTAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACa                                                            <  4:360523/149‑1 (MQ=255)
                                                                                        aGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACa                                                            <  1:257129/149‑1 (MQ=255)
                                                                                            tttGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGttt                                                        >  5:141088/1‑149 (MQ=255)
                                                                                             ttGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTccacc                                                                        <  1:8532/132‑1 (MQ=255)
                                                                                             ttGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTccacc                                                                        >  2:8532/1‑132 (MQ=255)
                                                                                               gCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTttctt                                                     >  1:217324/1‑149 (MQ=255)
                                                                                                ccAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTttctt                                                     <  7:30808/148‑1 (MQ=255)
                                                                                                 cAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCg                                                   <  8:236238/149‑1 (MQ=255)
                                                                                                    tatCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGgtg                                                >  3:363883/1‑149 (MQ=255)
                                                                                                        gcgcAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGtt                                              >  3:88716/1‑147 (MQ=255)
                                                                                                          gcaggcagGTCAGCACATTCATATTGCCGAACCAGTCATACTTGTTCGCCGCGATAGCAGAGCTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGCACAGTTTCTTCGGTGTTCAGc                                          <  4:249013/149‑1 (MQ=255)
                                                                                                            aGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCt                                         >  6:168710/1‑148 (MQ=255)
                                                                                                                   tCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTc                                                                            <  8:153612/106‑1 (MQ=255)
                                                                                                                   tCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTc                                                                            >  7:153612/1‑106 (MQ=255)
                                                                                                                     aGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCgg                                <  6:149940/148‑1 (MQ=255)
                                                                                                                       cacaTTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAg                             <  3:96875/149‑1 (MQ=255)
                                                                                                                              atatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTg                        >  3:183833/1‑147 (MQ=255)
                                                                                                                              atatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTg                        >  2:295554/1‑147 (MQ=255)
                                                                                                                              atatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGa                       <  4:88716/148‑1 (MQ=255)
                                                                                                                               tatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAg                     <  4:28321/149‑1 (MQ=255)
                                                                                                                               tatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAg                     >  7:241860/1‑149 (MQ=255)
                                                                                                                               tatTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTGTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAg                     >  8:114346/1‑149 (MQ=255)
                                                                                                                                atTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGc                                                                                                 <  3:154704/72‑1 (MQ=255)
                                                                                                                                 ttGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGc                    <  8:122402/148‑1 (MQ=255)
                                                                                                                                   gCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTc                  <  6:141932/148‑1 (MQ=255)
                                                                                                                                   gCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTc                  <  5:168710/148‑1 (MQ=255)
                                                                                                                                   gCCGAACCAGTAATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTc                  <  3:12853/148‑1 (MQ=255)
                                                                                                                                     cGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGTCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTGCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCggg               <  6:115547/149‑1 (MQ=255)
                                                                                                                                          cAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGc           >  2:190625/1‑148 (MQ=255)
                                                                                                                                             tCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAAt       >  5:323803/1‑149 (MQ=255)
                                                                                                                                               aTAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAAt       >  8:333294/1‑147 (MQ=255)
                                                                                                                                               aTAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATc      <  5:54334/148‑1 (MQ=255)
                                                                                                                                                  gttgttCGCAGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCATACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATCAGAc  >  8:111314/1‑149 (MQ=255)
                                                                                                                                                   ttgttCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATCAGAc  <  4:356384/148‑1 (MQ=255)
                                                                                                                                                   ttgttCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATCAGAc  >  3:356384/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGGAAAACTCAGTGAACGAAATCCCCTGATCTTCACGGTTGAGACGCTGCTTAACCGCTTCTTTGTTGATCATCTGGTTAACGGAGAAGTGTTTGCCAATATCGCGCAGGAAGGTCAGCACATTCATATTGCCGAACCAGTCATAGTTGTTCGCCGCGATAGCAGAGTTTTCTCCACAGTCGAAATCGAGGAACGGGGCAACCTGCTTACGGATTTTGTCCACCCACTCCTGAACAGTTTCTTCGGTGTTCAGCTTACGCTCGGCAGCTTTGAAGCTCGGGTCGCCAATCAGAC  >  NC_000913/1716699‑1716993

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