Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,424,884 T→C 5.0% H20R (CAT→CGT)  PP_5647 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,424,8840TC 5.0% 233.1 / 2.5 79H20R (CAT→CGT) PP_5647hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base C (2/2);  total (42/37)
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.77e-01

TTTCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATC  >  NC_002947/4424738‑4425030
                                                                                                                                                    |                                                                                                                                                  
tttCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGtgatga                                                                                                                                                  >  8:266598/1‑149 (MQ=255)
tttCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGtgatga                                                                                                                                                  >  8:138225/1‑149 (MQ=255)
        ggTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGt                                                                                                                                           <  7:215537/148‑1 (MQ=255)
           cTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCATCACGACGaaatccttg                                                                                                                                       <  3:127868/149‑8 (MQ=255)
                   gTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGTGCTTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTgg                                                                                                                                 >  1:323092/1‑149 (MQ=255)
                   gTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGt                                                                                                                               >  2:215749/1‑149 (MQ=255)
                    tAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTc                                                                                                                              >  2:94042/1‑149 (MQ=255)
                        ccTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGcc                                                                                                                          >  4:160312/1‑149 (MQ=255)
                              ctACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGAc                                                                                                                                                 >  8:97819/1‑120 (MQ=255)
                              ctACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGAc                                                                                                                                                 <  7:97819/120‑1 (MQ=255)
                                aCAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTc                                                                                                                    >  3:392250/1‑147 (MQ=255)
                                aCAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGACTGTACGCGAAAGACGGGTATGCACGTCTATACTACTGCCTTCCTCAAGGATTTCGTCGTGATGACGAAGGAGAAGCGATGGGTCTGCCCAGGTcc                                                                                                                   >  5:181538/1‑148 (MQ=255)
                                 cAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTc                                                                                                                    >  5:254881/1‑146 (MQ=255)
                                    ttCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAAt                                                                                                              <  1:94042/149‑1 (MQ=255)
                                         tCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTc                                                                                                         <  3:61290/149‑1 (MQ=255)
                                          ccccGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcg                                                                                                        >  4:197759/1‑149 (MQ=255)
                                           cccGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcgc                                                                                                       >  7:360890/1‑149 (MQ=255)
                                              gCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcgccgc                                                                                                    <  4:392250/149‑1 (MQ=255)
                                                  ttGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                                  <  7:138225/147‑1 (MQ=255)
                                                       gTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                           <  1:215749/149‑1 (MQ=255)
                                                       gTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                           <  1:58656/149‑1 (MQ=255)
                                                        tACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                           >  5:165073/1‑148 (MQ=255)
                                                          cacCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCatat                                                                                        <  7:266598/149‑1 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAgg                                                                                                                      <  7:235581/118‑1 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAgg                                                                                                                      >  8:235581/1‑118 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAg                                                                                                                       <  3:24884/117‑1 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAg                                                                                                                       >  4:24884/1‑117 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACAtctc                                                                                            <  2:325192/144‑1 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACAtctc                                                                                            >  1:325192/1‑144 (MQ=255)
                                                           acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCatat                                                                                        >  7:91660/1‑148 (MQ=255)
                                                             cccGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGATGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCaa                                                                                     >  5:152830/1‑149 (MQ=255)
                                                              ccGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCaa                                                                                     >  3:329737/1‑148 (MQ=255)
                                                                  tGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGCCATCTCATATCAATctcc                                                                                >  3:113124/1‑149 (MQ=255)
                                                                    aaGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcct                                                                               <  3:160312/148‑1 (MQ=255)
                                                                         cTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTcttact                                                                         >  3:393393/1‑149 (MQ=255)
                                                                              t‑‑ttCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATcc                                                                    >  3:270179/1‑149 (MQ=255)
                                                                                    gCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                                 <  8:91660/148‑1 (MQ=255)
                                                                                    gCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                                 <  4:270179/148‑1 (MQ=255)
                                                                                     ccTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                                 >  1:6104/1‑147 (MQ=255)
                                                                                      cTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAg                                                              <  6:165073/149‑1 (MQ=255)
                                                                                      cTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAg                                                              >  5:178872/1‑149 (MQ=255)
                                                                                         tCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                                  <  7:289978/110‑1 (MQ=255)
                                                                                         tCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                                  >  8:289978/1‑110 (MQ=255)
                                                                                           cgcgAAAGACGGGTTTGCACTTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATTAAGAAGGTGAAGCTCTGGGTCTGCACATGTCCAGAATTTTTCTCCTCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGtcgc                                                          <  6:181538/148‑1 (MQ=255)
                                                                                           cgcgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGcgcg                                                         >  4:392211/1‑149 (MQ=255)
                                                                                           cgcgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGcgcg                                                         <  4:329737/149‑1 (MQ=255)
                                                                                             cgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGt                                                       <  6:178872/149‑1 (MQ=255)
                                                                                                aaGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTaaa                                                    <  2:323092/149‑1 (MQ=255)
                                                                                                aaGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTaaa                                                    <  3:197759/149‑1 (MQ=255)
                                                                                                   aCGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAg                                                                                                                                             <  7:241590/57‑1 (MQ=255)
                                                                                                   aCGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAg                                                                                                                                             >  8:241590/1‑57 (MQ=255)
                                                                                                     gggTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCACAGTAGGCGCGGTAAagagag                                               >  8:110809/1‑149 (MQ=255)
                                                                                                      ggTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAgg                                              <  1:342793/149‑1 (MQ=255)
                                                                                                        tttGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGATGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGga                                            <  6:152830/149‑1 (MQ=255)
                                                                                                             aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgatagat                                       >  1:216737/1‑149 (MQ=255)
                                                                                                             aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgatagat                                       <  3:392211/149‑1 (MQ=255)
                                                                                                             aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgataga                                        >  6:270149/1‑148 (MQ=255)
                                                                                                             aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgataga                                        >  7:311941/1‑148 (MQ=255)
                                                                                                               gTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGatat                                     >  7:135429/1‑149 (MQ=255)
                                                                                                                   aTACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCAt                                 <  4:248339/149‑1 (MQ=255)
                                                                                                                       tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcctc                                                                              >  4:218416/1‑100 (MQ=255)
                                                                                                                       tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcctc                                                                              <  3:218416/100‑1 (MQ=255)
                                                                                                                       tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                             >  5:381354/1‑149 (MQ=255)
                                                                                                                             cGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGGCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                       <  5:82974/149‑1 (MQ=255)
                                                                                                                              gTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCTCTTAAAt                      <  5:270149/149‑1 (MQ=255)
                                                                                                                                  tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAagagag                                               <  7:361536/120‑1 (MQ=255)
                                                                                                                                  tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAagagag                                               >  8:361536/1‑120 (MQ=255)
                                                                                                                                  tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGa                  >  5:257460/1‑149 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATc                                                                     >  7:183106/1‑96 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATc                                                                     <  8:183106/96‑1 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                       >  5:323862/1‑142 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                       <  6:323862/142‑1 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGAcc                >  2:107026/1‑149 (MQ=255)
                                                                                                                                    aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGAcc                >  5:186363/1‑149 (MQ=255)
                                                                                                                                       gATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGt             <  2:6104/149‑1 (MQ=255)
                                                                                                                                       accttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTctta                                                                           <  4:334261/84‑1 (MQ=255)
                                                                                                                                       accttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTctta                                                                           >  3:334261/4‑87 (MQ=255)
                                                                                                                                        ccttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTg            >  3:376498/3‑149 (MQ=255)
                                                                                                                                        aTTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTg            <  6:254881/149‑1 (MQ=255)
                                                                                                                                            cgtcgtGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGaa         >  5:243821/1‑148 (MQ=255)
                                                                                                                                             gtcgtGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGaaaa       >  3:326193/1‑149 (MQ=255)
                                                                                                                                                 tgatgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                             >  2:264577/1‑123 (MQ=255)
                                                                                                                                                 tgatgaCGAAGGTGAAGAGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                             <  1:264577/123‑1 (MQ=255)
                                                                                                                                                   atgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATg     <  6:409558/145‑1 (MQ=255)
                                                                                                                                                   atgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATg     >  5:409558/1‑145 (MQ=255)
                                                                                                                                                   atgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATc  <  8:135429/148‑1 (MQ=255)
                                                                                                                                                   atgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATc  <  8:360890/148‑1 (MQ=255)
                                                                                                                                                   atgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATc  <  4:113124/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATC  >  NC_002947/4424738‑4425030

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: