Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 585,866 G→A 5.4% intergenic (‑233/‑191) ompT ← / → pauD DLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,8660GA 5.4% 409.1 / 6.9 129intergenic (‑233/‑191)ompT/pauDDLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base G (58/62);  new base A (5/2);  total (63/66)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 137 total aligned reads displayed.

GTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAG  >  NC_000913/585718‑586015
                                                                                                                                                    |                                                                                                                                                     
gTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTagag                                                                                                                                                     <  2:901429/150‑1 (MQ=255)
    atatTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTATCTATTTAAGTGCCATGCCTTTAGAGGCat                                                                                                                                                 <  1:1091180/150‑1 (MQ=255)
      atTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                                <  1:711399/150‑1 (MQ=255)
      atTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                                >  2:1354457/1‑150 (MQ=255)
        tttGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCAtata                                                                                                                                              >  2:412905/1‑150 (MQ=255)
           ttttAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                           <  2:1117466/149‑1 (MQ=255)
           gTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                           <  1:554746/150‑1 (MQ=255)
           gTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                           >  2:148284/1‑150 (MQ=255)
           gTTTAGTATTGAAAGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                           <  2:1339194/150‑1 (MQ=255)
            tttAGTATTGAATGATTATAACAATGCGGTATCTGTTTTTGGTTTAATTATCTGTTTTTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCg                                                                                                                                          <  1:173879/150‑1 (MQ=255)
             ttAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATATGCGa                                                                                                                                         >  1:174164/1‑150 (MQ=255)
               aGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaa                                                                                                                                       >  1:1562572/1‑150 (MQ=255)
                gTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaaa                                                                                                                                      >  1:934707/1‑150 (MQ=255)
                   ttGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAg                                                                                                                                   <  2:902965/150‑1 (MQ=255)
                     gAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATTTTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGTTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa                                                                                                                                 <  1:1618541/150‑1 (MQ=255)
                     gAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGATTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa                                                                                                                                 >  1:1113902/1‑150 (MQ=255)
                         gATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAg                                                                                                                             >  2:564465/1‑150 (MQ=255)
                           ttATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGt                                                                                                                           >  1:1657292/1‑150 (MQ=255)
                             aTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGttt                                                                                                                         >  1:151295/1‑150 (MQ=255)
                              tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTa                                                                                                                        >  1:1407951/1‑150 (MQ=255)
                                acaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATc                                                                                                                      <  1:1038743/150‑1 (MQ=255)
                                acaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATc                                                                                                                      <  2:641237/150‑1 (MQ=255)
                                 caATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTTTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATcc                                                                                                                     <  1:656892/150‑1 (MQ=255)
                                 caATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATcc                                                                                                                     <  1:1680799/150‑1 (MQ=255)
                                   aTGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTc                                                                                                                   >  2:1369109/1‑150 (MQ=255)
                                      cgcTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATTCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAAt                                                                                                                <  1:844007/150‑1 (MQ=255)
                                             tGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATg                                                                                                         <  2:505838/150‑1 (MQ=255)
                                                ttttGGTTTAATTTTCTGTTTTTGTTTCATATTTCGGTTTTTCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTTTCTATTTAATTGCCATTTCTTTAGAGTCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGttt                                                                                                      <  1:261302/150‑1 (MQ=255)
                                                  ttGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGttttt                                                                                                    <  2:1593985/150‑1 (MQ=255)
                                                    ggTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                                  >  1:434850/1‑150 (MQ=255)
                                                      tttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTAAATTACTATGTTTTTTAGt                                                                                                >  2:887406/1‑150 (MQ=255)
                                                      tttAATTATCTGTAATAGGTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTTAGTGCCATGCCTTTAGATGCATATAAGCGAAAATATCAAGAGGTTTATCCTCAATTACTATGTTTTTTAGt                                                                                                <  1:312582/150‑1 (MQ=255)
                                                         aaTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACa                                                                                             <  2:1286752/150‑1 (MQ=255)
                                                          aTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaa                                                                                            >  2:1047859/1‑150 (MQ=255)
                                                            tATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaa                                                                                          >  2:1645567/1‑150 (MQ=255)
                                                            tATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaa                                                                                          >  2:349715/1‑150 (MQ=255)
                                                                 gTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                     >  2:569199/1‑150 (MQ=255)
                                                                  ttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAgg                                                                                    <  2:919834/150‑1 (MQ=255)
                                                                   tATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTTCTTTTATCTATTTAAGTGCCATGCCTTTAGATGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAggg                                                                                   <  1:440358/150‑1 (MQ=255)
                                                                   tATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAggg                                                                                   <  2:1618560/150‑1 (MQ=255)
                                                                      tGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACa                                                                                <  2:1227862/150‑1 (MQ=255)
                                                                      tGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACa                                                                                <  2:636211/150‑1 (MQ=255)
                                                                      tGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACa                                                                                <  2:625867/150‑1 (MQ=255)
                                                                       ttttCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACaa                                                                               <  1:549559/149‑1 (MQ=255)
                                                                              atTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                        <  1:126886/150‑1 (MQ=255)
                                                                               tttCGGTTTTGCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgaga                                                                       <  1:636932/150‑1 (MQ=255)
                                                                               tttCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgaga                                                                       <  1:746362/150‑1 (MQ=255)
                                                                                ttCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATAACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTTAGac                                                                      >  1:598656/1‑150 (MQ=255)
                                                                                 tCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGaca                                                                     >  1:382697/1‑150 (MQ=255)
                                                                                  cGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTTGTGCTTATATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGacac                                                                    <  1:850422/150‑1 (MQ=255)
                                                                                   ggTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGAcaca                                                                   >  1:1080357/1‑150 (MQ=255)
                                                                                    gTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCAATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACAt                                                                  <  2:665323/150‑1 (MQ=255)
                                                                                        tACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAgg                                                              <  2:382697/150‑1 (MQ=255)
                                                                                        tACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCAGATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAgg                                                              >  2:848581/1‑150 (MQ=255)
                                                                                         aCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGc                                                             >  2:538357/1‑150 (MQ=255)
                                                                                          cTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGcc                                                            <  2:787338/150‑1 (MQ=255)
                                                                                          cTGTGAGGTTTTATTATGCTATTGTGGCGCTTTTATCAATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGcc                                                            <  2:44446/150‑1 (MQ=255)
                                                                                            gtgtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTc                                                          >  1:1016934/1‑150 (MQ=255)
                                                                                                 gTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAAt                                                     >  1:600974/1‑150 (MQ=255)
                                                                                                  tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATg                                                    <  2:1362806/150‑1 (MQ=255)
                                                                                                  tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATg                                                    <  1:1369109/150‑1 (MQ=255)
                                                                                                  tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATg                                                    >  1:556759/1‑150 (MQ=255)
                                                                                                   ttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGAGAAAATAGCATGAGGTTTATCATCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAAttg                                                   >  1:1290792/1‑148 (MQ=255)
                                                                                                    tttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGc                                                  <  2:1015694/150‑1 (MQ=255)
                                                                                                     ttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCt                                                 <  1:968633/150‑1 (MQ=255)
                                                                                                            cTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                        <  1:361101/120‑1 (MQ=255)
                                                                                                            cTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAgg                                          <  1:1139022/150‑1 (MQ=255)
                                                                                                            cTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAgg                                          >  1:1458755/1‑150 (MQ=255)
                                                                                                               ttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCtt                                       <  1:1001992/150‑1 (MQ=255)
                                                                                                                tgttgtGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCttt                                      <  1:24373/150‑1 (MQ=255)
                                                                                                                tgtgGTGCTTTTATCTATTTAAGTGCCATTCCTTTTGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAATAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCttt                                      <  2:929271/150‑1 (MQ=255)
                                                                                                                tgtgGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCttt                                      <  2:151295/150‑1 (MQ=255)
                                                                                                                 gtggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTa                                     >  2:450329/1‑150 (MQ=255)
                                                                                                                  tggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTAc                                    >  2:854600/1‑150 (MQ=255)
                                                                                                                  tggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTAc                                    >  2:383557/1‑150 (MQ=255)
                                                                                                                   ggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACa                                   >  2:471159/1‑150 (MQ=255)
                                                                                                                   ggtgCTTTTATCTAATTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATATCATGATGTTTATCCTCAATTACTATGTTTTTTATTACAAAAAAGAGGTACAAAACTGAGACACATAAGGCCTCTCAATGACTTGCAAGGCATTACa                                   >  1:106979/1‑150 (MQ=255)
                                                                                                                     tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATg                                 >  1:872583/1‑150 (MQ=255)
                                                                                                                       cTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGAGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGtt                               >  1:1655708/1‑150 (MQ=255)
                                                                                                                          ttATCTTTTTAAGTGCCATTCCTTTATATGGATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTg                            <  2:241100/150‑1 (MQ=255)
                                                                                                                          ttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttttt                            >  1:958602/1‑149 (MQ=255)
                                                                                                                          ttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTg                            <  2:173447/150‑1 (MQ=255)
                                                                                                                          ttATATATTTCAGTGCCATGCCTTTAGAGGCATATAAGCGAAAAGAGCCTGAGGGTTATCCTCAATTAATATGTTTTTTAGTACAAAAAAGAGGGACAAAAATGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTGTTg                            <  2:311176/150‑1 (MQ=255)
                                                                                                                            aTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAg                          <  2:485791/150‑1 (MQ=255)
                                                                                                                            aTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAg                          <  1:1098142/150‑1 (MQ=255)
                                                                                                                                    cgcttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGATTTAAATGTTTTTAGTTAGTggg                  >  2:117259/8‑150 (MQ=255)
                                                                                                                                     aGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGa                 >  2:1200653/1‑150 (MQ=255)
                                                                                                                                      cttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGAc                >  1:165165/6‑150 (MQ=255)
                                                                                                                                       ttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGAGGGAct               >  2:35105/5‑149 (MQ=255)
                                                                                                                                       tGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACg               <  2:181134/150‑1 (MQ=255)
                                                                                                                                        gCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt              >  1:491111/1‑150 (MQ=255)
                                                                                                                                        gCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt              >  1:1505965/1‑150 (MQ=255)
                                                                                                                                         ccATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGTGACAAAACTGAGACTCATAAGGGCTCGCTATGGCTTGCAAGGCTTTTCATGTTTTGAGGTAGTGGGACgtg             <  1:1316966/150‑1 (MQ=255)
                                                                                                                                          cATGCATTTAGAGGAATATAAGCTAAAATAGCATGATGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACTAAACTGAGACAGATAATGCCTCGCAATGGCTTGAAATGCTTTACATGTTTTGAGGTAGTGGGACgtgt            <  2:1290794/150‑1 (MQ=255)
                                                                                                                                                 ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCa     >  2:232311/1‑150 (MQ=255)
                                                                                                                                                  tAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCag    >  2:85643/1‑150 (MQ=255)
                                                                                                                                                   agagGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCaga   <  2:1187374/150‑1 (MQ=255)
                                                                                                                                                   agagGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACGATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATATGGCCTCGCAATGGCTTGCAAGCCTTTAAATGTTTTGAGTTACTGGGACGTGTGAGCGCaga   >  1:89220/1‑150 (MQ=255)
                                                                                                                                                   agagGAATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTGTTTTAGTACAAAAAAGAGGGACAAAACTGAGACCCATACGACCTCGCAATGGCTTGAAACGCTTGACATTTTTTTAGGTATTGGGACTTTTGAGCGCata   >  2:131311/1‑148 (MQ=255)
                                                                                                                                                    gagGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCagag  >  2:1575742/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAG  >  NC_000913/585718‑586015

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