Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,616 C→G 5.1% 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 5.1% 438.9 / 10.5 148G145A (GGG→GCG) yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (67/71);  new base G (5/3);  total (72/76)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 153 total aligned reads displayed.

TCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAAC  >  NC_000913/2449467‑2449763
                                                                                                                                                      |                                                                                                                                                   
tcatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcc                                                                                                                                                     <  2:619527/150‑1 (MQ=255)
  atcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcccc                                                                                                                                                   >  1:901538/1‑150 (MQ=255)
  atcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcccc                                                                                                                                                   >  1:1826236/1‑150 (MQ=255)
   tcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCt                                                                                                                                                  <  2:1977516/150‑1 (MQ=255)
    caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAACCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                                                                 >  1:1786036/1‑150 (MQ=255)
    caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                                                                 >  1:1670416/1‑150 (MQ=255)
     aaaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGAAATTATTCTCATTACATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCGCCTGa                                                                                                                                                >  2:121835/1‑150 (MQ=255)
      aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                                                                               <  1:1416562/150‑1 (MQ=255)
      aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                                                                               >  2:1049431/1‑150 (MQ=255)
      aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                                                                               <  2:1765365/150‑1 (MQ=255)
          cGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCAGATGGTTGAGCTGAACCGGTTGGACAAATCGTGGCGTGAGTATCGGGCGCCTGATACgg                                                                                                                                           >  2:1959532/1‑149 (MQ=255)
                 atcGTTTTCTCAACCACACCGCACGTCGACTTCAGGGTGATCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGTCGTATCAGGGGCCTGATAGTCACgaca                                                                                                                                    <  2:1140962/148‑4 (MQ=255)
                 aGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCCTCTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                                                                    >  2:128016/1‑150 (MQ=255)
                 aGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATACCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgaca                                                                                                                                    <  1:699809/150‑4 (MQ=255)
                  gCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATACGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                                                                   <  1:415589/150‑5 (MQ=255)
                        tctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                                                             <  2:495932/150‑1 (MQ=255)
                         ctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGGCATAAGGTTGAGCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgagc                                                                                                                            >  2:1316281/1‑149 (MQ=255)
                           cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                                                          >  1:1661952/1‑150 (MQ=255)
                             accaccCCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                        <  2:1487010/150‑1 (MQ=255)
                             aCCACACCGCCCGTCGACTGCATGGTGAGCAGTTTTGACTTACTGACATAATTCTCTTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                        <  1:933575/150‑1 (MQ=255)
                             aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAGTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                                        >  1:324178/1‑150 (MQ=255)
                                 cacCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                                                                    <  2:1884476/150‑1 (MQ=255)
                                         gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                            <  1:1176353/150‑1 (MQ=255)
                                         gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                            <  1:1986406/150‑1 (MQ=255)
                                         gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                            >  1:362951/1‑150 (MQ=255)
                                           cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACTGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                          <  1:219817/150‑1 (MQ=255)
                                           cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                          <  1:1871324/150‑1 (MQ=255)
                                            gACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCata                                                                                                         >  1:1936983/1‑150 (MQ=255)
                                             aCTGCAGTGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                                                                                        <  2:1826235/150‑1 (MQ=255)
                                              cTGCATGTTGAGCAGTTTTGACTTACTGACATAATTCTCATACAATGAAACGTCATACGGTTGAACAGAAACGGATGGACAAATCGGGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGGCATTCataat                                                                                                       <  2:1119087/150‑1 (MQ=255)
                                                  aGGGTGAGCAGTTTTTACTTACTGACATAATTCTCATTCCATGAAACGGCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGGCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATCATcacc                                                                                                   <  2:144706/150‑1 (MQ=255)
                                                    ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaa                                                                                                 >  2:914251/1‑150 (MQ=255)
                                                         gCAGTTTTGACTTTCTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                            <  1:986694/150‑1 (MQ=255)
                                                         gCAGTTTTGAATTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                            <  2:166343/150‑1 (MQ=255)
                                                         ccAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                            <  1:128099/149‑1 (MQ=255)
                                                           aGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATAGTTTAGTGACTATCAGGCCCCTGTAAGGTGACGCCCGCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                                                          <  1:252050/150‑1 (MQ=255)
                                                            gTTTTGACTTACTGACATAATTCTTATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                                                         <  1:2041314/150‑1 (MQ=255)
                                                            gTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                                                         <  1:420379/150‑1 (MQ=255)
                                                            gTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGATGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                                                         <  2:1150651/150‑1 (MQ=255)
                                                            gTTTTGACTAACTGACATCATGCTCATTCCAGGAAACGTCCTAAGGTTGAACAGAAACGGATGGACAAATCGGGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATCGAc                                                                                         <  2:448467/150‑1 (MQ=255)
                                                                 gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                                                                    >  1:1508613/1‑150 (MQ=255)
                                                                   cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                                                                  >  1:1768139/1‑150 (MQ=255)
                                                                    ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGTCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                                                                 <  1:372668/150‑1 (MQ=255)
                                                                        tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                                                             >  1:528562/1‑150 (MQ=255)
                                                                         gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                            >  2:1856797/1‑150 (MQ=255)
                                                                         gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGGGTCGTGACTATCTGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                            <  2:1278750/150‑1 (MQ=255)
                                                                           cATAATTCTCATTCCCTGAAACGTCATAAGGTTGAACAGAAACGGATGGACATATAGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                          <  2:1372896/150‑1 (MQ=255)
                                                                           cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACCAGAAAAGTAAGAtt                                                                          >  1:90394/1‑150 (MQ=255)
                                                                           cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAACAGTAAGCtt                                                                          >  1:1213084/1‑150 (MQ=255)
                                                                           cATAATTCTCATACCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGTCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                          >  1:331311/1‑150 (MQ=255)
                                                                           aaTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                          <  2:851605/149‑1 (MQ=255)
                                                                            aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtta                                                                         >  2:682629/1‑149 (MQ=255)
                                                                             tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                                                        >  1:1392791/1‑150 (MQ=255)
                                                                              aaTTCTCATTCCATGAAACGTCATAAGGTTGACCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCACTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                       >  1:1226102/1‑150 (MQ=255)
                                                                              aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                       <  1:2021756/150‑1 (MQ=255)
                                                                                 tctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtca                                                                    >  2:727772/1‑149 (MQ=255)
                                                                                  ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtctc                                                                   >  2:1735000/1‑150 (MQ=255)
                                                                                    cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATGGACAAGAAAAGAAAGATTCAATCTCTg                                                                 <  2:1383286/150‑1 (MQ=255)
                                                                                      ttCCATGAAACGTCATAAGGTTGAACAGAAACGTATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                               <  2:1261540/150‑1 (MQ=255)
                                                                                      ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                               <  2:26682/150‑1 (MQ=255)
                                                                                        ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                                                             <  2:501806/150‑1 (MQ=255)
                                                                                         cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                            >  1:1153524/1‑150 (MQ=255)
                                                                                          aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga                                                           >  2:538612/1‑150 (MQ=255)
                                                                                           tGAAACGTCATAAGGTTGAACAGAAACGGATGGCCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgagaa                                                          >  2:1084655/1‑149 (MQ=255)
                                                                                            gAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGGCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAAAAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                                                         <  2:1889555/150‑1 (MQ=255)
                                                                                              aaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt                                                       <  1:1263796/150‑1 (MQ=255)
                                                                                                cGTCATAAGGTTGAACAGAAACGGATGGAAAAATCGTGTCTTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa                                                     <  1:1363285/150‑1 (MQ=255)
                                                                                                 gTCATAAGGTTGAACAGAAACGGATGGACACATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCCATCTCTGCAGTGACAATAAc                                                    >  1:1004595/1‑150 (MQ=255)
                                                                                                 gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                                    >  1:1113508/1‑150 (MQ=255)
                                                                                                  tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                   >  1:677409/1‑150 (MQ=255)
                                                                                                   cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTTACTATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt                                                  <  1:285129/150‑1 (MQ=255)
                                                                                                     tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                                >  1:276973/1‑150 (MQ=255)
                                                                                                        ggTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAAt                                             <  1:1980230/150‑1 (MQ=255)
                                                                                                         gctcaaCAGAAACGGCTGGACAACTCGGCTCGTGACTATCCGCCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                            <  2:1369518/146‑1 (MQ=255)
                                                                                                         gTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                            <  2:1972338/150‑1 (MQ=255)
                                                                                                           tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                          <  1:1529722/150‑1 (MQ=255)
                                                                                                           tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                          <  1:1473228/150‑1 (MQ=255)
                                                                                                            gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                         <  2:430437/150‑1 (MQ=255)
                                                                                                             aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCAGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtata                                        >  2:350514/1‑150 (MQ=255)
                                                                                                               cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                                      <  2:362951/150‑1 (MQ=255)
                                                                                                                      ggATGGACAAATCGTGTCGTGACTATCAGGCCACTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCATATGTATAATACCCGCa                               >  2:212134/1‑150 (MQ=255)
                                                                                                                       gATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAg                              <  2:1121387/150‑1 (MQ=255)
                                                                                                                         tGGACAAATCGTGTCGTGACTATCAGTAACATGATACTTCACGCTCACTGGAGAGTCTTCAGAGTCATTCATAATAACAAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCATATGTATAATACCCGCAGtt                            <  1:264718/150‑1 (MQ=255)
                                                                                                                             cAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAc                        >  2:1348630/1‑150 (MQ=255)
                                                                                                                               aaTCTTGTCGTGACTATCAGGCACCTGATACGTCACGCCCCCTGGAGATGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAccc                      <  1:142129/150‑1 (MQ=255)
                                                                                                                               aaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAATTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGGATAATACGCGCCGTTCCAacc                      >  2:371036/1‑150 (MQ=255)
                                                                                                                               aaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAccc                      <  1:1959533/150‑1 (MQ=255)
                                                                                                                                 tcgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCCACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAg                    >  1:1567919/1‑150 (MQ=255)
                                                                                                                                    gggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAATCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACg                 >  2:188797/4‑150 (MQ=255)
                                                                                                                                     gtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                <  2:1786035/150‑1 (MQ=255)
                                                                                                                                     ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                >  2:1878297/3‑150 (MQ=255)
                                                                                                                                     ggcgtgACTATCAGGCCCCTGATACGTCACGCCACCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg                <  2:331311/148‑1 (MQ=255)
                                                                                                                                       cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa              >  1:564601/1‑150 (MQ=255)
                                                                                                                                           aCTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATcc          >  1:1655405/1‑150 (MQ=255)
                                                                                                                                            cTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCa         >  2:861697/1‑150 (MQ=255)
                                                                                                                                              tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAg        <  1:1735001/150‑1 (MQ=255)
                                                                                                                                                tCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCAACCCAGACGGTAATACAGGa      >  2:1482466/1‑150 (MQ=255)
                                                                                                                                                 cAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGAc     <  1:914252/150‑1 (MQ=255)
                                                                                                                                                 cAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGATCCACCCAGACGGTAATCCAGGAc     >  2:1319450/1‑150 (MQ=255)
                                                                                                                                                   ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACaa   >  1:1601133/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAAC  >  NC_000913/2449467‑2449763

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: