Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,709,445 A→T 6.5% T22S (ACC→TCC) ‡ ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,709,4450AT 6.5% 180.7 / 2.5 61T22S (ACC→TCC) ‡ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGCCCGTAGC  >  NC_002947/1709299‑1709587
                                                                                                                                                   |                                                                                                                                              
ggCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGaaa                                                                                                                                                <  4:30207/147‑1 (MQ=255)
 gCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCa                                                                                                                                             >  8:4165/1‑149 (MQ=255)
     tgGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGTATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCAt                                                                                                                                            <  1:98776/147‑1 (MQ=255)
         cACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCaa                                                                                                                                     <  2:180767/149‑1 (MQ=255)
           ctctTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACa                                                                                                                                   <  1:27346/149‑1 (MQ=255)
                 ggggATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATg                                                                                                                              <  1:17453/148‑1 (MQ=255)
                     aTACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGAt                                                                                                                         <  2:54629/149‑1 (MQ=255)
                        ccGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGc                                                                                                                       <  8:116112/148‑1 (MQ=255)
                        ccGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCg                                                                                                                      >  2:292702/1‑149 (MQ=255)
                         cgctcgctCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAATGGTTTCGCCaag                                                                                                                                    <  1:187793/134‑2 (MQ=255)
                         cgctcgctCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAATGGTTTCGCCaag                                                                                                                                    >  2:187793/1‑133 (MQ=255)
                         cgctcgctCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGt                                                                                                                     <  3:196384/149‑1 (MQ=255)
                              gctcCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTcc                                                                                                                  <  1:292702/147‑1 (MQ=255)
                              gctcCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCt                                                                                                                 <  8:231353/148‑1 (MQ=255)
                                   tGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGa                                                                                                            >  2:130918/1‑148 (MQ=255)
                                    gCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTc                                                                                                                   <  1:106972/140‑1 (MQ=255)
                                    gCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTc                                                                                                                   >  2:106972/1‑140 (MQ=255)
                                        aGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAATGGTTTCGCCaag                                                                                                                                    <  2:138712/119‑2 (MQ=255)
                                        aGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAATGGTTTCGCCaag                                                                                                                                    >  1:138712/1‑118 (MQ=255)
                                        aGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATc                                                                                                       >  6:97758/1‑148 (MQ=255)
                                         gggCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGa                                                                                                     <  3:258261/149‑1 (MQ=255)
                                         gggCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGa                                                                                                     <  5:97758/149‑1 (MQ=255)
                                          ggCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAg                                                                                                    <  1:81705/149‑1 (MQ=255)
                                          ggCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAg                                                                                                    >  2:81705/1‑149 (MQ=255)
                                          ggCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAg                                                                                                    <  5:93276/149‑1 (MQ=255)
                                          ggCGCGCCGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTTGTTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAg                                                                                                    <  7:219139/149‑1 (MQ=255)
                                               gcTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACa                                                                                                >  4:311902/1‑148 (MQ=255)
                                                     gCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCca                                                                                         >  6:18249/1‑149 (MQ=255)
                                                               aTGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAg                                                                               <  1:91239/149‑1 (MQ=255)
                                                                   aaCGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTACTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCAGCCACAGTGCCAAGGcc                                                                            >  4:292648/1‑148 (MQ=255)
                                                                     cGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGAc                                                                         <  2:301680/149‑1 (MQ=255)
                                                                      ggAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGAcc                                                                        >  1:85456/1‑149 (MQ=255)
                                                                       gagagaGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCg                                                                       <  4:158557/149‑1 (MQ=255)
                                                                                       gcGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGTTGAACTACTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCAGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGag                                                       <  3:292648/148‑1 (MQ=255)
                                                                                       ccGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGag                                                       >  7:68200/1‑149 (MQ=255)
                                                                                        ggATTTCGATGTGGGTTTGCGTGAAGATGTCGATGTGTTGGGTGTACTGCTTGGCGAAACCTTTGGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAggt                                                                                                  <  1:310598/104‑2 (MQ=255)
                                                                                        cGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGa                                                                                                  <  1:305908/105‑1 (MQ=255)
                                                                                        cGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGa                                                                                                  >  2:305908/1‑105 (MQ=255)
                                                                                        cGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGagc                                                      <  1:130918/149‑1 (MQ=255)
                                                                                        cGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGag                                                       >  5:55301/1‑148 (MQ=255)
                                                                                            atCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCt                                                  >  7:275007/1‑149 (MQ=255)
                                                                                              cGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTc                                                                                                                                          >  8:284669/1‑59 (MQ=255)
                                                                                               gATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAg                                               <  8:275007/149‑1 (MQ=255)
                                                                                                        tttGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGAccc                                       >  1:139302/1‑149 (MQ=255)
                                                                                                           gCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTg                                     <  2:85456/148‑1 (MQ=255)
                                                                                                             gTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCg                                  >  1:135863/1‑149 (MQ=255)
                                                                                                              tGAAGATGTCCATGTGTTTGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCg                                  >  8:303510/1‑148 (MQ=255)
                                                                                                                 aGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGAcc                              >  3:321153/1‑149 (MQ=255)
                                                                                                                     gTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGc                          >  6:41336/1‑149 (MQ=255)
                                                                                                                      tCCATGTGTTTGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGc                          <  7:303510/148‑1 (MQ=255)
                                                                                                                      tCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGcc                         >  7:187416/1‑149 (MQ=255)
                                                                                                                      tCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGcc                         >  2:56931/1‑149 (MQ=255)
                                                                                                                         aTGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCgaa                      <  8:187416/149‑1 (MQ=255)
                                                                                                                          tgtgTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCg                        >  4:320637/1‑146 (MQ=255)
                                                                                                                            tgtTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCgaaga                    <  7:4165/148‑1 (MQ=255)
                                                                                                                              ttGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGcc                         >  6:248973/1‑141 (MQ=255)
                                                                                                                              ttGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGcc                         <  5:248973/141‑1 (MQ=255)
                                                                                                                              ttGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGaagaag                  <  5:41336/148‑1 (MQ=255)
                                                                                                                                  gTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACctg             <  1:56931/149‑1 (MQ=255)
                                                                                                                                   tGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACctgc            >  3:15147/1‑149 (MQ=255)
                                                                                                                                        tGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGccc       >  2:35086/1‑149 (MQ=255)
                                                                                                                                               cGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGCCCGTAGc  >  3:189976/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGCCCGTAGC  >  NC_002947/1709299‑1709587

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: