Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,857 A→G 6.3% V69A (GTG→GCG)  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,8570AG 6.3% 188.1 / 2.5 67V69A (GTG→GCG) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/25);  new base G (2/2);  total (39/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 = 8.17e-01

GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAA  >  NC_002947/1010711‑1011005
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                                    >  7:199290/1‑123 (MQ=255)
                          aGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACg                                                                                                                                                    <  8:199290/123‑1 (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                                                            >  1:272498/1‑144 (MQ=255)
                                                                                             cGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                            <  2:272498/144‑1 (MQ=255)
                                                                                                 tttGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                                   >  3:309533/1‑149 (MQ=255)
                                                                                                   tGTCGACCCGGTGCGGTCGCCTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCGTGCAGGTAGAACCCAAGATCGTCGCGAGACCGCGCGATCGTGTCGATGCCGCTc                                                 >  6:362244/1‑149 (MQ=255)
                                                                                                   tGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                                   <  1:171597/147‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                 <  4:309533/148‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                 <  3:64839/148‑1 (MQ=255)
                                                                                                    gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                 <  6:366292/148‑1 (MQ=255)
                                                                                                     tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                      >  5:340454/1‑142 (MQ=255)
                                                                                                     tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                      <  6:340454/142‑1 (MQ=255)
                                                                                                       gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                      <  1:87709/140‑1 (MQ=255)
                                                                                                       gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                      >  2:87709/1‑140 (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                             >  1:153960/1‑148 (MQ=255)
                                                                                                                        ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                             >  6:7075/1‑148 (MQ=255)
                                                                                                                        ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCACCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                             >  8:219155/1‑148 (MQ=255)
                                                                                                                          aaCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                          >  4:232201/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                       >  2:201951/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                       >  4:148967/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)
                                                                                                                                                ccACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc    <  8:14343/149‑1 (MQ=255)
                                                                                                                                                ccACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATcc     >  7:14343/1‑148 (MQ=255)
                                                                                                                                                  aCGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCaa  <  3:232201/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAA  >  NC_002947/1010711‑1011005

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