Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,141,380 A→T 5.2% L281H (CTC→CAC) ‡ ptsA ← putative PTS multiphosphoryl transfer protein PtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,141,3800AT 5.2% 507.4 / 11.6 154L281H (CTC→CAC) ‡ptsAputative PTS multiphosphoryl transfer protein PtsA
Reads supporting (aligned to +/- strand):  ref base A (71/75);  new base T (6/2);  total (77/77)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 157 total aligned reads displayed.

AAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAGCTG  >  NC_000913/4141231‑4141529
                                                                                                                                                     |                                                                                                                                                     
aaaTCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGa                                                                                                                                                       <  1:1191181/150‑1 (MQ=255)
  aTCGTTTGTTGCAGCCACGGAGTAACGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAgg                                                                                                                                                     <  2:896237/150‑1 (MQ=255)
     gTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGGGTTGGAATGTTGCACGCACGGGCAAGGATCACCGTATGGGAGGTGGTGCCGCCGCTTTTGAGCACCAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGaag                                                                                                                                                  >  2:651806/1‑149 (MQ=255)
     gTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAAc                                                                                                                                                  >  2:658389/1‑150 (MQ=255)
     gTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAAc                                                                                                                                                  <  1:1241676/150‑1 (MQ=255)
     gTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAAc                                                                                                                                                  >  2:163164/1‑150 (MQ=255)
            gctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggctgg                                                                                                                                           <  1:130668/150‑1 (MQ=255)
             ctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggg                                                                                                                                          <  1:481534/150‑1 (MQ=255)
               gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGCACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                        >  1:813284/1‑150 (MQ=255)
               gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                        >  1:877194/1‑150 (MQ=255)
               gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                        >  2:701777/1‑150 (MQ=255)
               gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                        >  1:24207/1‑150 (MQ=255)
                 cACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTc                                                                                                                                      <  2:483303/150‑1 (MQ=255)
                   cGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAg                                                                                                                                    <  2:413083/150‑1 (MQ=255)
                      aGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtc                                                                                                                                 >  1:474744/1‑150 (MQ=255)
                       caaaGGGCATCAATATCCACACCAACCAGCGTTGGAATCTTGAACGAACCGGCAAGGATCACCGTATCTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtca                                                                                                                                >  2:498750/3‑150 (MQ=255)
                        tAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcat                                                                                                                               <  1:959703/150‑1 (MQ=255)
                         aaGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcatc                                                                                                                              <  1:160599/150‑1 (MQ=255)
                           gggCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAg                                                                                                                            >  2:248555/1‑150 (MQ=255)
                             gCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGcc                                                                                                                          <  2:839922/150‑1 (MQ=255)
                               aTCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCAt                                                                                                                        <  2:552535/150‑1 (MQ=255)
                                 cAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATAc                                                                                                                      <  2:1060672/150‑1 (MQ=255)
                                 cAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATAc                                                                                                                      <  1:100869/150‑1 (MQ=255)
                                       ccACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATgg                                                                                                                >  2:769246/1‑150 (MQ=255)
                                       ccACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATgg                                                                                                                >  1:840465/1‑150 (MQ=255)
                                            ccaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGc                                                                                                           >  1:998474/1‑150 (MQ=255)
                                             caaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGGCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCt                                                                                                          <  2:971211/150‑1 (MQ=255)
                                              aaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGA‑TTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCCGCCATACAAATGGCGGGCTGc                                                                                                        >  2:1140419/1‑150 (MQ=255)
                                                 caGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGt                                                                                                      >  2:785984/1‑150 (MQ=255)
                                                  aGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTATGCGGTGAGTTTTATCCAGTTCGAGGAACTGGCTGCGGGTCAGTTCATCAGCCATACAAATCGCGGGCTGCGTc                                                                                                     >  2:1424247/1‑150 (MQ=255)
                                                       tGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGttt                                                                                                <  2:1303004/150‑1 (MQ=255)
                                                           aTGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGccc                                                                                            <  2:300105/150‑1 (MQ=255)
                                                             gTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCgg                                                                                          <  2:951517/150‑1 (MQ=255)
                                                                gaacgaacGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGc                                                                                       >  1:966670/1‑150 (MQ=255)
                                                                   cgaacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGGTGGGTGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCgg                                                                                    <  2:950091/150‑1 (MQ=255)
                                                                    gaacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCggg                                                                                   <  2:468101/150‑1 (MQ=255)
                                                                     aacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGa                                                                                  <  2:314403/150‑1 (MQ=255)
                                                                      acgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGaa                                                                                 <  1:163164/150‑1 (MQ=255)
                                                                      acgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGaa                                                                                 >  1:813611/1‑150 (MQ=255)
                                                                         ggCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCAGGGAAGCg                                                                              >  2:1040172/1‑150 (MQ=255)
                                                                          gCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGt                                                                             >  1:748012/1‑150 (MQ=255)
                                                                            aaGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTg                                                                           <  2:1144383/150‑1 (MQ=255)
                                                                             aGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGc                                                                          <  2:1048287/150‑1 (MQ=255)
                                                                             aGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGc                                                                          <  1:994540/150‑1 (MQ=255)
                                                                             aGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGc                                                                          >  1:130097/1‑150 (MQ=255)
                                                                             aGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGc                                                                          <  2:268523/150‑1 (MQ=255)
                                                                                 tCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCAc                                                                      >  2:51648/1‑150 (MQ=255)
                                                                                       tATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGa                                                                >  1:1386985/1‑150 (MQ=255)
                                                                                          gtgAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGAttt                                                             >  1:765607/1‑150 (MQ=255)
                                                                                          gtgAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGAttt                                                             >  1:683086/1‑150 (MQ=255)
                                                                                            gAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGc                                                           <  1:1411058/150‑1 (MQ=255)
                                                                                            gAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAAAGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGc                                                           >  2:1384180/1‑150 (MQ=255)
                                                                                              ggtggtGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTg                                                         <  2:50082/150‑1 (MQ=255)
                                                                                              ggtggtGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTg                                                         <  1:1234430/150‑1 (MQ=255)
                                                                                                 ggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAg                                                      <  1:1057656/150‑1 (MQ=255)
                                                                                                   tgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTa                                                    >  2:1372240/1‑150 (MQ=255)
                                                                                                    gccgccGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTaa                                                   >  1:1403268/1‑150 (MQ=255)
                                                                                                      cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                                 >  2:270430/1‑150 (MQ=255)
                                                                                                      cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                                 >  2:419185/1‑150 (MQ=255)
                                                                                                      cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                                 >  1:104556/1‑150 (MQ=255)
                                                                                                      cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGAGTTGCTGGAGTAACt                                                 >  1:578333/1‑150 (MQ=255)
                                                                                                       gccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTg                                                >  1:1154906/1‑150 (MQ=255)
                                                                                                         cgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGa                                              <  2:877194/150‑1 (MQ=255)
                                                                                                         cgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGa                                              <  1:403314/150‑1 (MQ=255)
                                                                                                          gcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGaa                                             <  1:419185/150‑1 (MQ=255)
                                                                                                           cTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGTAGTAACTGGAAg                                            <  2:813611/150‑1 (MQ=255)
                                                                                                               tGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGa                                        <  1:1339091/150‑1 (MQ=255)
                                                                                                                    aacaaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcg                                   <  2:557219/150‑1 (MQ=255)
                                                                                                                      caaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcgcg                                 <  2:1184102/150‑1 (MQ=255)
                                                                                                                      caaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcgcg                                 <  2:130097/150‑1 (MQ=255)
                                                                                                                       aaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGt                                >  1:875246/1‑150 (MQ=255)
                                                                                                                       aaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGt                                <  1:677125/150‑1 (MQ=255)
                                                                                                                             ttGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGAGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTcc                          >  1:1111126/1‑150 (MQ=255)
                                                                                                                                aGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAgg                       <  1:979372/150‑1 (MQ=255)
                                                                                                                                 ggTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAggg                      >  2:408177/1‑150 (MQ=255)
                                                                                                                                   tGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAg                                                                 <  1:479936/105‑1 (MQ=255)
                                                                                                                                   tGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAg                                                                 >  2:479936/1‑105 (MQ=255)
                                                                                                                                   tGATTTTTATCCACTTCGATTAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATCGAGGGCTGCGTCAGTTTGCCCGGTGACGCGAAGCGTTGATCACAGTAGATTTGCTGGAGTAACTTGAAGCAGACGACCCGTACGTCCAGGCCa                    >  2:357573/1‑150 (MQ=255)
                                                                                                                                      tttttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGt                 >  1:273612/1‑150 (MQ=255)
                                                                                                                                      cccccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGt                 >  2:890021/8‑150 (MQ=255)
                                                                                                                                         ccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCt              >  1:941480/5‑150 (MQ=255)
                                                                                                                                         ccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCt              >  1:69573/5‑150 (MQ=255)
                                                                                                                                         ccaTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCt              >  1:783890/3‑150 (MQ=255)
                                                                                                                                          tATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCtt             >  2:507516/1‑150 (MQ=255)
                                                                                                                                           aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTg            <  1:886324/150‑1 (MQ=255)
                                                                                                                                           aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTg            >  1:108571/1‑150 (MQ=255)
                                                                                                                                            tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc           <  1:1017722/150‑1 (MQ=255)
                                                                                                                                            tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc           <  2:638297/150‑1 (MQ=255)
                                                                                                                                            gccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc           <  2:476292/149‑1 (MQ=255)
                                                                                                                                             ccAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCa          <  2:1126877/150‑1 (MQ=255)
                                                                                                                                              cAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAg         >  2:235433/1‑150 (MQ=255)
                                                                                                                                              cAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAg         <  2:1404227/150‑1 (MQ=255)
                                                                                                                                              cAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAg         >  2:865605/1‑150 (MQ=255)
                                                                                                                                               aGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAgg        >  1:655757/1‑150 (MQ=255)
                                                                                                                                                   cGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgc    >  1:856120/1‑150 (MQ=255)
                                                                                                                                                    gAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgct   >  1:148686/1‑150 (MQ=255)
                                                                                                                                                    gAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgct   >  1:222314/1‑150 (MQ=255)
                                                                                                                                                    gAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgct   >  1:914516/1‑150 (MQ=255)
                                                                                                                                                    gAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgct   >  1:47389/1‑150 (MQ=255)
                                                                                                                                                     aGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgctg  >  1:1159122/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAGCTG  >  NC_000913/4141231‑4141529

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