Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,223,929 T→C 5.3% intergenic (+22/+350) ycgI → / ← minE protein YcgI/Z‑ring positioning protein MinE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,223,9290TC 5.3% 468.1 / 11.0 153intergenic (+22/+350)ycgI/minEprotein YcgI/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base T (76/67);  new base C (5/3);  total (81/72)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

Only 100 of 154 total aligned reads displayed.

TATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACA  >  NC_000913/1223780‑1224077
                                                                                                                                                     |                                                                                                                                                    
  taCCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTc                                                                                                                                                    >  1:1006734/1‑150 (MQ=255)
   aCCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGGGGCTGGTGTGAACTGGTCGTGCTGACCATCAACGATTAAACTGCGCTTCg                                                                                                                                                   <  1:357163/150‑1 (MQ=255)
    ccAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTTGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCgg                                                                                                                                                  >  1:1589379/1‑150 (MQ=255)
    ccAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCgg                                                                                                                                                  >  1:631595/1‑150 (MQ=255)
           gCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtt                                                                                                                                           <  1:2389343/150‑1 (MQ=255)
           gCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtt                                                                                                                                           <  1:1038340/150‑1 (MQ=255)
            cTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGttt                                                                                                                                          >  1:2228808/1‑150 (MQ=255)
             tATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtttt                                                                                                                                         >  1:752177/1‑150 (MQ=255)
              aTTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTc                                                                                                                                        >  1:1378954/1‑150 (MQ=255)
              aTTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGAGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTc                                                                                                                                        >  2:2020079/1‑150 (MQ=255)
                tAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGt                                                                                                                                      >  1:762874/1‑150 (MQ=255)
                 aaaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGTGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGtt                                                                                                                                     <  1:451861/150‑1 (MQ=255)
                 aaaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGtt                                                                                                                                     >  1:409612/1‑150 (MQ=255)
                 aaaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtttaatcg                                                                                                                                     <  1:1034918/150‑6 (MQ=17)
                      cTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAg                                                                                                                                >  1:449636/1‑150 (MQ=255)
                           tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGt                                                                                                                           <  2:777991/150‑1 (MQ=255)
                            cagcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGtt                                                                                                                          >  1:1832746/1‑150 (MQ=255)
                             agcagTACGTTTAGCGGACATCTCGGGGGGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTg                                                                                                                         <  2:921310/150‑1 (MQ=255)
                                 gTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAg                                                                                                                     <  1:1962778/150‑1 (MQ=255)
                                 gTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAg                                                                                                                     <  1:1972527/150‑1 (MQ=255)
                                  tACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAgg                                                                                                                    <  2:1743844/150‑1 (MQ=255)
                                     gTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGaa                                                                                                                 <  1:726101/150‑1 (MQ=255)
                                           cGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCt                                                                                                           >  1:314486/1‑150 (MQ=255)
                                            ggACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTg                                                                                                          <  2:1365382/150‑1 (MQ=255)
                                               cATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTgggg                                                                                                       <  2:1772058/150‑1 (MQ=255)
                                               cATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTgggg                                                                                                       <  2:1464905/150‑1 (MQ=255)
                                                aTCTCGGGGTGGGGTATAGCCATGGTGGCGGTGGGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGTGGc                                                                                                      <  2:585790/150‑1 (MQ=255)
                                                 tctcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCa                                                                                                     >  1:380587/1‑150 (MQ=255)
                                                  ctcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaa                                                                                                    >  2:2265652/1‑150 (MQ=255)
                                                    cGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaaa                                                                                                  >  2:97272/1‑150 (MQ=255)
                                                    cGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGCaaaaa                                                                                                 >  1:1308629/1‑150 (MQ=255)
                                                       ggtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaaaaaa                                                                                               >  2:691509/1‑150 (MQ=255)
                                                        gtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAg                                                                                              >  1:159567/1‑150 (MQ=255)
                                                            ggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCg                                                                                          >  2:754109/1‑150 (MQ=255)
                                                              tataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGcc                                                                                        <  2:1814504/150‑1 (MQ=255)
                                                              tataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGcc                                                                                        >  2:2344173/1‑150 (MQ=255)
                                                              tataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGcc                                                                                        <  2:306408/150‑1 (MQ=255)
                                                               ataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCa                                                                                       >  1:2042866/1‑150 (MQ=255)
                                                                taGCCATGGGGCCGGTGTGGAATCCCGGTGGAAAGCGGTGGCTGGTGTGAACTGGGCGTTCTGACCATCAACGATTATACTGCGCGGCGGCGCAGTTTTCGTTTACAGGATGTGGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAg                                                                                      <  2:1576818/150‑1 (MQ=255)
                                                                  gCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGtt                                                                                    >  2:220446/1‑150 (MQ=255)
                                                                  gCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGtt                                                                                    >  1:1010347/1‑150 (MQ=255)
                                                                      tGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcgg                                                                                >  1:1196341/1‑150 (MQ=255)
                                                                       ggTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCGGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcggc                                                                               >  1:2262821/1‑150 (MQ=255)
                                                                        gTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcggcg                                                                              <  2:2472023/150‑1 (MQ=255)
                                                                          gCCGGTGTGGAATCCCCGTGGAACGCTGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCggg                                                                            <  2:433562/150‑1 (MQ=255)
                                                                          gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCggg                                                                            >  2:341807/1‑150 (MQ=255)
                                                                          gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCggg                                                                            <  1:960091/150‑1 (MQ=255)
                                                                           ccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGa                                                                           <  2:2153786/150‑1 (MQ=255)
                                                                                 gggAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTTGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTc                                                                     <  1:1750425/149‑1 (MQ=255)
                                                                                  ggAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCa                                                                    >  1:699439/1‑150 (MQ=255)
                                                                                  ggAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCa                                                                    >  2:1375047/1‑150 (MQ=255)
                                                                                       ccccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTa                                                               >  1:1240337/1‑150 (MQ=255)
                                                                                        cccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTAt                                                              >  1:1243874/1‑150 (MQ=255)
                                                                                        cccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTAt                                                              >  2:435146/1‑150 (MQ=255)
                                                                                          cGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATgg                                                            <  2:450781/150‑1 (MQ=255)
                                                                                            tGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGga                                                          >  2:1590508/1‑150 (MQ=255)
                                                                                              gAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgaga                                                        <  2:180577/150‑1 (MQ=255)
                                                                                               aaCGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGaa                                                       >  2:725858/1‑150 (MQ=255)
                                                                                                ccgcgGTGGCTGGGGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCTTTTACAGGATGTTGAAAGGGAAAATTCTGGGCCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAAc                                                      <  2:1326496/149‑1 (MQ=255)
                                                                                                acccgctGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCAAAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAAc                                                      >  2:274746/7‑150 (MQ=255)
                                                                                                aCGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAAc                                                      <  1:2344173/150‑1 (MQ=255)
                                                                                                  gcgGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACaa                                                    <  2:1918616/150‑1 (MQ=255)
                                                                                                  gcgGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACaa                                                    <  2:875470/150‑1 (MQ=255)
                                                                                                   cgGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAAt                                                   >  1:273960/1‑150 (MQ=255)
                                                                                                    ggTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATg                                                  <  1:2202783/150‑1 (MQ=255)
                                                                                                        gCTGGTGTTAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaata                                              <  2:2375751/150‑1 (MQ=255)
                                                                                                             tgtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaa                                         >  2:373914/1‑150 (MQ=255)
                                                                                                             tgtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaa                                         <  2:1783819/150‑1 (MQ=255)
                                                                                                              gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaaa                                        >  1:2472337/1‑150 (MQ=255)
                                                                                                              gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaaa                                        >  1:1181558/1‑150 (MQ=255)
                                                                                                               tgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAt                                       >  1:1686685/1‑150 (MQ=255)
                                                                                                                  aCTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGc                                    <  2:228297/150‑1 (MQ=255)
                                                                                                                   cTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCg                                   >  2:1420090/1‑150 (MQ=255)
                                                                                                                     ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCggg                                 <  1:2020079/150‑1 (MQ=255)
                                                                                                                     ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCggg                                 <  1:1117664/150‑1 (MQ=255)
                                                                                                                       tCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGt                               <  2:1859378/150‑1 (MQ=255)
                                                                                                                       tCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGt                               <  1:552444/150‑1 (MQ=255)
                                                                                                                           tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatc                           >  1:376315/1‑150 (MQ=255)
                                                                                                                            cTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatca                          >  1:182538/1‑150 (MQ=255)
                                                                                                                             tGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatcat                         >  2:1461021/1‑150 (MQ=255)
                                                                                                                              taCCATCAACGATTAAACTGCGCTTAGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatcatc                        <  2:2058828/149‑1 (MQ=255)
                                                                                                                              gACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatcatc                        <  2:808013/150‑1 (MQ=255)
                                                                                                                               aCCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtct                       <  2:122994/150‑1 (MQ=255)
                                                                                                                                ccATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtctc                      >  1:111197/1‑150 (MQ=255)
                                                                                                                                ccATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtctc                      >  1:373060/1‑150 (MQ=255)
                                                                                                                                  aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTccc                    <  2:752177/150‑1 (MQ=255)
                                                                                                                                   tCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCa                   <  1:1322206/150‑1 (MQ=255)
                                                                                                                                     aaCGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGGTATCATCTCCCAGt                 <  2:728186/150‑1 (MQ=255)
                                                                                                                                      aaacgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGta                >  2:1282660/7‑150 (MQ=255)
                                                                                                                                      aaacgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGta                >  1:535824/7‑150 (MQ=255)
                                                                                                                                       cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtctc                      <  2:2009683/143‑1 (MQ=255)
                                                                                                                                       aacgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtat               >  2:139311/6‑150 (MQ=255)
                                                                                                                                        acgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata              >  1:2280622/5‑150 (MQ=255)
                                                                                                                                         aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTatat             <  2:373060/150‑1 (MQ=255)
                                                                                                                                          ttAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATc            <  2:1246962/150‑1 (MQ=255)
                                                                                                                                          ttAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATc            <  1:893634/150‑1 (MQ=255)
                                                                                                                                           tAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATcc           <  2:2228808/150‑1 (MQ=255)
                                                                                                                                            aaaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCa          <  1:2162966/150‑1 (MQ=255)
                                                                                                                                               cTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATAc       <  2:2042866/150‑1 (MQ=255)
                                                                                                                                                   gcTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCAGCCATTTAAGGCGGGAAACCTCATCCTATTGGAGAAAAATGAATAATTAAATTGCTTGGTTATCATCTCCCAGTATATCCATTCTaa  <  2:131891/150‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACA  >  NC_000913/1223780‑1224077

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