Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,071,291 G→A 5.1% 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.1% 2157.8 / 60.0 666A87A (GCC→GCTPP_0924hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (301/330);  new base A (17/17);  total (318/348)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 712 total aligned reads displayed.

GCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGAA  >  NC_002947/1071142‑1071440
                                                                                                                                                     |                                                                                                                                                     
gCCAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGctgctg                                                                                                                                                       >  1:1953745/1‑150 (MQ=255)
  cAGTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTggc                                                                                                                                                     >  2:1715576/1‑150 (MQ=255)
       cGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCt                                                                                                                                                <  2:7000099/150‑1 (MQ=255)
        gCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTg                                                                                                                                               <  1:8262324/150‑1 (MQ=255)
        gCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTg                                                                                                                                               <  1:3680128/150‑1 (MQ=255)
         ccccACAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                              <  2:5663272/150‑1 (MQ=255)
         ccccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTgg                                                                                                                                              >  2:7143503/1‑150 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCATGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTggg                                                                                                                                             <  2:6580047/150‑2 (MQ=255)
          cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGt                                                                                                                                             >  2:3046136/1‑150 (MQ=255)
           ccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTg                                                                                                                                            >  1:2117053/1‑150 (MQ=255)
           ccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTg                                                                                                                                            >  1:5071385/1‑150 (MQ=255)
           ccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTg                                                                                                                                            >  1:3356893/1‑150 (MQ=255)
           ccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgggc                                                                                                                                            <  1:935111/150‑3 (MQ=255)
            cAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGa                                                                                                                                           <  1:6910485/150‑1 (MQ=255)
              aaGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                         <  1:2617384/150‑1 (MQ=255)
              aaGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                         <  2:2733872/150‑1 (MQ=255)
                gCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGa                                                                                                                                       <  2:5869030/150‑1 (MQ=255)
                gCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGa                                                                                                                                       <  2:495805/150‑1 (MQ=255)
                 cACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCCGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAg                                                                                                                                      <  2:7606652/150‑1 (MQ=255)
                  aCGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGc                                                                                                                                     <  1:4328549/150‑1 (MQ=255)
                  aCGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGc                                                                                                                                     <  2:1059216/150‑1 (MQ=255)
                    ggCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAg                                                                                                                                   >  1:3858395/1‑150 (MQ=255)
                    ggCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAg                                                                                                                                   >  1:4942275/1‑150 (MQ=255)
                            tACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                           <  1:3191981/150‑1 (MQ=255)
                            tACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                           >  1:8677888/1‑150 (MQ=255)
                            tACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                           <  2:3598774/150‑1 (MQ=255)
                             aCCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTTCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGACTGCTGGCGGCCTGGTGATCGAGCAGGAGggtggt                                                                                                                          <  2:7431746/150‑1 (MQ=255)
                             aCCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCCGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtggt                                                                                                                          <  2:8657761/150‑1 (MQ=255)
                             aCCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtggt                                                                                                                          >  2:2269250/1‑150 (MQ=255)
                                gccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGc                                                                                                                       >  2:6147770/1‑150 (MQ=255)
                                   gTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAAc                                                                                                                    <  2:14928/150‑1 (MQ=255)
                                        ggTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGtt                                                                                                               >  1:5847881/1‑150 (MQ=255)
                                        ggTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGtt                                                                                                               >  2:23691/1‑150 (MQ=255)
                                           aGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGc                                                                                                            <  2:5774503/150‑1 (MQ=255)
                                            gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGcc                                                                                                           >  1:7969938/1‑150 (MQ=255)
                                            gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGcc                                                                                                           <  2:6380267/150‑1 (MQ=255)
                                             tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCg                                                                                                          <  1:4570019/150‑1 (MQ=255)
                                               gggggCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGt                                                                                                        >  1:6060484/1‑150 (MQ=255)
                                                  ggCGACGCAGGTACGGCCCAGTTCGAGGATTGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCt                                                                                                     >  2:5436738/1‑150 (MQ=255)
                                                   gCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTc                                                                                                    >  2:1408312/1‑150 (MQ=255)
                                                          aGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcgc                                                                                             <  1:6143156/150‑1 (MQ=255)
                                                          aGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcgc                                                                                             <  1:6531863/150‑1 (MQ=255)
                                                                ggCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       <  1:1715576/150‑1 (MQ=255)
                                                                ggCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       >  1:7887486/1‑150 (MQ=255)
                                                                ggCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       >  1:3565385/1‑150 (MQ=255)
                                                                  cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCa                                                                                     >  2:3257946/1‑150 (MQ=255)
                                                                  cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCa                                                                                     >  2:3815040/1‑150 (MQ=255)
                                                                   ccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCaa                                                                                    >  2:3776906/1‑150 (MQ=255)
                                                                   ccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCaa                                                                                    >  1:5560243/1‑150 (MQ=255)
                                                                   ccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCaa                                                                                    >  2:2518257/1‑150 (MQ=255)
                                                                      gTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAAtgt                                                                                 <  1:8574700/150‑1 (MQ=255)
                                                                       ttCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAAtgtg                                                                                <  2:6012911/150‑1 (MQ=255)
                                                                         cGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggc                                                                              <  2:1850302/150‑1 (MQ=255)
                                                                           aGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggcgg                                                                            <  2:526862/150‑1 (MQ=255)
                                                                            ggATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggcggc                                                                           >  2:4830351/1‑150 (MQ=255)
                                                                             gATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCa                                                                          >  2:8575949/1‑150 (MQ=255)
                                                                             gATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCa                                                                          <  1:1456264/150‑1 (MQ=255)
                                                                              aTCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAg                                                                         >  1:6753881/1‑150 (MQ=255)
                                                                               tCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGt                                                                        <  1:4446370/150‑1 (MQ=255)
                                                                                     cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcgt                                                                  >  1:286131/1‑150 (MQ=255)
                                                                                     cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcgt                                                                  <  2:6084031/150‑1 (MQ=255)
                                                                                         cAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                              <  1:1541685/150‑1 (MQ=255)
                                                                                         cAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                              <  1:728126/150‑1 (MQ=255)
                                                                                                aaCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCATACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTcc                                                       <  2:4111515/150‑1 (MQ=255)
                                                                                                   aGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATg                                                    >  2:557260/1‑150 (MQ=255)
                                                                                                   aGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATg                                                    <  1:569323/150‑1 (MQ=255)
                                                                                                     gCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTc                                                  >  1:2247736/1‑150 (MQ=255)
                                                                                                         tGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAgg                                              <  2:2604159/150‑1 (MQ=255)
                                                                                                          gCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGc                                             <  2:3989522/150‑1 (MQ=255)
                                                                                                            aGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGccc                                           >  1:1607164/1‑150 (MQ=255)
                                                                                                            aGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGccc                                           <  1:1001621/150‑1 (MQ=255)
                                                                                                              gCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTg                                         >  2:7245307/1‑150 (MQ=255)
                                                                                                                tGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCt                                       >  2:8095713/1‑150 (MQ=255)
                                                                                                                       ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                                >  1:7161678/1‑150 (MQ=255)
                                                                                                                         cttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCtt                              <  1:8452148/150‑1 (MQ=255)
                                                                                                                            cGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGa                           <  1:7368817/150‑1 (MQ=255)
                                                                                                                             gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                          >  1:5265553/1‑150 (MQ=255)
                                                                                                                               tGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCt                        >  1:7600330/1‑150 (MQ=255)
                                                                                                                                  aGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCCATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCGTTGAGCTTgg                     >  2:6681012/1‑150 (MQ=255)
                                                                                                                                   gAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc                    >  1:5181074/1‑150 (MQ=255)
                                                                                                                                    aaaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGcc                   <  1:1081918/150‑1 (MQ=255)
                                                                                                                                     aaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCt                  <  2:2117053/150‑1 (MQ=255)
                                                                                                                                      aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 >  1:4164457/1‑150 (MQ=255)
                                                                                                                                      aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 >  1:771832/1‑150 (MQ=255)
                                                                                                                                        ggCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGa               <  2:3412527/150‑1 (MQ=255)
                                                                                                                                         gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              <  2:7869522/150‑1 (MQ=255)
                                                                                                                                         caccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              >  2:4262433/3‑150 (MQ=255)
                                                                                                                                           ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAtt            <  2:1954023/150‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAtt            >  1:8085755/1‑150 (MQ=255)
                                                                                                                                            cAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTc           >  1:2644609/1‑150 (MQ=255)
                                                                                                                                            cAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTc           >  2:1252552/1‑150 (MQ=255)
                                                                                                                                             aGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt          <  1:6389331/150‑1 (MQ=255)
                                                                                                                                              ggCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTg         >  2:5101266/1‑150 (MQ=255)
                                                                                                                                              ggCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTg         <  1:1304514/150‑1 (MQ=255)
                                                                                                                                              ggCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTg         <  2:8343322/150‑1 (MQ=255)
                                                                                                                                               gctgctGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgc        <  2:2916583/150‑1 (MQ=255)
                                                                                                                                                ctgctgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgcg       <  1:238741/150‑1 (MQ=255)
                                                                                                                                                 tgctgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgcgc      >  1:5494885/1‑150 (MQ=255)
                                                                                                                                                    tgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGa                           <  2:4546376/126‑1 (MQ=255)
                                                                                                                                                    tgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGa   <  1:4783138/150‑1 (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: