Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 26,844 A→T 6.1% intergenic (+40/+145) PP_0021 → / ← PP_0022 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294726,8440AT 6.1% 575.9 / 23.0 198intergenic (+40/+145)PP_0021/PP_0022hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (105/81);  new base T (8/4);  total (113/85)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

Only 100 of 205 total aligned reads displayed.

GGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCTA  >  NC_002947/26695‑26991
                                                                                                                                                      |                                                                                                                                                   
  ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAgg                                                                                                                                                   >  1:286169/1‑149 (MQ=255)
  ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                                   >  2:2314857/1‑150 (MQ=255)
  ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAACTGCGATACGCTTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAgg                                                                                                                                                   >  2:64062/1‑149 (MQ=255)
       aTCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGttt                                                                                                                                              >  2:368072/1‑150 (MQ=255)
         cccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTAc                                                                                                                                            <  1:1472282/150‑1 (MQ=255)
          ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGttaca                                                                                                                                            <  1:1723004/150‑4 (MQ=255)
          ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                           >  1:1322392/1‑150 (MQ=255)
           cTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAg                                                                                                                                          >  1:218670/1‑150 (MQ=255)
            ttATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                         >  1:619890/1‑150 (MQ=255)
            ttATCCGGATGCCGGTTGATCATACTTGTCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                         >  1:1639752/1‑150 (MQ=255)
                  ggATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaa                                                                                                                                   >  2:824505/1‑150 (MQ=255)
                  ggATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaa                                                                                                                                   <  2:1424091/150‑1 (MQ=255)
                    aTGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                                 <  2:195567/150‑1 (MQ=255)
                    aTGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                                 <  2:265904/150‑1 (MQ=255)
                         ggTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTggg                                                                                                                            >  1:2298016/1‑149 (MQ=255)
                          gTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTg                                                                                                                                                  <  1:1829904/127‑1 (MQ=255)
                               tcatcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCCGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACaa                                                                                                                      >  1:1258433/1‑150 (MQ=255)
                               tcatcCTTATCGCAGAATGGATAGAGTGAAGATGGCCAGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACaa                                                                                                                      >  2:2197273/1‑150 (MQ=255)
                                     ttATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                                 >  1:2057092/1‑133 (MQ=255)
                                     ttATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGc                                                                                                                >  2:1988027/1‑150 (MQ=255)
                                        tCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                             >  2:1855419/1‑150 (MQ=255)
                                        tCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                             >  2:1950161/1‑150 (MQ=255)
                                          gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGa                                                                                                           >  2:1669631/1‑150 (MQ=255)
                                          gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTAACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAATGGTTTACAGTGAAAAAAAGAGGTTTACAACATAGCTAAGa                                                                                                           >  2:30007/1‑150 (MQ=38)
                                           cAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTgg                                                                                                                                                 >  2:2176903/1‑111 (MQ=255)
                                               aTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACAcgt                                                                                                      <  1:585489/150‑2 (MQ=255)
                                               aTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGGGGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACAcgc                                                                                                      <  1:1005234/150‑1 (MQ=255)
                                                      gagTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               <  2:2442894/150‑1 (MQ=255)
                                                       ggggAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGTCTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                              <  1:388663/147‑1 (MQ=255)
                                                       agTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGCCTACGGCGTTGGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                              <  2:1639752/150‑1 (MQ=38)
                                                         tGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCa                                                                                            <  2:2207534/150‑1 (MQ=255)
                                                              aTGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTCTTGTAAACCACTCTGCTGCACCGTAACCACAGTGGTTTACAGGGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTccc                                                                                       >  2:1263219/1‑150 (MQ=255)
                                                                   ccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTCACCACAGTGGTTTACAGTGAAAACAAGAGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                                  >  2:2335559/1‑150 (MQ=255)
                                                                   ccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                                  <  2:481819/150‑1 (MQ=255)
                                                                    cGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCTCTGTAACTACAGTGGTTTACAGTGAAAAAAATTGGTTTCAAACGTAGCTACGACACGCGCCCTGGCCACTTCCCAGGcc                                                                                 >  1:915967/1‑150 (MQ=39)
                                                                           aaCTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAgcgcgc                                                                          >  2:196726/1‑150 (MQ=255)
                                                                             cTGTGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAATCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGCTTTACACCGTAGCCACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctc                                                                        <  2:1016772/150‑1 (MQ=255)
                                                                              tGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCAAAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctcc                                                                       >  2:1112385/1‑150 (MQ=255)
                                                                                     tCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGa                                                                <  2:1825228/150‑1 (MQ=255)
                                                                                        gctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTc                                                             <  1:1950161/150‑1 (MQ=255)
                                                                                          tgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTTGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGt                                                           >  2:209234/1‑150 (MQ=255)
                                                                                           ggAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                          >  1:2416678/1‑150 (MQ=255)
                                                                                           ggAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                          <  2:2043188/150‑1 (MQ=255)
                                                                                           ggAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGGGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCggt                                                          >  1:9782/1‑148 (MQ=255)
                                                                                              gCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTgg                                                                                                                             <  1:2248728/80‑1 (MQ=255)
                                                                                              gCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTgg                                                                                                                             >  2:2248728/1‑80 (MQ=255)
                                                                                               cTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGt                                                      >  2:121649/1‑150 (MQ=255)
                                                                                                 gACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTTTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCg                                                    >  2:293155/1‑150 (MQ=255)
                                                                                                    tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAGAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCAATCGTTATATCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                                 >  2:1380278/1‑150 (MQ=37)
                                                                                                    tCAGGCCTGAGGAGGCCTGTTGTAAACCACCTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCATGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                                 <  2:1280036/150‑1 (MQ=255)
                                                                                                         ccTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATGGATTCGTTCAGTCGGGCGTCGa                                            >  2:1197617/1‑150 (MQ=255)
                                                                                                           tGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAac                                          <  1:2216199/150‑1 (MQ=255)
                                                                                                           tGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAac                                          >  2:561912/1‑150 (MQ=255)
                                                                                                             aggaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAg                                                                                                                 <  2:803361/77‑1 (MQ=255)
                                                                                                             aggaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCGGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                        <  1:2319366/150‑1 (MQ=255)
                                                                                                              ggaggCCTGGTGTAAACCACTTTTTTGCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACAcc                                       <  2:569423/150‑1 (MQ=255)
                                                                                                               gaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCt                                      <  1:1491479/150‑1 (MQ=255)
                                                                                                                 ggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                                    <  1:209234/150‑1 (MQ=255)
                                                                                                                  gCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGaa                                           >  2:908554/1‑142 (MQ=255)
                                                                                                                    cTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgct                                 >  2:771244/1‑150 (MQ=255)
                                                                                                                     tgttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGTGCGCTCCTCAtt                                                                  <  1:2201013/116‑1 (MQ=255)
                                                                                                                      gttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAgcg                                                                             >  2:1778541/1‑104 (MQ=255)
                                                                                                                      gttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAgcg                                                                             <  1:1778541/104‑1 (MQ=255)
                                                                                                                       ttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTCCAGTGACATAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATgctgct                              >  1:2004871/1‑150 (MQ=255)
                                                                                                                         gtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCg                            >  1:1388035/1‑150 (MQ=255)
                                                                                                                          tAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                                  >  1:781231/1‑111 (MQ=255)
                                                                                                                          tAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgg                           >  1:874856/1‑150 (MQ=255)
                                                                                                                            aaCCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgggg                         >  1:899031/1‑150 (MQ=37)
                                                                                                                            aaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctc                                                                        >  2:1248433/1‑103 (MQ=255)
                                                                                                                            aaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgggg                         <  1:1070186/150‑1 (MQ=255)
                                                                                                                            aaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgggg                         >  2:1580960/1‑150 (MQ=255)
                                                                                                                             aCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                          >  1:1235855/1‑116 (MQ=255)
                                                                                                                             aCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                          <  2:1235855/116‑1 (MQ=255)
                                                                                                                             aCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGc                        >  1:1119017/1‑150 (MQ=255)
                                                                                                                               cACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               <  1:1889174/77‑1 (MQ=14)
                                                                                                                               cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGTTACTACACGCGCCCTGGAAACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGAGTCGAACACCTATGCTTCTCGGGGCCt                      >  1:705002/1‑150 (MQ=38)
                                                                                                                               cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCa                                                                                >  1:1415743/1‑92 (MQ=255)
                                                                                                                               cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCgg                                                   >  1:2408120/1‑121 (MQ=255)
                                                                                                                               cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCt                      >  2:2348383/1‑150 (MQ=255)
                                                                                                                               cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCt                      >  2:192395/1‑150 (MQ=255)
                                                                                                                                aCTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGt                                                           >  1:81455/1‑112 (MQ=255)
                                                                                                                                 cTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTg                    <  1:835182/150‑1 (MQ=255)
                                                                                                                                  tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                                  <  1:857857/103‑1 (MQ=255)
                                                                                                                                  tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                                  >  2:857857/1‑103 (MQ=255)
                                                                                                                                  tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                              <  2:71951/123‑1 (MQ=255)
                                                                                                                                  tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                              >  1:71951/1‑123 (MQ=255)
                                                                                                                                      tttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcg               >  2:416262/1‑150 (MQ=255)
                                                                                                                                       ttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCgg                                                   <  2:1622270/113‑1 (MQ=255)
                                                                                                                                       ttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCgg                                                   >  1:1622270/1‑113 (MQ=255)
                                                                                                                                        tCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgcc             >  2:1892148/1‑150 (MQ=255)
                                                                                                                                        tCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgcc             >  1:15852/1‑150 (MQ=255)
                                                                                                                                          aCTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccgc           >  2:1155464/1‑150 (MQ=38)
                                                                                                                                            tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                        <  1:1219578/119‑1 (MQ=255)
                                                                                                                                            tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                        >  2:1219578/1‑119 (MQ=255)
                                                                                                                                                 aCCACAGTGGTTTACAGTGAAAAAAAGTGGTTTAAAGCGTCGCTTTTACACGGGCTCTGGCCACTCCCCAGGCCGTCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTc     <  1:532958/150‑1 (MQ=255)
                                                                                                                                                  ccACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACATAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgtg       >  1:2120427/1‑146 (MQ=255)
                                                                                                                                                   cacaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCt   >  1:753006/1‑150 (MQ=255)
                                                                                                                                                    acaGTGGTTTACAGTGAAAACAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCa                                                                                            <  1:1495401/60‑1 (MQ=255)
                                                                                                                                                    acaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCa                                                                                            >  2:1495401/1‑60 (MQ=255)
                                                                                                                                                    acaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCta  >  2:2063515/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCTA  >  NC_002947/26695‑26991

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