Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,062 T→C 8.6% intergenic (+14/‑149) yniC → / → ydjM 2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,810,0620TC 8.6% 293.1 / 13.9 107intergenic (+14/‑149)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (44/52);  new base C (4/5);  total (50/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 110 total aligned reads displayed.

CGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCA  >  NC_000913/1809916‑1810211
                                                                                                                                                  |                                                                                                                                                     
cGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAt                                                                                                                                                    >  2:540989/1‑150 (MQ=255)
  ccTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCCCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtct                                                                                                                                                  >  1:1438950/1‑150 (MQ=255)
  ccTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGATAATGACCAGGGGCAGTGAtca                                                                                                                                                  >  2:413775/1‑149 (MQ=255)
      tAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCt                                                                                                                                              <  1:462875/150‑1 (MQ=255)
           ccgcCCGCATGCGTTCCAGCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                         <  2:302653/148‑1 (MQ=255)
           cagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                                                                         <  1:1155492/150‑1 (MQ=255)
            agcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCGGccccg                                                                                                                                        >  1:620975/1‑149 (MQ=255)
            agcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCt                                                                                                                                        <  1:1014355/150‑1 (MQ=255)
            agcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCt                                                                                                                                        <  1:1143990/150‑1 (MQ=255)
            agcCCGCATGCGTTCCATCGTCGGTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTACTGCCCAGGGGCAGTGATCTCGCTGCCCCt                                                                                                                                        <  1:299304/150‑1 (MQ=255)
            agcCCGCAGGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGAAGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGGTAATGACCAGCGGCAGCGAGATCACTGCCCCt                                                                                                                                        <  1:271608/150‑1 (MQ=255)
             gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTg                                                                                                                                       >  1:625876/1‑150 (MQ=255)
              cccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                                                                      >  1:653400/1‑150 (MQ=255)
               ccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTGTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGt                                                                                                                                     <  1:531577/150‑1 (MQ=255)
                 gCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTc                                                                                                                                   >  2:601402/1‑150 (MQ=255)
                   aTGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGtcatt                                                                                                                                 <  2:775593/150‑5 (MQ=255)
                    tGCGTTCCATCGTCGGTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACCTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                                                                <  2:1088179/150‑1 (MQ=255)
                          ccATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                                                                          >  2:190942/1‑150 (MQ=255)
                           cATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaa                                                                                                                         >  1:464194/1‑150 (MQ=255)
                             tcgtcgTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                                                                       <  1:1258430/150‑1 (MQ=255)
                                  gtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGACCTCACCACAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                                                                  <  1:275649/150‑1 (MQ=255)
                                       tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                                                                             <  1:1287041/150‑1 (MQ=255)
                                           ccAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCCACTGAc                                                                                                         >  1:98736/1‑150 (MQ=255)
                                            cataatcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCTCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        <  1:866843/144‑1 (MQ=255)
                                            cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        <  1:559217/150‑1 (MQ=255)
                                              gAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                                                                                      >  2:1573162/1‑150 (MQ=255)
                                              gAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                                                                                      <  2:1413158/150‑1 (MQ=255)
                                               aaGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGTTTAATGACCAGGGGCAGTGATCTCGCTTCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                                                                     >  1:97347/1‑150 (MQ=255)
                                                aGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                                                                                                    >  1:326687/1‑150 (MQ=255)
                                                   gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                                                                                                 <  2:1161730/150‑1 (MQ=255)
                                                    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                                <  2:100732/150‑1 (MQ=255)
                                                       aaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                                                                             <  2:218504/150‑1 (MQ=255)
                                                          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                                                                          >  2:1620195/1‑150 (MQ=255)
                                                          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGATTAATccgtcagt                                                                                         <  2:1556351/150‑8 (MQ=255)
                                                           aaccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                                                                         <  2:714259/148‑1 (MQ=255)
                                                           aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACATTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCACTGGATCTTTATCTGAAATTGAATTCAACTGACGGATTAATCGTCAAtt                                                                                         >  1:336524/1‑150 (MQ=255)
                                                            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                                                                        >  1:779229/1‑150 (MQ=255)
                                                            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                                                                        >  1:772685/1‑150 (MQ=255)
                                                            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                                                                        <  2:776521/150‑1 (MQ=255)
                                                             ccACGTTTTGTATTAGCAGACGTCAAACTTTCTTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                                                                       <  2:376810/150‑1 (MQ=255)
                                                              cACGTTTTGTATTAGCAGACGTCAAACTTTCAACGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                                                                                      >  1:47224/1‑150 (MQ=255)
                                                               aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                                                                     <  1:1175220/150‑1 (MQ=255)
                                                                      gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCCCCGCATAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                              >  1:1420401/1‑150 (MQ=255)
                                                                      gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                              <  2:1378189/150‑1 (MQ=255)
                                                                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                                                                             <  1:1663720/150‑1 (MQ=255)
                                                                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                                                                             <  1:1573162/150‑1 (MQ=255)
                                                                          tAGCAGACGTCAAACTTTCATCGCTGACAGAACTAACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtgt                                                                          <  2:653400/150‑3 (MQ=255)
                                                                          tAGAAGACGTCAAACTTTCATCGCTGACTGAACTCACAGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                                                                          <  1:1065513/150‑1 (MQ=255)
                                                                                  gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACCCCGTcc                                                                  >  2:571163/1‑149 (MQ=255)
                                                                                   tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                                                                 <  2:569797/150‑1 (MQ=255)
                                                                                    cAAACTTTCATCGCTGACAGAACTCCCCGCAAAAGACCATCTAGGTTAATGACCAGGGGCAGAGATCTCGCTGCCCCTGGTAAAATAGCTGAAAGTGCATACAACTGCCGGATTAAGCGTCAACTTAAGAGATAGACTTCCACCGTcaca                                                                >  1:286755/1‑149 (MQ=255)
                                                                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc                                                                <  1:975667/150‑1 (MQ=255)
                                                                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  2:237933/1‑150 (MQ=255)
                                                                                         tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc                                                           <  1:1352153/150‑1 (MQ=255)
                                                                                            cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                        >  2:770380/1‑150 (MQ=255)
                                                                                            cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                        >  2:1544169/1‑150 (MQ=255)
                                                                                               cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa                                                     >  1:1651976/1‑150 (MQ=255)
                                                                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCCCCACTTCCGTGCAc                                                    >  1:873700/1‑150 (MQ=255)
                                                                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAc                                                    >  1:1731737/1‑150 (MQ=255)
                                                                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACt                                                   <  1:532597/150‑1 (MQ=255)
                                                                                                   gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGAGATCACGCTGCCCCTCGTTCTTTATCAGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGCGTTACACCGTCACCACTTCCGTCCACTGt                                                 >  2:1575987/1‑150 (MQ=255)
                                                                                                     cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGtat                                               <  2:98736/150‑1 (MQ=255)
                                                                                                      aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACAGtata                                              >  1:218991/1‑150 (MQ=255)
                                                                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAGCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGCAAGAGTTACACCGTCACCACTTCCGGGCCCTGTAGaa                                             >  2:284220/1‑150 (MQ=255)
                                                                                                        aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                                                                          <  1:1237832/104‑1 (MQ=255)
                                                                                                        aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                                                                          <  1:1651029/104‑1 (MQ=255)
                                                                                                        aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                                                                          >  2:1651029/1‑104 (MQ=255)
                                                                                                        aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                                                                          >  2:1237832/1‑104 (MQ=255)
                                                                                                          cTCACCGCAAAAGACCTTCTCCGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTTAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaaaaa                                          >  2:131590/1‑150 (MQ=255)
                                                                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGGTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCATTTTAAGAGAAAGAGTTACACAGTCCCCACTTACGTGGACTGTATAAAAATc                                        >  2:1556603/1‑150 (MQ=255)
                                                                                                             cccGCAAAAGCCCTTCTCGGTTAATGACCACGGGCCGTGCTCTCGCTGCCCCTGGTTCTTTAACTGAACTTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACCCCGTCCCCACTTCCGTGCACTGTATAAAAATcc                                       <  2:711127/149‑1 (MQ=255)
                                                                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGTCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCTTCACCACTTCCGTGCACTGTATAAAAATcc                                       <  1:56175/150‑1 (MQ=255)
                                                                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGTATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATcc                                       <  1:565471/150‑1 (MQ=255)
                                                                                                             aCCGCAAAAGACCTTCTCGGTTAAAGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGGGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATcc                                       <  2:1714395/150‑1 (MQ=255)
                                                                                                               cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTTCACCGTCACCACTTCCGTGCACTGTATAAAAATCCta                                     >  2:727792/1‑150 (MQ=255)
                                                                                                                gCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtat                                    <  2:464194/150‑1 (MQ=255)
                                                                                                                  aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATAc                                  <  2:934061/150‑1 (MQ=255)
                                                                                                                      gACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                              >  1:965780/1‑150 (MQ=255)
                                                                                                                      gACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                              <  1:1620195/150‑1 (MQ=255)
                                                                                                                       aCCTTCTGGGTTAATGACCAGGGGCAGTGATATCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGTTTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTAc                             <  2:670563/150‑1 (MQ=255)
                                                                                                                        ccTTCTCGGTTAATGACCAGGGGCAGTTATCTCGCTGCCCCTGGTTCATTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACg                            <  1:1560870/150‑1 (MQ=255)
                                                                                                                        ccTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACg                            >  2:1473626/1‑150 (MQ=255)
                                                                                                                          ttCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACGTCCGTGCACTGTATAATAATCCTATACTGTACGTa                          >  1:481539/1‑150 (MQ=255)
                                                                                                                            ctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTTCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATc                        <  1:533774/150‑1 (MQ=255)
                                                                                                                             tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCg                       >  1:129314/1‑150 (MQ=255)
                                                                                                                                  tAATGCCCAGGGGCCGTGATCTCGCGGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGt                  <  2:265001/150‑1 (MQ=255)
                                                                                                                                     tGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTa               >  1:347796/1‑150 (MQ=255)
                                                                                                                                      gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAg              >  1:1521628/1‑150 (MQ=255)
                                                                                                                                      gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAg              >  1:619180/1‑150 (MQ=255)
                                                                                                                                      gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAg              <  2:1588577/150‑1 (MQ=255)
                                                                                                                                       aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGt             >  1:641877/1‑150 (MQ=255)
                                                                                                                                          aGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAg          >  1:402579/1‑150 (MQ=255)
                                                                                                                                          aGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAg          >  2:622128/1‑150 (MQ=255)
                                                                                                                                           ggggCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGt         >  2:1634296/1‑150 (MQ=255)
                                                                                                                                             ggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGttt       <  2:1438950/150‑1 (MQ=255)
                                                                                                                                             ggCAGTGATATCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAATATTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGttt       >  1:798289/1‑150 (MQ=255)
                                                                                                                                                 gTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATc   >  2:604853/1‑150 (MQ=255)
                                                                                                                                                  tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGACAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  <  2:1077570/150‑1 (MQ=255)
                                                                                                                                                  tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  <  1:413775/150‑1 (MQ=255)
                                                                                                                                                  tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  <  2:396328/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCA  >  NC_000913/1809916‑1810211

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: