Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,071,291 G→A 5.3% A87A (GCC→GCT)  PP_0924 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,071,2910GA 5.3% 789.9 / 20.8 245A87A (GCC→GCT) PP_0924hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (118/114);  new base A (3/10);  total (121/124)
Fisher's exact test for biased strand distribution p-value = 8.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 251 total aligned reads displayed.

GCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGAA  >  NC_002947/1071142‑1071440
                                                                                                                                                     |                                                                                                                                                     
gCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGctgctg                                                                                                                                                       >  2:871648/1‑150 (MQ=255)
gCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGctgctg                                                                                                                                                       <  1:3093117/150‑1 (MQ=255)
       cGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCt                                                                                                                                                <  2:3326672/150‑1 (MQ=255)
       cGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCt                                                                                                                                                <  2:688285/150‑1 (MQ=255)
        gCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCACCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTg                                                                                                                                               <  2:4136995/150‑1 (MQ=255)
        gCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTg                                                                                                                                               <  2:4318074/150‑1 (MQ=255)
         ccccACAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                              <  1:3257027/150‑1 (MQ=255)
         ccccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                              <  2:948395/150‑1 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGt                                                                                                                                             >  2:1507786/1‑150 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGt                                                                                                                                             >  2:4491772/1‑150 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGt                                                                                                                                             >  2:2570521/1‑150 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGt                                                                                                                                             >  2:2261681/1‑150 (MQ=255)
           ccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTg                                                                                                                                            >  2:3242309/1‑150 (MQ=255)
            cAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGa                                                                                                                                           <  1:3991032/150‑1 (MQ=255)
            cAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGa                                                                                                                                           <  1:2255841/150‑1 (MQ=255)
              aaGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGACCAGGCCGCCAGCAGCCTgggccgt                                                                                                                                         <  1:4876940/150‑6 (MQ=255)
                      cGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGccgaa                                                                                                                                 <  2:4090541/150‑5 (MQ=255)
                           gTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtg                                                                                                                            >  1:3642644/1‑150 (MQ=255)
                            tACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                           >  1:5472856/1‑150 (MQ=255)
                              ccgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         <  2:3462102/150‑1 (MQ=255)
                              ccgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         <  2:2028590/150‑1 (MQ=255)
                              ccgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTACTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         <  2:3213616/150‑1 (MQ=255)
                               cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGtggtgg                                                                                                                        >  2:40606/1‑150 (MQ=255)
                                 ccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCa                                                                                                                      >  1:2705470/1‑150 (MQ=255)
                                        ggTAGTCGGCGGCGACTCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACCGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGtt                                                                                                               >  2:1101524/1‑150 (MQ=255)
                                            gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGcc                                                                                                           <  1:1590947/150‑1 (MQ=255)
                                            gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGcc                                                                                                           <  2:3625079/150‑1 (MQ=255)
                                             tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCg                                                                                                          >  2:483448/1‑150 (MQ=255)
                                                ggggCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTg                                                                                                       <  2:800850/150‑1 (MQ=255)
                                                     gACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAg                                                                                                  <  2:171083/150‑1 (MQ=255)
                                                     gACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAg                                                                                                  >  2:2265562/1‑150 (MQ=255)
                                                           ggTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcgcg                                                                                            <  1:579125/150‑1 (MQ=255)
                                                                ggCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       >  1:3946663/1‑150 (MQ=255)
                                                                ggCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAAGTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       <  2:2973955/150‑1 (MQ=255)
                                                                 gCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACAcc                                                                                      >  2:5499449/1‑150 (MQ=255)
                                                                    cAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAAt                                                                                   >  1:2367766/1‑150 (MQ=255)
                                                                    cAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAAt                                                                                   >  1:4419473/1‑150 (MQ=255)
                                                                           aGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggcgg                                                                            <  2:1983140/150‑1 (MQ=255)
                                                                                cGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTg                                                                       >  1:740769/1‑150 (MQ=255)
                                                                                 gggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGc                                                                      <  2:1171543/150‑1 (MQ=255)
                                                                                  ggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCa                                                                     >  2:3323441/1‑150 (MQ=255)
                                                                                  ggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCa                                                                     >  1:3480566/1‑150 (MQ=255)
                                                                                         cAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                              <  2:2363609/150‑1 (MQ=255)
                                                                                         cAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                              <  2:3596390/150‑1 (MQ=255)
                                                                                          aGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAg                                                             <  1:5175000/150‑1 (MQ=255)
                                                                                               cAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTACTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcgtc                                                        <  2:4857378/150‑1 (MQ=255)
                                                                                                aaCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTcc                                                       <  1:961437/150‑1 (MQ=255)
                                                                                                 aCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCa                                                      <  2:1320875/150‑1 (MQ=255)
                                                                                                  cAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCAt                                                     >  2:3362899/1‑150 (MQ=255)
                                                                                                     gCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTc                                                  <  1:4491772/150‑1 (MQ=255)
                                                                                                     gCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTc                                                  >  2:4220591/1‑150 (MQ=255)
                                                                                                       cGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCa                                                >  1:4728150/1‑150 (MQ=255)
                                                                                                        gTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAg                                               <  1:3891779/150‑1 (MQ=255)
                                                                                                         tGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAgg                                              >  1:4386253/1‑150 (MQ=255)
                                                                                                         tGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAgg                                              <  2:1387589/150‑1 (MQ=255)
                                                                                                            aGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGccc                                           >  2:5601612/1‑150 (MQ=255)
                                                                                                             ggCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCt                                          >  2:3884255/1‑150 (MQ=255)
                                                                                                                tGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCt                                       >  1:1026234/1‑150 (MQ=255)
                                                                                                                tGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCt                                       >  1:3847628/1‑150 (MQ=255)
                                                                                                                    ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGc                                   >  1:4992579/1‑150 (MQ=255)
                                                                                                                    ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGc                                   <  1:1294018/150‑1 (MQ=255)
                                                                                                                     tcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCa                                  >  2:4122075/1‑150 (MQ=255)
                                                                                                                      cttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAc                                 >  2:1523266/1‑150 (MQ=255)
                                                                                                                       ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                                >  2:2759601/1‑150 (MQ=255)
                                                                                                                       ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                                >  1:3523478/1‑150 (MQ=255)
                                                                                                                       ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                                >  1:2655500/1‑150 (MQ=255)
                                                                                                                        tcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCt                               <  1:3553803/150‑1 (MQ=255)
                                                                                                                          ttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCttt                             >  1:585080/1‑150 (MQ=255)
                                                                                                                          ttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCttt                             <  1:4599544/150‑1 (MQ=255)
                                                                                                                          ttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCttt                             <  1:3284954/150‑1 (MQ=255)
                                                                                                                          ttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCttt                             <  1:2261682/150‑1 (MQ=255)
                                                                                                                            cGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGa                           <  1:1079623/150‑1 (MQ=255)
                                                                                                                             gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                          >  1:306548/1‑150 (MQ=255)
                                                                                                                             gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                          >  2:4578289/1‑150 (MQ=255)
                                                                                                                             gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                          <  1:871648/150‑1 (MQ=255)
                                                                                                                              cTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGc                         >  1:1346630/1‑150 (MQ=255)
                                                                                                                                gTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCtt                       <  1:674702/150‑1 (MQ=255)
                                                                                                                                 tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTg                      >  2:3981679/1‑150 (MQ=255)
                                                                                                                                 tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTg                      <  2:585080/150‑1 (MQ=255)
                                                                                                                                   gAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc                    >  1:5483376/1‑150 (MQ=255)
                                                                                                                                   gAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc                    >  2:3426788/1‑150 (MQ=255)
                                                                                                                                    aaaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGcc                   >  2:3132486/1‑150 (MQ=255)
                                                                                                                                     aaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCt                  >  2:1971478/1‑150 (MQ=255)
                                                                                                                                     aaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCt                  >  2:51434/1‑150 (MQ=255)
                                                                                                                                      aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 <  1:4275636/150‑1 (MQ=255)
                                                                                                                                      aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 <  1:2570521/150‑1 (MQ=255)
                                                                                                                                      aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 >  1:2972387/1‑150 (MQ=255)
                                                                                                                                        tcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCtt                              >  1:5326629/4‑135 (MQ=255)
                                                                                                                                        tcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCtt                              <  2:5326629/132‑1 (MQ=255)
                                                                                                                                         gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              >  1:1945901/1‑150 (MQ=255)
                                                                                                                                         gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              >  2:3186783/1‑150 (MQ=255)
                                                                                                                                         gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTACACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              >  1:1703409/1‑150 (MQ=255)
                                                                                                                                          cccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt             >  1:1522789/1‑150 (MQ=255)
                                                                                                                                          accAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt             >  1:2894751/2‑150 (MQ=255)
                                                                                                                                           ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAtt            <  1:5570128/150‑1 (MQ=255)
                                                                                                                                             aGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt          >  2:336188/1‑150 (MQ=255)
                                                                                                                                             aGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt          <  2:624062/150‑1 (MQ=255)
                                                                                                                                             aGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt          >  2:5011691/1‑150 (MQ=255)
                                                                                                                                               gctgctGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgc        >  1:416465/1‑150 (MQ=255)
                                                                                                                                                     ggcggcCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGaa  >  2:121527/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGAA  >  NC_002947/1071142‑1071440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑