Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,616 C→G 7.3% G145A (GGG→GCG) ‡ yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6160CG 7.3% 332.3 / 15.5 109G145A (GGG→GCG) ‡yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (57/44);  new base G (2/6);  total (59/50)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-02

Only 100 of 113 total aligned reads displayed.

CATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTT  >  NC_000913/2449468‑2449765
                                                                                                                                                     |                                                                                                                                                     
catcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGacc                                                                                                                                                      >  2:23109/1‑150 (MQ=255)
 atcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcccc                                                                                                                                                     >  1:1637280/1‑150 (MQ=255)
   caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTg                                                                                                                                                   >  1:620352/1‑150 (MQ=255)
    aaaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCGGa                                                                                                                                                  >  1:860575/1‑150 (MQ=255)
       tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATAc                                                                                                                                               <  1:766027/150‑1 (MQ=255)
        aCGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACg                                                                                                                                              >  1:355276/1‑150 (MQ=255)
          gTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATagtc                                                                                                                                             <  1:1357759/150‑4 (MQ=255)
          gTTGAAAGCGTTTTCTCAACCACACCCCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTc                                                                                                                                            <  1:627452/150‑1 (MQ=255)
           ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                                                                           <  1:1496385/150‑1 (MQ=255)
           ttGAAAGCGTTATCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGCATGGACAAATCGTGTCGTGAC‑TATCAGGGGCCTGATA‑GTCAc                                                                                                                                          <  2:333949/150‑1 (MQ=255)
            tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGAGATAATTCTCAATCCATGAAACGTCATAAGGTTGTACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                                                                          <  1:405203/150‑1 (MQ=255)
            tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGGGCCTGATA‑GTCACg                                                                                                                                         <  2:182970/150‑1 (MQ=255)
                 gCGTTTTCTCAAACACACCGCCCTTAGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                                                                     <  2:60811/150‑5 (MQ=255)
                    ttttCTCAACCACACCGACCGACGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacgat                                                                                                                                  <  2:1171898/150‑8 (MQ=38)
                       tctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                                                               <  2:1113982/150‑1 (MQ=255)
                        ctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                                                                              >  1:1103960/1‑150 (MQ=255)
                         tcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                                                             <  1:1702447/150‑1 (MQ=255)
                            aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                          <  2:926162/150‑1 (MQ=255)
                            aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                          >  2:771039/1‑150 (MQ=255)
                            aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                          >  1:1104220/1‑150 (MQ=255)
                            aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                          >  1:1110022/1‑150 (MQ=255)
                            aCCACACCGCCCGTCGACTGCAGGGGGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                          <  2:1087909/150‑1 (MQ=255)
                             ccACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAAGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                                                                         <  2:630743/150‑1 (MQ=255)
                               acacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                                                                       <  1:1176777/150‑1 (MQ=255)
                               acacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                                                                       >  2:364300/1‑150 (MQ=255)
                               acacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTACATGAAACGTCATAAGTTTGAACAGAAACGGAGGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                                                                       <  1:1122168/150‑1 (MQ=255)
                                cacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                                                                      <  2:531920/150‑1 (MQ=255)
                                  ccGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                    >  2:578170/1‑150 (MQ=255)
                                  ccGCCCGTCGACTGCAGGGGGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                    <  1:818701/150‑1 (MQ=255)
                                   cGCCCGTCGCCTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                   >  1:344683/1‑150 (MQ=255)
                                   cGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                   <  2:1637280/150‑1 (MQ=255)
                                   cGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                   >  1:1558809/1‑150 (MQ=255)
                                   cGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGAACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                   <  2:325167/150‑1 (MQ=255)
                                       cGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                                                               >  2:545450/1‑150 (MQ=255)
                                       cGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                                                               <  1:1377210/150‑1 (MQ=255)
                                          cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                            >  2:1419048/1‑150 (MQ=255)
                                          cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                            >  1:191633/1‑150 (MQ=255)
                                            aCTGCAGTGTTAGCAGTTTTTACTTACTGACATAATTCTCTTTCCATGAAACGTCATAAGGTTGAACAGTAACGGATGGACAAATCGTGTCGTGAC‑TATCATGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                                                                                          <  1:312099/150‑1 (MQ=255)
                                             cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                         >  1:148584/1‑150 (MQ=255)
                                                  gggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                                                                    >  1:511380/1‑150 (MQ=255)
                                                   ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCTTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTAAAGAGTCATTCATAATaacaaa                                                                                                   >  1:816581/1‑150 (MQ=255)
                                                       aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGGTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                                                                               <  1:1371609/150‑1 (MQ=255)
                                                        gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                              >  2:1682517/1‑150 (MQ=255)
                                                         ctgTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGGCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATgg                                                                                             <  2:1702068/148‑1 (MQ=255)
                                                          aGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                                                            <  1:1205388/150‑1 (MQ=255)
                                                           gTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                                                           >  1:207391/1‑150 (MQ=255)
                                                                      cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACAAACCATGGACAAGAAAAGTaa                                                                                >  1:417524/1‑150 (MQ=255)
                                                                        gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCTTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                              >  2:1466483/1‑150 (MQ=255)
                                                                        gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                              >  2:1473155/1‑150 (MQ=255)
                                                                        gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                              >  1:1307022/1‑150 (MQ=255)
                                                                         aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                                                                             >  2:770790/1‑150 (MQ=255)
                                                                          cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTATGAtt                                                                            >  1:1075366/1‑150 (MQ=255)
                                                                          cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                            >  1:445284/1‑150 (MQ=255)
                                                                           aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTc                                                                           <  2:1103960/150‑1 (MQ=255)
                                                                            tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                                                          <  1:545450/150‑1 (MQ=255)
                                                                             aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGCTGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                         <  2:323019/150‑1 (MQ=255)
                                                                              aTTCTCATTACATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                                                                        <  1:599530/150‑1 (MQ=255)
                                                                               ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                                                                       <  1:357667/150‑1 (MQ=255)
                                                                                   cATTCCATGAAACGTCATAATGTTGAACAGAACCGGATGGACACATCGTGTCGTGAC‑TATAAGGCCCATGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                   <  1:873036/150‑1 (MQ=255)
                                                                                   cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGACCCCTGGAGAGGCTTCAGAGTCATTCTTAATAACCATCCATGGTCAAGAAAAGTAAGATTCAATCTCTg                                                                   >  1:111179/1‑150 (MQ=255)
                                                                                    aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                                                                  >  2:1690235/1‑150 (MQ=255)
                                                                                     ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                                 <  2:445284/150‑1 (MQ=255)
                                                                                          tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                            <  2:1307022/150‑1 (MQ=255)
                                                                                           gAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                                                           <  1:1315333/150‑1 (MQ=255)
                                                                                             aaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt                                                         >  2:1419950/1‑150 (MQ=255)
                                                                                                  cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt                                                    >  2:1107870/1‑150 (MQ=255)
                                                                                                       ggTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACTTCAAAt                                               >  2:1337467/1‑150 (MQ=255)
                                                                                                        gTTGAACAGAAACGGATTGACAAATCTTGTCGTTAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACAATGTACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                              <  1:204470/150‑1 (MQ=255)
                                                                                                        gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                              >  1:234045/1‑150 (MQ=255)
                                                                                                          tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATACCCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCATAGACAATAACCGCAAAAGta                                            >  2:275937/1‑150 (MQ=255)
                                                                                                          tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                            <  1:69723/150‑1 (MQ=255)
                                                                                                          tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                            >  2:159606/1‑150 (MQ=255)
                                                                                                           taaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                           <  1:1682517/149‑1 (MQ=255)
                                                                                                              cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                                        >  2:727908/1‑150 (MQ=255)
                                                                                                               aGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataata                                       >  1:1347693/1‑150 (MQ=255)
                                                                                                                 aaaCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAcc                                     <  2:191633/150‑1 (MQ=255)
                                                                                                                  aaCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAccc                                    >  2:1308557/1‑150 (MQ=255)
                                                                                                                    cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGc                                  >  2:840041/1‑150 (MQ=255)
                                                                                                                      gATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAAGCAACCATGGACAAGAAAAGTAAGCTTCAAGCTCTGCAGAGACAATCACCTCATATGTATAATACACGCAg                                >  2:296605/1‑150 (MQ=255)
                                                                                                                         ggACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGTAGAGACAATAACCTCAAATGTATAATACCCGCAGTTc                             >  1:337526/1‑150 (MQ=255)
                                                                                                                         ggACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTc                             <  1:159606/150‑1 (MQ=255)
                                                                                                                         ggACAAATCGTGTCGTGAC‑TATCAGGACCCTGATACGTCACTCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTc                             <  1:275937/150‑1 (MQ=255)
                                                                                                                          gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATATCCAACCATGGACAAGAAAAGTAAGATTGAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCAGCAGTAcc                            >  1:264382/1‑150 (MQ=255)
                                                                                                                           aCAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCa                           <  2:1558809/150‑1 (MQ=255)
                                                                                                                            cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAc                          <  1:1559427/150‑1 (MQ=255)
                                                                                                                                    gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                  >  1:1327492/1‑150 (MQ=255)
                                                                                                                                    gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                  >  1:957732/1‑150 (MQ=255)
                                                                                                                                    ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                  >  1:753957/3‑150 (MQ=255)
                                                                                                                                     tcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGt                 <  2:148584/150‑1 (MQ=255)
                                                                                                                                       gtgacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTTAGATTCAATCTCTGCAGAGTCATTATGCTCAAATGTATTATACCCGCAGTTCCACCCAGACGGTag               >  2:289344/7‑149 (MQ=255)
                                                                                                                                        tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAAt              >  1:603176/1‑150 (MQ=255)
                                                                                                                                         gAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATc             <  2:620352/150‑1 (MQ=255)
                                                                                                                                             tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATATCCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAg          <  2:1467846/150‑1 (MQ=255)
                                                                                                                                             tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAAGCATGGACAACAAAAGGAAGATTCAATCTCTGCAGAGAAAATAACCTCAATTGTAAAATACCCGAAGTACCACACAGCAGGTAATCCAg          >  2:83106/1‑150 (MQ=255)
                                                                                                                                              aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTGAGAGTCATTCATAATAACCAACCATGTACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTTAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCTgg         >  2:292186/1‑150 (MQ=255)
                                                                                                                                              aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAgg         >  1:1616727/1‑150 (MQ=255)
                                                                                                                                               tCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGa        <  1:1364079/150‑1 (MQ=255)
                                                                                                                                               tCAGGCCCCTGAAACGGCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGa        <  1:295513/150‑1 (MQ=255)
                                                                                                                                                  ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACaa     >  2:619854/1‑150 (MQ=255)
                                                                                                                                                   gCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGCGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAAc    >  1:546099/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTT  >  NC_000913/2449468‑2449765

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