Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,020 C→G 8.9% K68N (AAG→AAC)  grpE ← heat shock protein GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,0200CG 8.9% 241.0 / 11.7 79K68N (AAG→AAC) grpEheat shock protein GrpE
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base G (5/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

AAGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGC  >  NC_002947/5377873‑5378168
                                                                                                                                                   |                                                                                                                                                    
aaGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACt                                                                                                                                                     >  5:139887/1‑149 (MQ=255)
aaGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACt                                                                                                                                                     >  3:247308/1‑149 (MQ=255)
  ggTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTtc                                                                                                                                                   >  4:372828/1‑149 (MQ=255)
   gTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTtct                                                                                                                                                  <  3:133103/149‑1 (MQ=255)
      gTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTcc                                                                                                                                                <  3:292441/148‑1 (MQ=255)
       tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCTGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCa                                                                                                                                               <  1:95344/148‑1 (MQ=255)
       tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCa                                                                                                                                               <  5:87661/148‑1 (MQ=255)
         gcgcTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAGTTCTCCAgc                                                                                                                                             <  5:297154/148‑2 (MQ=255)
         gcgcTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAggg                                                                                                                                            <  5:109388/149‑1 (MQ=255)
         gcgcTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAgg                                                                                                                                             <  3:89201/148‑1 (MQ=255)
               gAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAAc                                                                                                                                       <  3:361826/148‑1 (MQ=255)
                     gtCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAggg                                                                                                                                            >  8:582/1‑137 (MQ=255)
                     gtCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAggg                                                                                                                                            <  7:582/137‑1 (MQ=255)
                      tCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGc                                                                                                                               <  6:139887/149‑1 (MQ=255)
                           gAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtct                                                                                                                          <  2:256380/149‑1 (MQ=255)
                           gAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtct                                                                                                                          >  5:48351/1‑149 (MQ=255)
                             aCATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTcc                                                                                                                        <  2:223556/149‑1 (MQ=255)
                              cATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCGAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCca                                                                                                                       <  7:241506/149‑1 (MQ=255)
                                tcctCAGGGTCATCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGTCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCTAACTTGTGGGCTTTCTCcaca                                                                                                                     <  8:151591/146‑1 (MQ=255)
                                  ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATc                                                                                                                   >  1:43186/1‑149 (MQ=255)
                                  ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATc                                                                                                                   >  5:235272/1‑149 (MQ=255)
                                  ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATc                                                                                                                   >  4:26427/1‑149 (MQ=255)
                                    cAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGTCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATc                                                                                                                   <  8:151711/147‑1 (MQ=255)
                                       ggTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGGCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCt                                                                                                              >  1:282041/1‑149 (MQ=255)
                                        gTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTc                                                                                                             <  6:235272/149‑1 (MQ=255)
                                              tCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTCGGCTTTCTCCACAACCTGCTCGGCac                                                                                                        >  5:355653/1‑148 (MQ=255)
                                                         gcgGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtacg                                                                                             <  6:88019/148‑1 (MQ=255)
                                                          cgGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGtt                                                                                           >  1:200715/1‑149 (MQ=255)
                                                           ggATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGtt                                                                                           <  8:167627/148‑1 (MQ=255)
                                                           ggATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGtt                                                                                           >  7:167627/1‑148 (MQ=255)
                                                           ggATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTc                                                                                          <  4:195090/149‑1 (MQ=255)
                                                             aTCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                        >  2:46173/1‑149 (MQ=255)
                                                             aTCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                        >  4:48047/1‑149 (MQ=255)
                                                              tCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGAGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                        >  8:54954/1‑148 (MQ=255)
                                                                    tGACGTGCTCATCGTCTGCACTGGAGTTAGCCAGTGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTc                                                                                  <  6:355653/148‑1 (MQ=255)
                                                                    tGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCg                                                                                 >  7:296697/1‑149 (MQ=255)
                                                                     gACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCgg                                                                                >  3:292232/1‑149 (MQ=255)
                                                                     gACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCg                                                                                 <  4:247308/148‑1 (MQ=255)
                                                                      aCGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAGTTCTCCAGAGCGAACATGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCgg                                                                                >  1:129920/1‑148 (MQ=255)
                                                                        gTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTa                                                                             <  3:372828/149‑1 (MQ=255)
                                                                              aTCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGgc                                                                         <  6:48351/147‑1 (MQ=255)
                                                                              aTCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGATCTGCAGGTCGGCTACGgcgc                                                                       >  3:196981/1‑149 (MQ=255)
                                                                               tcgtcgGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGAGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGgcgcg                                                                      >  4:79682/1‑149 (MQ=255)
                                                                                  tcgGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAg                                                                   >  3:365841/1‑149 (MQ=255)
                                                                                    ggCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAgc                                                                  <  7:131920/148‑1 (MQ=255)
                                                                                     gCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTATGGCGCGCAgcgc                                                                >  4:313160/1‑149 (MQ=255)
                                                                                      cACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGcc                                                               >  4:118610/1‑149 (MQ=255)
                                                                                              tGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTATGGCGCGCAGCGCCTGGTCCtt                                                       <  3:313160/149‑1 (MQ=255)
                                                                                                aGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTgg                                                     >  1:259250/1‑149 (MQ=255)
                                                                                                      ggCGAGTTGCAGGCTGTCGATCTCCGGCAGCAAGTCGGTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCAGATCCTGCTCGGGACGACGACGTACGTTCTGCAGGTCGGGTACGGCGCGCAGCGCCTGGTGGTTGGCTGCgg                                               <  4:196981/149‑1 (MQ=255)
                                                                                                       gCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGc                                              >  4:55891/1‑149 (MQ=255)
                                                                                                       gCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGc                                              <  2:259250/149‑1 (MQ=255)
                                                                                                           gTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGc                                          <  2:43186/149‑1 (MQ=255)
                                                                                                                   cTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTc                                  <  3:118610/149‑1 (MQ=255)
                                                                                                                    tGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTc                                  >  5:171231/1‑148 (MQ=255)
                                                                                                                      tCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCg                                 <  5:222493/147‑1 (MQ=255)
                                                                                                                      tCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAg                               <  8:296697/149‑1 (MQ=255)
                                                                                                                        gATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCt                             >  6:154170/1‑149 (MQ=255)
                                                                                                                            aCCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTc                            >  4:261021/1‑146 (MQ=255)
                                                                                                                            aCCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAg                         <  2:129920/149‑1 (MQ=255)
                                                                                                                            aCCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAg                         <  3:55891/149‑1 (MQ=255)
                                                                                                                             ccGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAg                         >  2:92217/1‑148 (MQ=255)
                                                                                                                             ccGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGa                        <  3:26427/149‑1 (MQ=255)
                                                                                                                                gcagcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAaca                     >  6:86101/1‑149 (MQ=255)
                                                                                                                                  agcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACg                   >  7:327842/1‑149 (MQ=255)
                                                                                                                                   gcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACgg                  >  7:244479/1‑149 (MQ=255)
                                                                                                                                     aaGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGt                <  2:282041/149‑1 (MQ=255)
                                                                                                                                      aGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTa               >  6:208475/1‑149 (MQ=255)
                                                                                                                                       gTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTTCGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTc                            <  5:316026/135‑1 (MQ=255)
                                                                                                                                       gTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTc                            >  6:316026/1‑135 (MQ=255)
                                                                                                                                       gTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTAt              >  1:210357/1‑149 (MQ=255)
                                                                                                                                        gcccTGGAGAAGTTCTCCA‑‑GCG‑ACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGc          >  2:241603/4‑149 (MQ=255)
                                                                                                                                         cccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTAt              >  4:59141/3‑147 (MQ=255)
                                                                                                                                          ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTcc                                    >  7:337117/2‑124 (MQ=255)
                                                                                                                                          ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTcc                                    <  8:337117/123‑1 (MQ=255)
                                                                                                                                          ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAg           >  4:288786/2‑149 (MQ=255)
                                                                                                                                            tGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCt         >  5:62277/1‑149 (MQ=255)
                                                                                                                                              gagaACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTg        >  3:362967/1‑148 (MQ=255)
                                                                                                                                                  aCTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCAtt    <  5:86101/148‑1 (MQ=255)
                                                                                                                                                   cTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGc  >  7:303725/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGC  >  NC_002947/5377873‑5378168

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