Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,667,913 G→T 5.4% E305D (GAG→GAT ffh → signal recognition particle GTPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,667,9130GT 5.4% 263.4 / 3.6 83E305D (GAG→GATffhsignal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base T (2/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667767‑1668054
                                                                                                                                                  |                                                                                                                                               
ggTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCa                                                                                                                                               <  5:317154/149‑1 (MQ=255)
 gTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAg                                                                                                                                              <  1:178534/149‑1 (MQ=255)
  tgctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAgg                                                                                                                                             >  6:342085/1‑149 (MQ=255)
  tgctgaGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGTTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATCGGCGACGTGCTCAGCCTGATCGAGCatg                                                                                                                                             >  7:39275/1‑147 (MQ=255)
    ctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGct                                                                                                                                           <  6:219232/149‑2 (MQ=255)
       cgcTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGa                                                                                                                                         <  3:76880/148‑1 (MQ=255)
            tCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGATAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGTCTGATCGAGCAGGCCGAGCAGa                                                                                                                                    <  7:247508/148‑1 (MQ=255)
             cGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGAGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGCacgtc                                                                                                                                  <  7:255047/149‑5 (MQ=255)
                     gCAATCACCGGCAAACCGATCAAGTTCAACGGTATGGGGGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGACGCTTCGAGCATCCTCGTCATGGGCCAAGTGCTCAGCCAGATACAGCAGGCCGAGCAGACCAACGAcaa                                                                                                                          >  1:326021/1‑149 (MQ=255)
                       aaTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaagg                                                                                                                        <  7:265816/149‑1 (MQ=255)
                         tCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                       >  7:265544/1‑148 (MQ=255)
                         tCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                       >  8:20796/1‑148 (MQ=255)
                           aCCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggcc                                                                                                                      >  4:381744/1‑147 (MQ=255)
                             cGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCCCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaagg                                                                                                                  <  8:138342/149‑1 (MQ=255)
                             cGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGTGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaagg                                                                                                                  <  8:3752/149‑1 (MQ=255)
                             cGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaagg                                                                                                                  >  3:131046/1‑149 (MQ=255)
                             acccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACTTGCTCAGCCTGATCGAGCAGCCCGAGCAGACCATCGAcaaggccaagg                                                                                                                  >  8:344285/7‑149 (MQ=255)
                               ccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACaaggccaaggcc                                                                                                                >  3:323349/5‑149 (MQ=255)
                               ccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACaaggccaaggcc                                                                                                                >  5:72132/5‑149 (MQ=255)
                                cAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCGCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCg                                                                                                              <  4:249035/149‑1 (MQ=255)
                                       aTCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCt                                                                                                      >  8:225942/1‑149 (MQ=255)
                                          aaGTTCATCCCAATCGCTGAGAAGACCGAGGCCCTCGACCCGATCCACCCCGACCGCCTCGCTTACCGCATCCTCGGCATGGGCGACGTGCTCAGCCAGATCGAGCAGGCCGAGCAGACCATCGACAAGCCCAAGGACGACTAGCTCGc                                                                                                   >  2:142844/1‑149 (MQ=255)
                                            gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                                  >  7:299293/1‑148 (MQ=255)
                                            gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                                  >  8:154746/1‑148 (MQ=255)
                                            gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                                  >  3:270610/1‑148 (MQ=255)
                                              tCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                                  >  3:118436/1‑146 (MQ=255)
                                                    gTATGGGTGAGAAGCCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACt                                                                                          >  6:86708/1‑148 (MQ=255)
                                                         ggTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaagaa                                                                                    <  7:225942/149‑1 (MQ=255)
                                                         ggTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaagaa                                                                                    <  5:342085/149‑1 (MQ=255)
                                                          gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGaagaag                                                                                   <  4:267893/149‑1 (MQ=255)
                                                             agaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaaggg                                                                                 >  7:255580/1‑148 (MQ=255)
                                                               aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaa                                                                              <  7:20796/149‑1 (MQ=255)
                                                               aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                                >  5:274816/1‑147 (MQ=255)
                                                               aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                                >  5:322520/1‑147 (MQ=255)
                                                                 gaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaagg                                                                            <  2:212975/149‑1 (MQ=255)
                                                                 gaCAGAGGCCCTCGATCCGTTCCAGCCCGACCGCGTCGCTTCGCTCATGCTGGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCAAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaagg                                                                            <  8:39275/149‑1 (MQ=255)
                                                                    cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCt                                                                         >  7:154937/1‑149 (MQ=255)
                                                                    cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCt                                                                         >  6:269126/1‑149 (MQ=255)
                                                                      aGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                                       <  6:248785/149‑1 (MQ=255)
                                                                        gCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGa                                                                     <  4:131046/149‑1 (MQ=255)
                                                                         cccTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGAc                                                                    <  7:154746/149‑1 (MQ=255)
                                                                            tCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTc                                                                 >  8:35336/1‑149 (MQ=255)
                                                                               aGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCTAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                              <  8:265544/149‑1 (MQ=255)
                                                                                 ccGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGa                                                            <  4:156351/149‑1 (MQ=255)
                                                                                  cGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGAc                                                           <  6:274816/149‑1 (MQ=255)
                                                                                    ttCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACACGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGATGACtt                                                         >  5:278857/1‑149 (MQ=255)
                                                                                      tcACCCCGACCGCCTCGCTTCGCGGATCCTGGGGATTGGCGACGGGCTCATGCTGATCGATCAGGCCGAGCAGACCATCGAAAATGGCAAGGCCAGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                  <  1:223536/122‑1 (MQ=255)
                                                                                      ccACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTcc                                                       >  4:297679/1‑149 (MQ=255)
                                                                                      ccACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTcc                                                       >  7:202199/1‑149 (MQ=255)
                                                                                      ccACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTc                                                        <  5:86708/148‑1 (MQ=255)
                                                                                      ccACCCCGACCGCGTCGCTTCGCGAATCCTCGGCATGGGCGACGTGCTCAGCCTGATCTAGCAGGCCGAGCAGACCATCGCCAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                  >  2:223536/1‑122 (MQ=255)
                                                                                         ccccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCcgcg                                                    >  7:158267/1‑149 (MQ=255)
                                                                                           ccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  <  8:299293/149‑1 (MQ=255)
                                                                                             gACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  >  4:74843/1‑147 (MQ=255)
                                                                                             gACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCa                                                <  6:322520/149‑1 (MQ=255)
                                                                                                           gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATg                                   <  4:323349/148‑1 (MQ=255)
                                                                                                           gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGCCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATg                                   <  6:72132/148‑1 (MQ=255)
                                                                                                            cgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaa                                 <  3:297679/149‑1 (MQ=255)
                                                                                                            cgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaa                                 <  5:269126/149‑1 (MQ=255)
                                                                                                               aTCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaagaa                              <  5:236849/149‑1 (MQ=255)
                                                                                                                   tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATg                          >  5:83136/1‑149 (MQ=255)
                                                                                                                   tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATg                          >  7:78541/1‑149 (MQ=255)
                                                                                                                        aTGTGCGACTTTGTCAGACTTATCGAGAATGAGGAGCAGACCATCGACAAGTCCAAGGACGAAAAGCTTGCCAAGAAACTTAATAATGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGggcg                      <  7:344285/148‑1 (MQ=255)
                                                                                                                         tGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGggcggc                    >  5:397438/1‑149 (MQ=255)
                                                                                                                            gCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTc                 >  5:83161/1‑149 (MQ=255)
                                                                                                                                  tGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGCCTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           >  2:252274/1‑149 (MQ=255)
                                                                                                                                  tGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGATCATGGGCGTCCTCggcggc           <  8:202199/149‑1 (MQ=255)
                                                                                                                                  tGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           <  6:83161/149‑1 (MQ=255)
                                                                                                                                    cTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt         >  3:231273/1‑149 (MQ=255)
                                                                                                                                    cTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt         <  5:289878/149‑1 (MQ=255)
                                                                                                                                    cTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt         <  8:158267/149‑1 (MQ=255)
                                                                                                                                      cggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGa       >  3:337009/3‑149 (MQ=255)
                                                                                                                                      cAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTgca                                          <  5:234501/114‑1 (MQ=255)
                                                                                                                                      cAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTgca                                          >  6:234501/1‑114 (MQ=255)
                                                                                                                                       aGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGAt      <  8:154937/149‑1 (MQ=255)
                                                                                                                                        gCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGCGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     <  4:270610/149‑1 (MQ=255)
                                                                                                                                         ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggc              >  7:240426/1‑139 (MQ=255)
                                                                                                                                         ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggc              <  8:240426/139‑1 (MQ=255)
                                                                                                                                         ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     >  7:113793/1‑148 (MQ=255)
                                                                                                                                            gATCGAGCAGGCCGAGCAGCCCATCGACAAGGCCAAGGCCGAGCAAGCTTGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGAACTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCCGCCTGATGGa   >  1:273743/1‑148 (MQ=255)
                                                                                                                                             aTCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc  >  8:132819/1‑147 (MQ=255)
                                                                                                                                             aTCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc  >  7:361909/1‑147 (MQ=255)
                                                                                                                                             aTCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc  >  5:154306/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667767‑1668054

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: