Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 226,663 G→A 6.5% intergenic (+70/‑63) PP_0174 → / → PP_0175 hypothetical protein/transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947226,6630GA 6.5% 923.3 / 35.7 294intergenic (+70/‑63)PP_0174/PP_0175hypothetical protein/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (158/116);  new base A (8/11);  total (167/127)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 313 total aligned reads displayed.

AGGAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTG  >  NC_002947/226514‑226812
                                                                                                                                                     |                                                                                                                                                     
  gAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTgacgac                                                                                                                                                     >  1:3096282/1‑150 (MQ=255)
    gTGGCGGGACCGGGCGGTGCAGTTGGGCAGCACATTGGCCTTGACCATGAATCTGGAAGACCAGATCCATCCCTGAGTATTTCTCGTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAg                                                                                                                                                   <  1:2285969/150‑1 (MQ=255)
     tGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAgg                                                                                                                                                  <  1:4277258/150‑1 (MQ=255)
      ggCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTCGTCAACAggg                                                                                                                                                 <  1:3616475/150‑1 (MQ=255)
       gCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAgggg                                                                                                                                                <  2:6452868/150‑1 (MQ=255)
       gCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAgggg                                                                                                                                                >  1:4844103/1‑150 (MQ=255)
        cGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAggggg                                                                                                                                               >  1:190502/1‑150 (MQ=255)
              cGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAATACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAAAAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACTAAAGGGGGATGTCg                                                                                                                                         >  1:5396896/1‑150 (MQ=255)
               gggCGGTGCAGCTGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCgg                                                                                                                                        <  2:4746150/150‑1 (MQ=255)
                ggCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCggg                                                                                                                                       >  2:6264304/1‑150 (MQ=255)
                    gTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGt                                                                                                                                   <  1:5217360/150‑1 (MQ=255)
                       cAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGtata                                                                                                                                <  2:4216495/150‑1 (MQ=255)
                        aGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTatat                                                                                                                               <  1:2182408/150‑1 (MQ=255)
                             ggCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATa                                                                                                                          <  2:1415766/150‑1 (MQ=255)
                               cagcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGt                                                                                                                        >  1:3284525/1‑150 (MQ=255)
                                 gcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTa                                                                                                                      <  2:2273189/150‑1 (MQ=255)
                                 gcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTa                                                                                                                      <  1:910871/150‑1 (MQ=255)
                                     aTTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCaa                                                                                                                  >  1:1794483/1‑150 (MQ=255)
                                      ttGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCaaa                                                                                                                 <  1:7304720/150‑1 (MQ=255)
                                        ggCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGt                                                                                                                                           >  1:1101265/1‑122 (MQ=255)
                                          ccTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTAt                                                                                                             >  2:2950026/1‑150 (MQ=255)
                                           cTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATg                                                                                                            <  1:6346772/150‑1 (MQ=255)
                                            ttGCCCATGAAGCTGGAAGACCAGATCCATCCCGGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGa                                                                                                           <  2:427801/150‑1 (MQ=255)
                                            ttGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGa                                                                                                           >  1:2925468/1‑150 (MQ=255)
                                               aCCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAata                                                                                                        <  2:2121910/150‑1 (MQ=255)
                                                 cATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAAtata                                                                                                      >  1:6993347/1‑150 (MQ=255)
                                                 cATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAAtata                                                                                                      >  2:4906803/1‑150 (MQ=255)
                                                  aTGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATAc                                                                                                     >  2:5297618/1‑150 (MQ=255)
                                                      aGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTgg                                                                                                 <  1:6157222/150‑1 (MQ=255)
                                                        cTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCCGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGtt                                                                                               >  1:157118/1‑150 (MQ=255)
                                                           gAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTtctc                                                                                            <  2:104905/150‑1 (MQ=255)
                                                            aaGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTctct                                                                                           >  2:3801511/1‑150 (MQ=255)
                                                             aGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTGGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCtt                                                                                          >  2:3044384/1‑150 (MQ=255)
                                                                 cAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCTTg                                                                                      >  2:2224723/1‑150 (MQ=255)
                                                                    atccatccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTcc                                                                                   >  2:550860/1‑150 (MQ=255)
                                                                     tccatccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTccc                                                                                  <  1:4627248/150‑1 (MQ=255)
                                                                      ccatccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCt                                                                                 >  1:1463521/1‑150 (MQ=255)
                                                                      ccatccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCt                                                                                 >  1:4397835/1‑150 (MQ=255)
                                                                        atccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTg                                                                               >  2:5739092/1‑150 (MQ=255)
                                                                        atccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTCGTCAACAgg                                                                                                                                                  >  1:1035744/1‑83 (MQ=255)
                                                                         tccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGa                                                                              <  2:6100429/150‑1 (MQ=255)
                                                                          cccTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGAt                                                                             <  1:1576317/150‑1 (MQ=255)
                                                                             tGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTc                                                                          <  1:4918831/150‑1 (MQ=255)
                                                                              gAGTATTTCTCGTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTcc                                                                         <  1:3357464/150‑1 (MQ=255)
                                                                               aGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCt                                                                        <  2:3151571/150‑1 (MQ=255)
                                                                               aGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCt                                                                        <  1:851112/150‑1 (MQ=255)
                                                                                gTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTg                                                                       <  2:3100222/150‑1 (MQ=255)
                                                                                   tttCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTgcgc                                                                    >  1:6740304/1‑150 (MQ=255)
                                                                                   tttCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTgcgc                                                                    >  1:4725626/1‑150 (MQ=255)
                                                                                    ttCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGcc                                                                   >  2:3593435/1‑150 (MQ=255)
                                                                                         ttCAGTGCTGACTGAAAAATGAACAACATTTAAACAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCGCtt                                                                                >  1:4244177/1‑132 (MQ=255)
                                                                                            aGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTc                                                           >  1:5063972/1‑150 (MQ=255)
                                                                                             gTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTca                                                          <  1:3513019/150‑1 (MQ=255)
                                                                                              tGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTcag                                                         >  2:2856460/1‑150 (MQ=255)
                                                                                                ctgactgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTcagca                                                       <  1:2612854/150‑1 (MQ=255)
                                                                                                     tgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGc                                                  <  1:3912022/150‑1 (MQ=255)
                                                                                                     tgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGc                                                  >  1:1342227/1‑150 (MQ=255)
                                                                                                     tgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGc                                                  >  2:1261301/1‑150 (MQ=255)
                                                                                                     tgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGc                                                  >  2:6715184/1‑150 (MQ=255)
                                                                                                       aaaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCt                                                >  2:1563668/1‑150 (MQ=255)
                                                                                                       aaaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCt                                                >  2:6038829/1‑150 (MQ=255)
                                                                                                       aaaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCt                                                <  1:4544699/150‑1 (MQ=255)
                                                                                                        aaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTg                                               <  1:120705/150‑1 (MQ=255)
                                                                                                        aaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTg                                               <  1:957744/150‑1 (MQ=255)
                                                                                                         aaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTgg                                              >  2:3394080/1‑150 (MQ=255)
                                                                                                           aTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTa                                            <  1:2817370/150‑1 (MQ=255)
                                                                                                           aTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTa                                            <  2:465388/150‑1 (MQ=255)
                                                                                                                caacaTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGc                                       >  2:4494736/1‑150 (MQ=255)
                                                                                                                  acaTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTg                                     >  1:5595879/1‑150 (MQ=255)
                                                                                                                    aTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGcc                                   >  2:350133/1‑150 (MQ=255)
                                                                                                                     tttAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGccc                                  <  1:1677007/150‑1 (MQ=255)
                                                                                                                           cAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACt                            >  1:3462440/1‑150 (MQ=255)
                                                                                                                           cAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACt                            >  1:2105550/1‑150 (MQ=255)
                                                                                                                            aGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTg                           >  2:6191426/1‑150 (MQ=255)
                                                                                                                             ggCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTgg                          >  1:596990/1‑150 (MQ=255)
                                                                                                                             ggCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTgg                          >  2:579606/1‑150 (MQ=255)
                                                                                                                             ggCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGACAGCAGTGGCCTGGTAGTCGCTGCCCGCCACggg                          >  1:5417949/1‑150 (MQ=255)
                                                                                                                              gCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTggc                         >  1:5487903/1‑150 (MQ=255)
                                                                                                                               cTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTggcg                        <  2:1713142/150‑1 (MQ=255)
                                                                                                                                tATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTggcgg                       >  2:6155743/1‑150 (MQ=255)
                                                                                                                                 aTCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTggcggc                      >  2:7191974/1‑150 (MQ=255)
                                                                                                                                  tCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGgcggcg                     >  2:824493/1‑150 (MQ=255)
                                                                                                                                  tCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGgcggcg                     >  1:6884018/1‑150 (MQ=255)
                                                                                                                                       atcccccTGTCGTCAACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCt                >  2:2771445/3‑150 (MQ=255)
                                                                                                                                        tcccccTGTCGTCAACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTg               >  2:1962833/2‑150 (MQ=255)
                                                                                                                                         cccccTGTTTACGACAGGTGGATGTCGGGGCGTATATTGATAGTTAGCAAAATATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTTGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGt              <  2:5027095/150‑1 (MQ=255)
                                                                                                                                         cccccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGt              >  2:1637669/1‑150 (MQ=255)
                                                                                                                                          ccccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCCCTGGCGGCGCATCTGTc             >  1:2501688/1‑150 (MQ=255)
                                                                                                                                           cccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCCTCTGTCa            >  1:739726/1‑150 (MQ=255)
                                                                                                                                           cccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCa            >  1:403150/1‑150 (MQ=255)
                                                                                                                                            ccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCAt           >  2:1181078/1‑150 (MQ=255)
                                                                                                                                              tGTTGACGCCAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGGTCTCTTCATGTCCCTTGAGTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATAc         <  2:6904117/150‑1 (MQ=255)
                                                                                                                                              tGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATAc         <  2:6562493/150‑1 (MQ=255)
                                                                                                                                                 tGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGc      <  2:4078815/150‑1 (MQ=255)
                                                                                                                                                 tGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTCGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGc      <  2:282465/150‑1 (MQ=255)
                                                                                                                                                 agacgacAGGGGGATGTCGGTGCTTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATCTCCCTCGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATAcctc      <  2:3858727/149‑3 (MQ=255)
                                                                                                                                                  gacgacAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGcc     <  2:4700162/150‑1 (MQ=255)
                                                                                                                                                  gacgacAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGcc     >  2:5641623/1‑150 (MQ=255)
                                                                                                                                                     gacAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTg  >  2:108496/1‑150 (MQ=255)
                                                                                                                                                     gacAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTg  >  2:7429445/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTG  >  NC_002947/226514‑226812

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: