Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,851 A→C 6.4% I71S (ATC→AGC)  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,8510AC 6.4% 193.7 / 4.3 63I71S (ATC→AGC) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base C (2/2);  total (36/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 = 6.56e-01

GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCC  >  NC_002947/1010711‑1010997
                                                                                                                                            |                                                                                                                                                  
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                                              >  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                     >  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               >  8:287105/1‑149 (MQ=255)
                                                                                                                             gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc               >  4:148967/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)
                                                                                                                                            |                                                                                                                                                  
GCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCC  >  NC_002947/1010711‑1010997

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