Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 542,182 G→A 6.2% V99I (GTC→ATC)  glxK → glycerate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913542,1820GA 6.2% 208.9 / 2.4 68V99I (GTC→ATC) glxKglycerate kinase II
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGC  >  NC_000913/542035‑542326
                                                                                                                                                   |                                                                                                                                                
gATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTgg                                                                                                                                                  <  5:358009/148‑1 (MQ=255)
   gCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTcg                                                                                                                                               >  4:11480/1‑148 (MQ=255)
      aagctCGCCGCGACTGGCGTCAACATCGTGACGCTTGAATACTTCGGGCCGATGGGCGAAAAAGTGAATGCTTTAAAAGGCCTAACCGGCGACGGGAAAACGGCGGCGATTGAGATGGCGGCACCAAGTGGCCTGATGCTGCCCgcgc                                                                                                                                            <  8:198101/144‑1 (MQ=255)
      aTGGTCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCgcgc                                                                                                                                            >  7:198101/1‑148 (MQ=255)
          tCGCCGCGCCGCGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGCAAAACGGCGGGGATTGAGATGGCGGCAGCAAGTGGCCTGATGCAGGTCGCGccag                                                                                                                                         >  2:284085/1‑145 (MQ=255)
          tCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTGTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTg                                                                                                                                         >  3:12930/1‑147 (MQ=255)
            gccgcGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGaaa                                                                                                                                      <  1:237897/148‑1 (MQ=255)
              cgcgACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAg                                                                                                                                    <  3:181358/148‑1 (MQ=255)
                 gACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCg                                                                                                                                  >  7:372764/1‑147 (MQ=255)
                 gACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGt                                                                                                                                 <  2:305902/148‑1 (MQ=255)
                 gACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGt                                                                                                                                 <  5:243862/148‑1 (MQ=255)
                   cGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTaa                                                                                                                               <  4:334800/148‑1 (MQ=255)
                         gcAACATCGTGACGCTTGAAGTCTGCGGGCCGAAGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGtt                                                                                                                         <  7:255446/148‑1 (MQ=255)
                                 gTGACGCTTGAAGGCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGTTTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGcc                                                                                                                  <  3:35727/147‑1 (MQ=255)
                                    aCGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTc                                                                                                                >  4:144799/1‑146 (MQ=255)
                                      gCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGt                                                                                                            <  2:271539/148‑1 (MQ=255)
                                         tGAAGTCTGCGGGCCGATGGGTGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGATCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTGCAGtttt                                                                                                         >  7:283237/1‑148 (MQ=255)
                                         tGAAGTCTGCGGGCCGATGGGTGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGATCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGtttt                                                                                                         >  7:276818/1‑148 (MQ=255)
                                          gAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGtttt                                                                                                         <  3:101621/147‑1 (MQ=255)
                                          gAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGttt                                                                                                          <  7:231587/146‑1 (MQ=255)
                                             gTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTa                                                                                                     <  3:144799/148‑1 (MQ=255)
                                                 gCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGTGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCa                                                                                                              >  7:342131/1‑135 (MQ=255)
                                                 gCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGTGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCa                                                                                                              <  8:342131/135‑1 (MQ=255)
                                                 gCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACggg                                                                                                 >  4:244036/1‑148 (MQ=255)
                                                  cGGGCCCATGGGCGAAAAACGGCATGCTTTTTATGGCCTTACCGGCGACGGGAAAACCGCGGGGATTGACCTGGCGGCACCAATTGGCCCGCAGCCGGGCGCGCCTGAAAAGCCTAATCCGTTGCTGGCCTCCACTTTTGGTACggg                                                                                                 <  8:7257/147‑1 (MQ=255)
                                                   gggCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACggggg                                                                                               <  3:244036/148‑1 (MQ=255)
                                                        gATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGtt                                                                                           >  6:315572/1‑147 (MQ=255)
                                                        gATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTa                                                                                          <  3:11480/148‑1 (MQ=255)
                                                          tGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAAt                                                                                        <  4:169785/148‑1 (MQ=255)
                                                            ggCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAAtt                                                                                       <  5:198734/147‑1 (MQ=255)
                                                            ggCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTc                                                                                      <  8:372764/148‑1 (MQ=255)
                                                             gCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCg                                                                                     <  6:33414/148‑1 (MQ=255)
                                                              cGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGt                                                                                    <  2:113193/148‑1 (MQ=255)
                                                              cGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGt                                                                                    <  4:73780/148‑1 (MQ=255)
                                                                aaaaaGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTTGTCGCGCCTGAAAAGCGt                                                                                                                                 <  5:54881/101‑1 (MQ=255)
                                                                aaaaaGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGt                                                                                                                                 >  6:54881/1‑101 (MQ=255)
                                                                aaaaaGTGAAAGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGCTGATTGAGATGGAGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGATGGCCTACAGTATTGGTACGGTGGAGTTAATTCGTCa                                                                                  >  4:255978/1‑148 (MQ=255)
                                                                 aaaaGTGAAAGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGCGATGGCGGCAGCAAGTGGCCTGATGCTGCTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCAt                                                                                 >  8:46631/1‑148 (MQ=255)
                                                                    aGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATgcg                                                                              <  2:136453/148‑1 (MQ=255)
                                                                            cTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGATGGTCGCGCCTGAAAAGCGTAATCCTTTGCTTGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATa                                                                       <  3:255978/147‑1 (MQ=255)
                                                                               tttATGGCCTTACGGGCGACGGGTAAACGGCGGTGTTTGAGTTGGCGGCAGCAAGTGGCCTGTTGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACg                                                                    <  3:305595/147‑1 (MQ=255)
                                                                                ttATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGa                                                                   <  6:82748/147‑1 (MQ=255)
                                                                                      ccTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGCTGCTGCTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGc                                                            >  3:315656/1‑148 (MQ=255)
                                                                                         tACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCa                                                          >  8:410660/1‑147 (MQ=255)
                                                                                          aCCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCata                                                        >  3:73319/1‑148 (MQ=255)
                                                                                             ggCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGAc                                                                  <  2:254839/135‑1 (MQ=255)
                                                                                             ggCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGAc                                                                  >  1:254839/1‑135 (MQ=255)
                                                                                                 aCGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGATCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTTCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtct                                                 <  8:283237/148‑1 (MQ=255)
                                                                                                 aCGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGATCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtct                                                 <  8:276818/148‑1 (MQ=255)
                                                                                                  cGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtctc                                                >  4:126348/1‑148 (MQ=255)
                                                                                                  cGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCGGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtctc                                                <  1:284085/148‑1 (MQ=255)
                                                                                                    ggAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCgg                                              >  4:33928/1‑148 (MQ=255)
                                                                                                        aaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCAtt                                          >  5:273443/1‑148 (MQ=255)
                                                                                                                 gATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGt                                 >  6:345257/1‑148 (MQ=255)
                                                                                                                   ttGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGc                               >  1:181349/1‑148 (MQ=255)
                                                                                                                                agcaAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcggc                  <  2:181349/148‑1 (MQ=255)
                                                                                                                                   aaGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTACGCCGTTGCTGGCCTCCCGTTTTGGTACGGGGGAGGTACTTCGTCATGCGCTGGATAAGGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGt               <  4:12930/148‑1 (MQ=255)
                                                                                                                                   aaGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGt               <  6:273443/148‑1 (MQ=255)
                                                                                                                                    aGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAg                                                                                             <  7:307996/69‑1 (MQ=255)
                                                                                                                                    aGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAg                                                                                             >  8:307996/1‑69 (MQ=255)
                                                                                                                                     gTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAt             >  5:299204/1‑148 (MQ=255)
                                                                                                                                     gTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGCTAt             >  1:101612/1‑148 (MQ=255)
                                                                                                                                           tGATGCTGGTCGCGCCTGAAACGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCAATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCAt       >  1:209731/1‑148 (MQ=255)
                                                                                                                                            gATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCAtatt                                                      >  6:295166/1‑100 (MQ=255)
                                                                                                                                            gATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCAtatt                                                      <  5:295166/100‑1 (MQ=255)
                                                                                                                                            gATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGCCGGCGGTATGGGCATg      >  6:191202/1‑148 (MQ=255)
                                                                                                                                             aTGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATgg     >  3:249231/1‑148 (MQ=255)
                                                                                                                                              tGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGgc    >  3:357917/1‑148 (MQ=255)
                                                                                                                                                cTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGgcgc  <  4:357917/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGCTTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGC  >  NC_000913/542035‑542326

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