Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,676,160 C→G 8.5% intergenic (+22/‑38) 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,1600CG 8.5% 184.3 / 7.1 59intergenic (+22/‑38)hom/thrChomoserine dehydrogenase/threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base G (3/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GTGCCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCC  >  NC_002947/1676013‑1676304
                                                                                                                                                   |                                                                                                                                                
gTGCCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGccgccc                                                                                                                                                 >  4:265922/1‑149 (MQ=255)
   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                                                                                                                           >  4:269140/1‑149 (MQ=255)
                      aTGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATc                                                                                                                           <  3:9244/149‑1 (MQ=255)
                         gCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                        <  5:136710/149‑1 (MQ=255)
                         gCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGaa                                                                                                                        <  5:739/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                                                                                                       >  3:303169/1‑149 (MQ=255)
                                          aTCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATgcgc                                                                                                       >  5:15376/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                                                                                                >  3:201343/1‑134 (MQ=255)
                                                                tGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCa                                                                                                <  4:201343/134‑1 (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                                                                    <  2:29867/127‑1 (MQ=255)
                                                                                                   cgcATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGaa                                                                    >  1:29867/1‑127 (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                                                               <  2:16029/113‑1 (MQ=255)
                                                                                                                      tCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCg                                                               >  1:16029/1‑113 (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)
                                                                                                                                                   |                                                                                                                                                
GTGCCAATGATCCTGGTGACCCATGGCGTGGTCGAGCAGCGTATCAATGACGCCATCGTCGCTCTGGAAGCCCTGCAGGACGTGGTCGGCAAGGTTGTGCGCATTCGCGTCGAGCAGCTCAACTAACAAATATTGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCC  >  NC_002947/1676013‑1676304

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: