Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,831,463 G→C 6.2% R221G (CGC→GGC)  PP_3384 ← gluconate 2‑dehydrogenase gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,831,4630GC 6.2% 209.6 / 4.3 67R221G (CGC→GGC) PP_3384gluconate 2‑dehydrogenase gamma subunit
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (2/2);  total (35/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 = 5.68e-01

GATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCGAAAACCAGCTCTCC  >  NC_002947/3831317‑3831607
                                                                                                                                                  |                                                                                                                                                
gATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCtgg                                                                                                                                                >  4:3325/1‑149 (MQ=255)
 aTCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCCGTCCATGAAATCGGCGCGGgc                                                                                                                                               >  2:314991/1‑149 (MQ=255)
    gCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGcc                                                                                                                                            >  2:354207/1‑149 (MQ=255)
      aCCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCgg                                                                                                                                          >  1:207785/1‑149 (MQ=255)
          gCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGaa                                                                                                                                      <  1:147604/149‑1 (MQ=255)
           cccAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAg                                                                                                                                     <  8:350419/149‑1 (MQ=255)
            ccAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGc                                                                                                                                    >  5:202771/1‑149 (MQ=255)
            ccAGCAGAAGACGACGATCACCGAATCCACTTTGTTAAACACAGTCGCCATGACGCAGGCCCTCTCACCACGTATCGATACAGACGGGAACGTGTAAGGATCATTGCGATCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGc                                                                                                                                    >  8:241905/1‑149 (MQ=255)
                 cgaagccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATc                                                                                                                               <  1:12014/149‑1 (MQ=255)
                  gaagccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCTCTCGACCAAGTCCATGAAATCGGCGCGGGAGCCGGGGAAGCCGATCt                                                                                                                              <  8:162592/149‑1 (MQ=255)
                    agccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGa                                                                                                                            <  7:8284/149‑1 (MQ=255)
                    agccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGa                                                                                                                            <  1:314991/149‑1 (MQ=255)
                     gccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAg                                                                                                                           <  8:23007/149‑1 (MQ=255)
                         acgaTCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCa                                                                                                                       <  4:41262/149‑1 (MQ=255)
                                   cATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCAtt                                                                                                             >  7:89140/1‑149 (MQ=255)
                                       cACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCtt                                                                                                         >  3:309557/1‑149 (MQ=255)
                                       cACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCtt                                                                                                         <  1:354207/149‑1 (MQ=255)
                                        tatttGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACTGGTAAGGCTCATTGAGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                                                        >  7:242900/3‑149 (MQ=255)
                                        aCTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                                                        <  2:18210/149‑1 (MQ=255)
                                        aCTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                                                        <  6:202771/149‑1 (MQ=255)
                                        aCTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                                                        >  8:102845/1‑149 (MQ=255)
                                        aCTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                                                        >  1:189256/1‑149 (MQ=255)
                                          tttGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGGGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttgt                                                                                                      >  5:317719/1‑149 (MQ=255)
                                              ttCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCTGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgc                                                                                                   <  3:3325/148‑1 (MQ=255)
                                              ttCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgct                                                                                                  <  6:317719/149‑2 (MQ=255)
                                              ttCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgcc                                                                                                  >  2:180067/1‑149 (MQ=255)
                                                 aaCACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCTATACAGCCTGGAACGGGTAAGGCTCCTTGCGCTCGACCCAGTCCATGAAATAGGCGCGGGCGCAGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgtcgcc                                                                                               <  7:241905/149‑1 (MQ=255)
                                                 aaCACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgccgc                                                                                                >  8:333173/1‑148 (MQ=255)
                                                 aaCACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgccgc                                                                                                >  6:31229/1‑148 (MQ=255)
                                                  acacCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCAGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCa                                                                                              >  8:105797/1‑149 (MQ=255)
                                                   cacCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCATCCGACCATTCCTTTGTTGCCGCCAt                                                                                             <  4:309557/149‑1 (MQ=255)
                                                            cATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGt                                                                                    <  2:189256/149‑1 (MQ=255)
                                                            cATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCCCACCACCCTGTCAATGAATTCGGCGCGGGCGCCCGGGAAGCCGATCTGAGACCACCCGACCATACCTTTGTTGACGCCTTGAATCGGGt                                                                                    >  2:259122/1‑149 (MQ=255)
                                                             aTGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTc                                                                                   >  2:323392/1‑149 (MQ=255)
                                                                  tCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCag                                                                               >  7:247882/1‑148 (MQ=255)
                                                                        ccTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGcc                                                                                                                    <  8:332863/105‑1 (MQ=255)
                                                                        ccTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGcc                                                                                                                    >  7:332863/1‑105 (MQ=255)
                                                                              accacGGATCGATACACCCGGGATCGGGTATGGGTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCa                                                                   <  8:247882/148‑1 (MQ=255)
                                                                                        gATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgccgcc                                                                                               <  8:151051/110‑1 (MQ=255)
                                                                                        gATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTgccgcc                                                                                               >  7:151051/1‑110 (MQ=255)
                                                                                        gATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGtttt                                                         <  2:207785/148‑1 (MQ=255)
                                                                                         aTACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTg                                                        >  2:77878/1‑148 (MQ=255)
                                                                                             aGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGttt                                                          >  8:302286/1‑142 (MQ=255)
                                                                                             aGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGttt                                                          <  7:302286/142‑1 (MQ=255)
                                                                                                       gggTAAGGCTCGTTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTACCTCGTGTTTTGCACCAGCAGGCTg                                           >  8:331932/1‑147 (MQ=255)
                                                                                                       gggTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTgaa                                         <  8:89140/149‑1 (MQ=255)
                                                                                                            aGGCTCATTGCGCTGCACCGAGTCGATGAAATCGCCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGATTCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACGAGCAGGCTGAAGAaca                                    <  1:259122/149‑1 (MQ=255)
                                                                                                             ggCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAacac                                   >  8:86759/1‑149 (MQ=255)
                                                                                                               cTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCt                                 <  1:77878/149‑1 (MQ=255)
                                                                                                               cTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCt                                 >  4:207645/1‑149 (MQ=255)
                                                                                                                tCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTgaa                                         >  3:100749/1‑140 (MQ=255)
                                                                                                                tCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTgaa                                         <  4:100749/140‑1 (MQ=255)
                                                                                                                 cATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTg                               <  7:333173/149‑1 (MQ=255)
                                                                                                                  gttGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTg                               <  7:331932/147‑1 (MQ=255)
                                                                                                                         tCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGACCACCTTGGCCGGAAc                       >  5:309761/1‑149 (MQ=255)
                                                                                                                         tCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAAc                       >  2:223567/1‑149 (MQ=255)
                                                                                                                            aCCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGa                    <  1:323392/149‑1 (MQ=255)
                                                                                                                               cAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTc                 <  1:223567/149‑1 (MQ=255)
                                                                                                                                aGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTc                 >  7:81763/1‑148 (MQ=255)
                                                                                                                                aGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCg                <  5:31229/149‑1 (MQ=255)
                                                                                                                                      tGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCGAAAAc           >  3:158731/1‑148 (MQ=255)
                                                                                                                                       gAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCGAAAAc           >  3:188021/1‑147 (MQ=255)
                                                                                                                                         ccccGGCGCCCGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGAtt                  <  2:264511/135‑1 (MQ=255)
                                                                                                                                         ccccGGCGCCCGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGAtt                  >  1:264511/4‑138 (MQ=255)
                                                                                                                                         aaTCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCGAAAACCAGc       >  1:356795/1‑149 (MQ=255)
                                                                                                                                            cGGCGCCCGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTg                               >  7:93726/1‑122 (MQ=255)
                                                                                                                                            cGGCGCCCGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTg                               <  8:93726/122‑1 (MQ=255)
                                                                                                                                              gcgcGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACCCCTTGGCCGGAACGGATTCGAAAACCAGCTCTcc  >  6:120618/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGCCGGAACGGATTCGAAAACCAGCTCTCC  >  NC_002947/3831317‑3831607

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: ‑