Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,280 G→A 5.9% intergenic (+15/‑37) efp → / → ecnA protein chain elongation factor EF‑P/entericidin A lipoprotein, antidote to entericidin B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,2800GA 5.9% 876.2 / 23.1 272intergenic (+15/‑37)efp/ecnAprotein chain elongation factor EF‑P/entericidin A lipoprotein, antidote to entericidin B
Reads supporting (aligned to +/- strand):  ref base G (129/125);  new base A (7/9);  total (137/135)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GATACCGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTC  >  NC_000913/4376131‑4376428
                                                                                                                                                     |                                                                                                                                                    
  tACCGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGc                                                                                                                                                    <  1:1541939/150‑1 (MQ=255)
  tACCGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGc                                                                                                                                                    >  1:774174/1‑150 (MQ=255)
   aCCGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGcc                                                                                                                                                   <  1:1848107/150‑1 (MQ=255)
    ccGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCt                                                                                                                                                  >  1:1890011/1‑150 (MQ=255)
      gCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGc                                                                                                                                                <  2:1586161/150‑1 (MQ=255)
          gTACTGGGGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGc                                                                                                                                            <  2:850286/150‑1 (MQ=255)
           tACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCt                                                                                                                                           >  1:1187556/1‑150 (MQ=255)
             cTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGc                                                                                                                                         <  2:349786/150‑1 (MQ=255)
              tggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCa                                                                                                                                        <  2:625791/150‑1 (MQ=255)
              tggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCa                                                                                                                                        <  1:916788/150‑1 (MQ=255)
               ggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCAc                                                                                                                                       >  2:1145509/1‑150 (MQ=255)
               ggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCAc                                                                                                                                       <  2:1081952/150‑1 (MQ=255)
               ggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCAc                                                                                                                                       >  1:1388972/1‑150 (MQ=255)
               ggtggCAAACCGGCTACCCTGTCTACTGGCGATGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCAc                                                                                                                                       >  2:429440/1‑150 (MQ=255)
                gtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCAcc                                                                                                                                      <  2:1731706/150‑1 (MQ=255)
                 tggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCACCa                                                                                                                                     <  2:1912672/150‑1 (MQ=255)
                    cAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCa                                                                                                                                  >  2:408158/1‑150 (MQ=255)
                     aaaCCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCACCacaac                                                                                                                                 <  2:688198/150‑5 (MQ=255)
                         cGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTAt                                                                                                                             >  1:350627/1‑150 (MQ=255)
                          ggCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTAtt                                                                                                                            >  1:1427754/1‑150 (MQ=255)
                          ggCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTAtt                                                                                                                            >  1:466417/1‑150 (MQ=255)
                            cTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCa                                                                                                                                  <  2:1124613/142‑1 (MQ=255)
                            cTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCa                                                                                                                                  >  1:1124613/1‑142 (MQ=255)
                            cTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                                          >  2:154176/1‑150 (MQ=255)
                                 cTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCa                                                                                                                     >  2:125746/1‑150 (MQ=255)
                                    tCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCagag                                                                                                                  <  2:578174/150‑1 (MQ=255)
                                       aCTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGc                                                                                                                                                <  2:613132/117‑1 (MQ=255)
                                         tGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAgatgat                                                                                                             <  2:1187556/150‑1 (MQ=255)
                                         tGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAgatgat                                                                                                             <  1:610364/150‑1 (MQ=255)
                                            cgcTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGc                                                                                                                                            >  2:1314532/1‑116 (MQ=255)
                                                  ggTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCt                                                                                                    >  2:550449/1‑150 (MQ=255)
                                                  ggTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCt                                                                                                    >  2:582827/1‑150 (MQ=255)
                                                       aaGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCg                                                                                               <  2:1991935/150‑1 (MQ=255)
                                                       aaGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCg                                                                                               >  1:28872/1‑150 (MQ=255)
                                                             cGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTg                                                                                         >  2:379502/1‑150 (MQ=255)
                                                                tGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGttt                                                                                      <  2:1012529/150‑1 (MQ=255)
                                                                tGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGttt                                                                                      <  2:106095/150‑1 (MQ=255)
                                                                 gTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGtttt                                                                                     >  1:1350605/1‑150 (MQ=255)
                                                                      tACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgcttg                                                                                <  2:1399423/150‑1 (MQ=255)
                                                                      tACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgcttg                                                                                <  1:814286/150‑1 (MQ=255)
                                                                       aCAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgcttgc                                                                               >  1:1004754/1‑150 (MQ=255)
                                                                       aCAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgcttgc                                                                               <  1:1793119/150‑1 (MQ=255)
                                                                        cAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCTGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGcc                                                                              <  1:679542/150‑1 (MQ=255)
                                                                        cAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGcc                                                                              >  2:1800950/1‑150 (MQ=255)
                                                                        cAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGcc                                                                              >  1:444759/1‑150 (MQ=255)
                                                                        cAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGcc                                                                              >  2:1653555/1‑150 (MQ=255)
                                                                          aaTCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTTCCAg                                                                            <  2:517426/150‑1 (MQ=255)
                                                                          aaTCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAg                                                                            <  1:1868486/150‑1 (MQ=255)
                                                                           aTCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGc                                                                           <  2:28872/150‑1 (MQ=255)
                                                                              ggCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACg                                                                        >  2:326934/1‑150 (MQ=255)
                                                                                cGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgct                                                                      <  1:1913304/150‑1 (MQ=255)
                                                                                  aaGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgctgc                                                                    <  2:92248/150‑1 (MQ=255)
                                                                                  aaGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgctgc                                                                    <  2:466417/150‑1 (MQ=255)
                                                                                       atcaAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACg                                                               <  2:791039/150‑1 (MQ=255)
                                                                                         caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttg                                                                                    >  2:344812/1‑127 (MQ=255)
                                                                                         caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACggg                                                             >  1:1537041/1‑150 (MQ=255)
                                                                                          aaaGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGc                                                            >  2:767409/1‑150 (MQ=255)
                                                                                          aaaGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGc                                                            <  1:1580237/150‑1 (MQ=255)
                                                                                          aaaGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGc                                                            <  2:1354429/150‑1 (MQ=255)
                                                                                          aaaGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGc                                                            <  1:1318541/150‑1 (MQ=255)
                                                                                             gTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTACTCACGGGCTGt                                                         <  2:533973/150‑1 (MQ=255)
                                                                                              tGGATAGCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTa                                                        <  2:1563786/150‑1 (MQ=255)
                                                                                              tGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTa                                                        >  1:548868/1‑150 (MQ=255)
                                                                                              tGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTa                                                        <  1:379502/150‑1 (MQ=255)
                                                                                               ggATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTaa                                                       >  2:1925801/1‑150 (MQ=255)
                                                                                                     ccGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGc                                                 >  1:1693792/1‑150 (MQ=255)
                                                                                                       gCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTc                                               <  1:951320/150‑1 (MQ=255)
                                                                                                       gCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTc                                               >  1:644176/1‑150 (MQ=255)
                                                                                                           tGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCgg                                           <  1:51767/150‑1 (MQ=255)
                                                                                                           tGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCgg                                           <  1:326934/150‑1 (MQ=255)
                                                                                                           tGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCgg                                           <  1:1266210/150‑1 (MQ=255)
                                                                                                            ggTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGt                                          >  2:1478498/1‑150 (MQ=255)
                                                                                                                aaTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCg                                      >  2:1873777/1‑150 (MQ=255)
                                                                                                                 aTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCgg                                     <  1:582827/150‑1 (MQ=255)
                                                                                                                  tACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCGAGCGGtgttaca                                    <  1:73345/150‑7 (MQ=255)
                                                                                                                   aCGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCg                                   >  1:1196457/1‑150 (MQ=255)
                                                                                                                     gTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGaa                                 <  1:1983170/150‑1 (MQ=255)
                                                                                                                     gTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGaa                                 <  2:652911/150‑1 (MQ=255)
                                                                                                                      tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAg                                >  1:1483639/1‑150 (MQ=255)
                                                                                                                      tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAg                                <  1:1655794/150‑1 (MQ=255)
                                                                                                                           cgcgTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATc                           >  1:1404176/1‑150 (MQ=255)
                                                                                                                           cgcgTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTACACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATc                           <  1:149578/150‑1 (MQ=255)
                                                                                                                               tGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAAc                       <  2:1529881/150‑1 (MQ=255)
                                                                                                                                  aGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACAtc                    >  2:963307/1‑150 (MQ=255)
                                                                                                                                   gTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACAtct                   >  1:1740838/1‑150 (MQ=255)
                                                                                                                                    tAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACAtctc                  >  1:1074122/1‑150 (MQ=255)
                                                                                                                                       tGCGGTTGTGGTGCGGCCTGCAGGCTGTACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGc               >  2:257369/1‑150 (MQ=255)
                                                                                                                                       tGCGGTTGTGGGGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGc               <  1:1213726/150‑1 (MQ=255)
                                                                                                                                        gCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCa              >  1:1312445/1‑150 (MQ=255)
                                                                                                                                          ggTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAAc            <  2:1498798/150‑1 (MQ=255)
                                                                                                                                          ggTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAAc            <  1:429440/150‑1 (MQ=255)
                                                                                                                                              gtggtgCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCa        <  1:172264/150‑1 (MQ=255)
                                                                                                                                              atggtgCAGCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCa        >  2:1460881/2‑150 (MQ=255)
                                                                                                                                               tggtgCAGCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCAt       >  2:1214124/1‑150 (MQ=255)
                                                                                                                                                ggtgCAGCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCAtc      >  1:1037297/1‑150 (MQ=255)
                                                                                                                                                 gtgCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTcg                                              <  1:560165/109‑1 (MQ=255)
                                                                                                                                                  tgCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCAtctc    >  1:1625907/1‑150 (MQ=255)
                                                                                                                                                    cGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCAtctctc  >  2:1278103/1‑150 (MQ=255)
                                                                                                                                                     ggCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAAc                       >  2:1977932/1‑128 (MQ=255)
                                                                                                                                                     ggCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAAc                       <  1:1977932/128‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATACCGCAGGTACTGGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTC  >  NC_000913/4376131‑4376428

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: