Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 887,948 C→A 8.8% R328R (CGG→AGG)  PP_0768 → TPR domain‑containing response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947887,9480CA 8.8% 183.2 / 2.6 69R328R (CGG→AGG) PP_0768TPR domain‑containing response regulator
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base A (3/4);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGA  >  NC_002947/887803‑888095
                                                                                                                                                    |                                                                                                                                                   
cGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCggg                                                                                                                                                  >  6:281680/1‑149 (MQ=255)
 gAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGt                                                                                                                                                 >  7:86506/1‑149 (MQ=255)
  aaGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTCGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTc                                                                                                                                                >  1:85636/1‑149 (MQ=255)
   aGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCg                                                                                                                                               <  3:53599/149‑1 (MQ=255)
    gACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCga                                                                                                                                              <  8:86506/149‑1 (MQ=255)
       ttCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTc                                                                                                                                                >  4:165055/1‑144 (MQ=255)
       ttCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCgagac                                                                                                                                           >  3:126659/1‑148 (MQ=255)
       ttCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCgaga                                                                                                                                            >  3:43902/1‑148 (MQ=255)
         cGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTTGAGAt                                                                                                                                           <  4:258885/147‑1 (MQ=255)
              gcgcATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAacac                                                                                                                                     >  8:100870/1‑148 (MQ=255)
                gcATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAGAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGAAGGCTTCGGCCTGGACGCCCGTACCAGGGTCGAGATCAAcacca                                                                                                                                   >  1:101442/1‑149 (MQ=255)
                 cATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACg                                                                                                                                 >  3:224126/1‑149 (MQ=255)
                           gTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAg                                                                                                                       >  4:37410/1‑149 (MQ=255)
                              ccGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGTTTTTTTTTTTATGATAGTTCGAGGTTCGAGATCTACACCACGCTCAGCGAAGTg                                                                                                                     <  4:64478/148‑1 (MQ=255)
                                   aTGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGTGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGcc                                                                                                                  >  1:116856/1‑147 (MQ=255)
                                    tGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAgg                                                                                                              <  3:314211/149‑1 (MQ=255)
                                        gTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGaaaa                                                                                                          >  7:388875/1‑149 (MQ=255)
                                           aaCCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAAccc                                                                                                       >  6:37824/1‑149 (MQ=255)
                                           aaCCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAAc                                                                                                         <  5:281680/147‑1 (MQ=255)
                                              cAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAAccct                                                                                                      >  8:115691/1‑146 (MQ=255)
                                               aGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTATTTATTTATGAACGTACGAGGTTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAg                                                                                                   <  2:101442/149‑1 (MQ=255)
                                                   gtagagTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCTGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGa                                                                                                    <  6:50483/142‑1 (MQ=255)
                                                   gtagagTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGa                                                                                                    >  5:50483/3‑144 (MQ=255)
                                                   gCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAAccc                                                                                                       <  1:217204/141‑1 (MQ=255)
                                                   gCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAAccc                                                                                                       >  2:217204/1‑141 (MQ=255)
                                                      gagTTCGCGCTAGAAGGATGCCGAAAACAAGCTTGGCCTGGTGCAGGCGTTGAAGAGCAAGATTGCAGGTTTTGGGCTGGATGCACGTTCGCGGTTCGAGATCTAGACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAg                                                                                            <  2:85636/149‑1 (MQ=255)
                                                             cgcTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTg                                                                                       >  3:118810/1‑147 (MQ=255)
                                                             cgcTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCa                                                                                     <  5:37824/149‑1 (MQ=255)
                                                                tACAAGGATGCCGAAAACAACCTAGGCCTGGTGCAGGAGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGCAAAACCCCTAGGACCAAGGCCTGCAGgt                                                                                  >  6:271552/1‑149 (MQ=255)
                                                                  cAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAAACAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCAGGGTCGATAACAAAACAACTATCAGCTAAGTGGCAAAGTaaaa                                                                                                          >  1:84833/1‑124 (MQ=255)
                                                                     aGGATGCCTAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGTCCTTGACGCCCGTACCCGGGTAGAGATCAACACCTCGCTCAGCGAAGTGGCCATGGAAAACCCCGAGGACCAAGGCCTGCAGgtgcgt                                                                              <  5:271552/149‑1 (MQ=255)
                                                                              aaaaCAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGgatg                                                                     <  3:37410/149‑1 (MQ=255)
                                                                                aacaacCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGTCTTTTTATTTTATGATAGTTCGAGGTTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGgatgat                                                                   <  4:43902/149‑1 (MQ=255)
                                                                                aacaacCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGgatgat                                                                   >  7:245577/1‑149 (MQ=255)
                                                                                   aacCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGatgatg                                                                  >  3:90501/1‑147 (MQ=255)
                                                                                     ccTTGGCCTGGTGCAGGTGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAgg                                                              >  8:13810/1‑149 (MQ=255)
                                                                                         gtccTGGTTCAGGCGTTGATGAGCAAGAATGCAGGCTTCGTCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCt                                                            <  8:388875/145‑1 (MQ=255)
                                                                                            cTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAg                                                       >  1:88651/1‑149 (MQ=255)
                                                                                                gggCAGGCGTTGATGAGCAAGAATGCAGGCTTTGTTCTTTATGATCGTTCGAGGTGCGAGATCTACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGc                                                      <  8:245577/145‑1 (MQ=255)
                                                                                                gTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCTTCCTTGACGCTAGTACGCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGTAAAACCCCGAGGACCAAGGCCTGCa                                                                                     <  4:282274/116‑1 (MQ=255)
                                                                                                gTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCa                                                                                     >  3:282274/1‑116 (MQ=255)
                                                                                                gTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGc                                                      <  4:118810/147‑1 (MQ=255)
                                                                                                  gCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTg                                                   >  2:74531/1‑148 (MQ=255)
                                                                                                   cAGGCGTTGATGAGCAAGAATGCAGGCTTCGTCTTGTACGCTCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTgc                                                  <  2:88651/148‑1 (MQ=255)
                                                                                                   cAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTgc                                                  <  7:100870/148‑1 (MQ=255)
                                                                                                     ggCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTgcag                                                >  1:28952/1‑148 (MQ=255)
                                                                                                        gTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAgg                                            >  1:92793/1‑149 (MQ=255)
                                                                                                           gatgaGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGcc                                          <  3:165055/148‑1 (MQ=255)
                                                                                                           gatgaGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCg                                         >  3:281585/1‑149 (MQ=255)
                                                                                                            atgaGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCgg                                        <  1:74531/149‑1 (MQ=255)
                                                                                                                  cAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGAcc                                   >  7:346394/1‑149 (MQ=255)
                                                                                                                   aaGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGAccc                                  <  8:282776/149‑1 (MQ=255)
                                                                                                                   aaGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGAccc                                  >  7:282776/1‑149 (MQ=255)
                                                                                                                     gAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTgc                                                  <  4:15141/131‑1 (MQ=255)
                                                                                                                     gAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTgc                                                  >  3:15141/1‑131 (MQ=255)
                                                                                                                         gCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGgatg                                                                     <  3:49737/108‑1 (MQ=255)
                                                                                                                         gCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGgatg                                                                     >  4:49737/1‑108 (MQ=255)
                                                                                                                           aGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAAcgg                          >  5:137389/1‑149 (MQ=255)
                                                                                                                             gCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAAcggcc                        >  6:75332/1‑149 (MQ=255)
                                                                                                                             cctacGGACTGGACGCCCGAACCAGGTTCGAGATCAACACCACGCTCAGCTAAGTGGCCAAGGAAAACCCCGAGGAACAAGGCCTGCACTTGCGTGCGCGCATGAAGAAGGCTGCCACCCTGCAGCAGGCCGGCGACCCTGAAAcgacc                        >  6:413984/5‑149 (MQ=255)
                                                                                                                               ttCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAAcggccgg                      >  7:278258/1‑149 (MQ=255)
                                                                                                                                   gCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAg                  >  4:53078/1‑149 (MQ=255)
                                                                                                                                    ccTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAg                                                                                    >  3:91478/1‑82 (MQ=255)
                                                                                                                                    ccTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAg                                                                                    <  4:91478/82‑1 (MQ=255)
                                                                                                                                    ccTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACGCGGAAACGGCCGGCAAGt                 <  6:137389/149‑1 (MQ=255)
                                                                                                                                         aCGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACa            <  8:278258/149‑1 (MQ=255)
                                                                                                                                         aCGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACa            <  4:90501/149‑1 (MQ=255)
                                                                                                                                           gCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGa          <  5:75332/149‑1 (MQ=255)
                                                                                                                                               gttcgcGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACTGAAc        <  2:195507/142‑1 (MQ=255)
                                                                                                                                                tACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAgg     <  7:115691/149‑1 (MQ=255)
                                                                                                                                                   gcGGTTCGAGATCTACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGa  <  2:28952/148‑1 (MQ=255)
                                                                                                                                                   ccGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCACGTGCGTGCGCGGATGATGAAGGCTGCCCGCCTGCAGCAGGACGGCGCCCCGGAAACGGCCGGCAAGTTTACAGAACAGGCGa  >  2:309622/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAGCAAGAATGCAGGCTTCGGCCTGGACGCCCGTACCCGGGTCGAGATCAACACCACGCTCAGCGAAGTGGCCAAGGAAAACCCCGAGGACCAAGGCCTGCAGGTGCGTGCGCGGATGATGAAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGA  >  NC_002947/887803‑888095

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: