Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,676,163 C→G 8.9% intergenic (+25/‑35) hom → / → thrC homoserine dehydrogenase/threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,676,1630CG 8.9% 166.2 / 7.2 56intergenic (+25/‑35)hom/thrChomoserine dehydrogenase/threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base G (3/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

CCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTT  >  NC_002947/1676016‑1676310
                                                                                                                                                   |                                                                                                                                                   
ccAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcg                                                                                                                                                    >  5:277066/1‑149 (MQ=255)
 cAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcgc                                                                                                                                                   <  5:78453/149‑1 (MQ=255)
  aaTGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCcgcgcg                                                                                                                                                  >  3:39133/1‑149 (MQ=255)
       tCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggcccc                                                                                                                                             <  6:106829/149‑1 (MQ=255)
         cTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCGCGGGcggcccgtc                                                                                                                                           <  3:121589/149‑4 (MQ=255)
            gTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGTTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggccccggc                                                                                                                                          <  6:284557/147‑1 (MQ=255)
            gTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGcggccccggc                                                                                                                                          <  6:277066/147‑1 (MQ=255)
               ccccAAGGCGGGGGCGAGCCGCGTATCAATGACGCCATCGTCGCCCTGGACGCCCTGCAGGACGTGGTCGGCACGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCa                                                                                                                                      <  8:50048/147‑1 (MQ=255)
               aCCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCa                                                                                                                                      <  2:12914/148‑1 (MQ=255)
                  cATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCAt                                                                                                                                  <  4:39133/149‑1 (MQ=255)
                   aTGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATc                                                                                                                                 <  3:9244/149‑1 (MQ=255)
                   aTGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATc                                                                                                                                 >  4:269140/1‑149 (MQ=255)
                      gCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                              <  5:739/149‑1 (MQ=255)
                      gCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                              <  5:136710/149‑1 (MQ=255)
                          ggTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGtt                                                                                                                          <  5:277128/149‑1 (MQ=255)
                           gTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGttt                                                                                                                         <  3:265922/149‑1 (MQ=255)
                                gcagcGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCAcc                                                                                                                    >  3:207472/1‑149 (MQ=255)
                                       aTCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATgcgc                                                                                                             >  5:15376/1‑149 (MQ=255)
                                       aTCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATgcgc                                                                                                             >  3:303169/1‑149 (MQ=255)
                                          aaTGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACCTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCAACCAAGGTTTGCACCCATgcccaa                                                                                                           >  7:57568/1‑144 (MQ=255)
                                           aTGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCtata                                                                                                         <  5:180525/149‑1 (MQ=255)
                                                   aTCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTAcc                                                                                                 >  3:122733/1‑149 (MQ=255)
                                                             tGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                      <  4:201343/134‑1 (MQ=255)
                                                             tGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                      >  3:201343/1‑134 (MQ=255)
                                                              ggAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGc                                                                                      >  3:8048/1‑149 (MQ=255)
                                                               gAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAg                                                                                        <  3:269140/146‑1 (MQ=255)
                                                                aaGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCAc                                                                                    <  3:161119/149‑1 (MQ=255)
                                                                    ccTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGc                                                                                >  2:239504/1‑149 (MQ=255)
                                                                     cTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCa                                                                                                                                   <  7:122247/97‑1 (MQ=255)
                                                                     cTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTCCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCa                                                                                                                                   >  8:122247/1‑97 (MQ=255)
                                                                           gACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAAt                                                                         >  1:143887/1‑149 (MQ=255)
                                                                              gTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTc                                                                      <  4:122733/149‑1 (MQ=255)
                                                                                 gTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGaa                                                                   <  4:207472/149‑1 (MQ=255)
                                                                                 gTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGaa                                                                   <  6:15376/149‑1 (MQ=255)
                                                                                 gTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGa                                                                    <  5:233648/148‑1 (MQ=255)
                                                                                    ggCAAGGTTGTGCGCATTCGCTTCGATCAGCTCAACTAACAAATTTTGCCGACGGGCCGCCCGCGGGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCtat                                                                                                          <  6:268792/107‑1 (MQ=255)
                                                                                    ggCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATAGAAGGTTTGAACCCATGCGCtat                                                                                                          >  5:268792/1‑107 (MQ=255)
                                                                                      cAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGGACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGt                                                              >  2:272028/1‑149 (MQ=255)
                                                                                      cAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGt                                                              >  3:203571/1‑149 (MQ=255)
                                                                                          gTTGGGCGCATTCGGGTCTAGCAGACAAATTAAAAAAAATTGCCTCGGGGCCGCCGGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgctg                                                          <  2:143887/149‑1 (MQ=255)
                                                                                              tgCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTgctgg                                                      >  7:196913/1‑149 (MQ=255)
                                                                                                cgcATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGaa                                                                          >  1:29867/1‑127 (MQ=255)
                                                                                                cgcATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGaa                                                                          <  2:29867/127‑1 (MQ=255)
                                                                                                 gcATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTgctggctg                                                   >  4:243788/1‑149 (MQ=255)
                                                                                                       gcgTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGGACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGctggcctg                                              <  1:272028/148‑1 (MQ=255)
                                                                                                           cGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAg                                         >  7:212341/1‑149 (MQ=255)
                                                                                                                 gcTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGAcgg                                   >  7:105898/1‑149 (MQ=255)
                                                                                                                   tCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCg                                                                     >  1:16029/1‑113 (MQ=255)
                                                                                                                   tCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCg                                                                     <  2:16029/113‑1 (MQ=255)
                                                                                                                   tCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGAcggcg                                 >  3:141236/1‑149 (MQ=255)
                                                                                                                    cAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGAcggcgg                                >  4:178999/1‑149 (MQ=255)
                                                                                                                          aCAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgc                                                            <  6:73544/115‑1 (MQ=255)
                                                                                                                          aCAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGtgc                                                            >  5:73544/1‑115 (MQ=255)
                                                                                                                                     cGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaa                >  6:4523/1‑148 (MQ=255)
                                                                                                                                     cGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGaa                >  5:45636/1‑148 (MQ=255)
                                                                                                                                       cggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAAc             >  3:69391/3‑149 (MQ=255)
                                                                                                                                        ggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCa                                          >  3:200133/2‑119 (MQ=255)
                                                                                                                                        ggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCa                                          <  4:200133/118‑1 (MQ=255)
                                                                                                                                        ggggCCGCGCGGGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAAcc            >  6:292824/2‑149 (MQ=255)
                                                                                                                                            ccgcccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGcc        >  7:182688/1‑149 (MQ=255)
                                                                                                                                            ccgcccgcGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGc         >  4:54346/1‑148 (MQ=255)
                                                                                                                                                  gcgcgGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCtt  >  1:60648/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTT  >  NC_002947/1676016‑1676310

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: