Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,854 C→G 6.2% G70A (GGC→GCC)  PP_0871 ← glycine betaine/carnitine/choline ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,010,8540CG 6.2% 206.9 / 3.6 65G70A (GGC→GCC) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (36/25);  new base G (2/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATC  >  NC_002947/1010711‑1011001
                                                                                                                                               |                                                                                                                                                   
gCCCGGCGCCGATGCCCAGGTTACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGCATAGGGTGGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCTCg                                                                                                                                                >  5:258504/1‑149 (MQ=255)
   cGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGc                                                                                                                                              >  7:34767/1‑148 (MQ=255)
     gcgcCGATGCCCGGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAg                                                                                                                                           >  2:213508/1‑149 (MQ=255)
          gATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGcgccg                                                                                                                                      >  4:353781/1‑149 (MQ=255)
                  ggATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATc                                                                                                                              >  1:5801/1‑149 (MQ=255)
                        tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGGCGGCGGCGCTGGCCGTGGGCATCCCCGCCGGCAttgc                                                                                                                            <  1:126742/145‑4 (MQ=37)
                        tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGGCGGCGGCGCTGGCCGTGGGCATCCCCGCCGGCAttgc                                                                                                                            >  2:126742/1‑142 (MQ=37)
                        tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGc                                                                                                                         >  1:262107/1‑148 (MQ=255)
                          aGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACg                                                                                                                                                <  8:199290/123‑1 (MQ=255)
                          aGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACg                                                                                                                                                >  7:199290/1‑123 (MQ=255)
                           gggCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAg                                                                                                                      >  3:74283/1‑148 (MQ=255)
                            ggCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCa                                                                                                                    <  5:238596/149‑1 (MQ=255)
                             gCGTTGGCCAGAACGGCCAGAGGGGGGAGGGGGGTGCCAACGGTGCAGAACTGCCTGACGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCGCTCGGCCAGCGCCGCCGCCATCGAAGCCAGCAc                                                                                                                   <  6:258504/149‑1 (MQ=255)
                             gCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCAc                                                                                                                   <  2:262107/149‑1 (MQ=255)
                                  ggCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATa                                                                                                              >  8:295774/1‑149 (MQ=255)
                                    ccAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATAcc                                                                                                            <  7:295774/149‑1 (MQ=255)
                                    ccAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATAcc                                                                                                            <  3:353781/149‑1 (MQ=255)
                                        aaCGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGt                                                                                                        >  4:203738/1‑149 (MQ=255)
                                          cGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTgc                                                                                                      <  3:203738/149‑1 (MQ=255)
                                             ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCt                                                                                                   >  3:105010/1‑149 (MQ=255)
                                             ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGCGCTGGCCGTGGGCATccccgccgg                                                                                                                                  <  1:114381/118‑2 (MQ=37)
                                             ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGCGCTGGCCGTGGGCATccccgccgg                                                                                                                                  >  2:114381/1‑117 (MQ=37)
                                                 agGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCa                                                                                               <  2:5801/149‑1 (MQ=255)
                                                    ggggATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGt                                                                                            >  8:326897/1‑149 (MQ=255)
                                                                 ccAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAg                                                                                              >  2:353251/1‑134 (MQ=255)
                                                                 ccAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAg                                                                                              <  1:353251/134‑1 (MQ=255)
                                                                  cAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTa                                                                                           >  1:108876/1‑136 (MQ=255)
                                                                  cAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGAtcg                                                                              >  4:64839/1‑149 (MQ=255)
                                                                  cAACGTTGAAGAACTGCATGAACCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTa                                                                                           <  2:108876/136‑1 (MQ=255)
                                                                   aaCGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGAtcgt                                                                             >  5:366292/1‑149 (MQ=255)
                                                                           agaaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCATTACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATaa                                                                     <  1:249023/149‑1 (MQ=255)
                                                                                  cATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGa                                                              >  2:171597/1‑149 (MQ=255)
                                                                                        gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                        <  4:105010/149‑1 (MQ=255)
                                                                                         cgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTc                                                       <  7:92977/149‑1 (MQ=255)
                                                                                          gcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCCCCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTc                                                       >  8:69081/1‑148 (MQ=255)
                                                                                             cGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                        <  2:272498/144‑1 (MQ=255)
                                                                                             cGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                        >  1:272498/1‑144 (MQ=255)
                                                                                                 tttGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                               >  3:309533/1‑149 (MQ=255)
                                                                                                   tGTCGACCCGGTGCGGTCGCCTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCGTGCAGGTAGAACCCAAGATCGTCGCGAGACCGCGCGATCGTGTCGATGCCGCTc                                             >  6:362244/1‑149 (MQ=255)
                                                                                                   tGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                               <  1:171597/147‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                             <  6:366292/148‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                             <  3:64839/148‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                             <  4:309533/148‑1 (MQ=255)
                                                                                                     tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                  <  6:340454/142‑1 (MQ=255)
                                                                                                     tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                  >  5:340454/1‑142 (MQ=255)
                                                                                                       gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                  >  2:87709/1‑140 (MQ=255)
                                                                                                       gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                  <  1:87709/140‑1 (MQ=255)
                                                                                                       gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAg                                          <  7:326897/148‑1 (MQ=255)
                                                                                                               cgggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCCGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACcagc                                 <  7:189030/147‑1 (MQ=255)
                                                                                                                cGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACcagc                                 <  1:213508/148‑1 (MQ=255)
                                                                                                                  gTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAgg                              <  4:74283/149‑1 (MQ=255)
                                                                                                                   tCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCAAGTGGACCAGCAggg                             >  1:242480/1‑149 (MQ=255)
                                                                                                                        ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                         >  6:7075/1‑148 (MQ=255)
                                                                                                                        ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                         >  1:153960/1‑148 (MQ=255)
                                                                                                                        ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCACCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                         >  8:219155/1‑148 (MQ=255)
                                                                                                                          aaCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                      >  4:232201/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                   >  4:148967/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                   >  2:201951/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                   >  8:287105/1‑149 (MQ=255)
                                                                                                                               aaTGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGCTGCCGATCCAGTGGACCAGCAGGGCCAGGCTGAc                   >  2:154763/1‑147 (MQ=255)
                                                                                                                               aaTGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGa                 <  5:7075/149‑1 (MQ=255)
                                                                                                                                     ggcggGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCAcgg           >  3:106752/1‑149 (MQ=255)
                                                                                                                                          ggATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGcc      >  2:94466/1‑149 (MQ=255)
                                                                                                                                              gCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATc  >  3:252502/1‑149 (MQ=255)
                                                                                                                                              gCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATc  >  3:247968/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATC  >  NC_002947/1010711‑1011001

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