Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,062 T→C 5.1% 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 5.1% 667.4 / 16.7 196intergenic (+14/‑149)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (86/100);  new base C (4/6);  total (90/106)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 202 total aligned reads displayed.

GATCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCA  >  NC_000913/1809913‑1810211
                                                                                                                                                     |                                                                                                                                                     
  tCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGa                                                                                                                                                     <  1:4047724/150‑1 (MQ=255)
    gCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                                                                   >  2:3708007/1‑150 (MQ=255)
    gCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGGTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                                                                   <  2:650886/150‑1 (MQ=255)
       tctAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCg                                                                                                                                                >  2:4077274/1‑150 (MQ=255)
       tctAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCAGCGCCAGAAGAGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACATTAATAGCTGACAGAACTCACCGCCAAAGACCTTCTCGGTTAATTACCAGGGGCAGTGATCTCg                                                                                                                                                >  1:1916863/1‑150 (MQ=255)
         tAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCt                                                                                                                                              <  2:3399360/150‑1 (MQ=255)
           aaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGc                                                                                                                                            >  2:856089/1‑150 (MQ=255)
             gcagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                                                                          <  1:3357601/150‑1 (MQ=255)
             gcagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGccc                                                                                                                                          <  2:1225155/150‑1 (MQ=255)
              cagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                         <  2:1272384/150‑1 (MQ=255)
               agcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                                                                        <  1:626670/150‑1 (MQ=255)
                gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAACATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACGCACCGCAAACGACCTTCTCGGTTAAGGACCAGGGGCAGTGATCTCGCTGCCCCTg                                                                                                                                       <  2:223172/150‑1 (MQ=255)
                gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTg                                                                                                                                       <  2:1154829/150‑1 (MQ=255)
                gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                                                                       <  1:3102955/150‑1 (MQ=255)
                gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                                                                       <  1:3777838/150‑1 (MQ=255)
                gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                                                                       <  2:551288/150‑1 (MQ=255)
                 cccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                                                                      <  2:2115878/150‑1 (MQ=255)
                 cccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                                                                      <  2:3725827/150‑1 (MQ=255)
                  ccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCAACTTTTTTATTATCAGACGTCAAAATTTAAACCCTGACAGAACACACCGCAAAAGAACATCTCGGTTAATGACCAGGGGCAGTGATCTAGCGGCCCCAGGt                                                                                                                                     >  1:2944324/1‑150 (MQ=255)
                      aTGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                                                                 <  2:1122691/150‑1 (MQ=255)
                           ttCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCt                                                                                                                            <  1:1537419/150‑1 (MQ=255)
                           ttCCATCGTCGTTCCAGCGCCAGAAGCGCAAAAATATACACCTTTTTTATTAGCAGAAGTCAAAAATTCAACGCTGACAGAACACACCCCAAAAGACCATCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCAGGttatttatat                                                                                                                            >  1:2382596/1‑148 (MQ=255)
                                 cgtcgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                                                                      <  2:472257/150‑1 (MQ=255)
                                 cgtcggTCCTGCGCCAGAAGCGCAAAATGCTCCCCGTTTTGTATTAGCCGACGTCAAACTGTCATCGCGGACAGAACTCCCCGCAACAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                                                                      <  2:2548926/150‑1 (MQ=255)
                                   tcgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                                                                    <  1:2608861/150‑1 (MQ=255)
                                    cgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                                                                                                   <  2:3677294/150‑1 (MQ=255)
                                     gtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                                                                  <  2:2143255/150‑1 (MQ=255)
                                          tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                                                                      >  2:3735573/1‑141 (MQ=255)
                                            cgcCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTg                                                                                                           >  2:1611609/1‑150 (MQ=255)
                                              ccAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAGTCAACTGAc                                                                                                         >  2:3218188/1‑150 (MQ=255)
                                               cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCCTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        <  2:543134/150‑1 (MQ=255)
                                               cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        >  1:3270974/1‑150 (MQ=255)
                                                 gAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                                                                                      >  1:3989799/1‑150 (MQ=255)
                                                       cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                                >  2:1576287/1‑150 (MQ=255)
                                                           aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                                                                            <  1:3313391/150‑1 (MQ=255)
                                                            tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                                                                                           <  2:3665602/150‑1 (MQ=255)
                                                              aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                                                                         >  2:1678692/1‑150 (MQ=255)
                                                               tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                                                                             <  2:76899/145‑1 (MQ=255)
                                                               tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                                                                             >  1:76899/1‑145 (MQ=255)
                                                                 cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                                                                                      <  2:2431916/150‑1 (MQ=255)
                                                                   cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                                                                    <  1:61215/150‑1 (MQ=255)
                                                                     ttttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga                                                                                  >  2:312130/1‑150 (MQ=255)
                                                                       ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa                                                                                <  2:1506589/150‑1 (MQ=255)
                                                                       ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa                                                                                <  1:3021087/150‑1 (MQ=255)
                                                                       ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa                                                                                <  2:397681/150‑1 (MQ=255)
                                                                         gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                              <  1:1869343/150‑1 (MQ=255)
                                                                                cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                                                                       >  1:3854189/1‑150 (MQ=255)
                                                                                cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                                                                       >  2:3244164/1‑150 (MQ=255)
                                                                                 aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                                                                      >  1:2742452/1‑150 (MQ=255)
                                                                                  gACGTCAAACTTTCCTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCg                                                                     <  2:28626/150‑1 (MQ=255)
                                                                                      tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                                                                 >  2:1306269/1‑150 (MQ=255)
                                                                                      tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                                                                 >  2:3450517/1‑150 (MQ=255)
                                                                                      acAAACTTACATCGCTGACAGAACTCACCGCAAAACACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCATTAACTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                                                                 <  2:1075860/149‑1 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  1:862302/1‑150 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  2:2366195/1‑150 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  1:255861/1‑150 (MQ=255)
                                                                                             ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTcc                                                          <  2:2746174/150‑1 (MQ=255)
                                                                                               cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                        <  1:492757/150‑1 (MQ=255)
                                                                                                aTCGCTGAGAGAACTCACCGCAAAAGACCTTCTCGGTTTATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGTAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTg                                                       <  2:263575/150‑1 (MQ=255)
                                                                                                aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTg                                                       >  2:3901678/1‑150 (MQ=255)
                                                                                                 tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGc                                                      >  1:68547/1‑150 (MQ=255)
                                                                                                  cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa                                                     <  1:3326234/150‑1 (MQ=255)
                                                                                                  cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa                                                     <  1:1899864/150‑1 (MQ=255)
                                                                                                   gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAc                                                    <  1:2715290/150‑1 (MQ=255)
                                                                                                   gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGAATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAc                                                    >  1:2254822/1‑150 (MQ=255)
                                                                                                          gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaa                                             <  1:2178373/150‑1 (MQ=255)
                                                                                                          gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaa                                             >  2:3848941/1‑150 (MQ=255)
                                                                                                          gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaa                                             <  1:3645285/150‑1 (MQ=255)
                                                                                                          gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaa                                             >  2:7063/1‑150 (MQ=255)
                                                                                                                aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGTGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATATAAATcc                                       >  2:1554969/1‑150 (MQ=255)
                                                                                                                 ccGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGGTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCt                                      >  1:3485709/1‑150 (MQ=255)
                                                                                                                 ccGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCt                                      >  1:4097228/1‑150 (MQ=255)
                                                                                                                  cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCta                                     <  2:2042355/150‑1 (MQ=255)
                                                                                                                    cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtata                                   <  1:2012063/150‑1 (MQ=255)
                                                                                                                         gACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                              >  2:1065131/1‑150 (MQ=255)
                                                                                                                         gACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                              <  2:1787947/150‑1 (MQ=255)
                                                                                                                          aCCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGACACCACTTCCGTGCACTGTATAAAAATCCTATACTGTAc                             >  1:3689866/1‑150 (MQ=255)
                                                                                                                             ttCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTa                          >  2:1550238/1‑150 (MQ=255)
                                                                                                                              tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGATATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTAt                         <  1:2627121/150‑1 (MQ=255)
                                                                                                                              tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTATGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTAt                         <  1:947506/150‑1 (MQ=255)
                                                                                                                              tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTAt                         >  1:2974433/1‑150 (MQ=255)
                                                                                                                                tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCg                       >  2:3497925/1‑150 (MQ=255)
                                                                                                                                tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCg                       >  2:1351524/1‑150 (MQ=255)
                                                                                                                                tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCg                       >  2:2105903/1‑150 (MQ=255)
                                                                                                                                   gTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACa                    <  1:3708007/150‑1 (MQ=255)
                                                                                                                                   gTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACa                    >  1:634741/1‑150 (MQ=255)
                                                                                                                                    taaagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAg                   >  2:2949074/7‑150 (MQ=39)
                                                                                                                                          tacAGGGGGAGTGATATCGGAGGCCCTTGTTATTTTTCTGAAATTGGATTCAAACTACGGATTAATCGTAAATTTTAGAGTAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGt             <  2:3270974/148‑1 (MQ=255)
                                                                                                                                            cAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGa           <  2:194782/150‑1 (MQ=255)
                                                                                                                                             aGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACATCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAg          >  1:3115958/1‑150 (MQ=255)
                                                                                                                                               gggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtt        <  1:1983783/150‑1 (MQ=255)
                                                                                                                                               gggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtt        <  1:2699753/150‑1 (MQ=255)
                                                                                                                                                 gCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtttt      >  1:352927/1‑150 (MQ=255)
                                                                                                                                                 gCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtttt      >  1:332039/1‑150 (MQ=255)
                                                                                                                                                  ccgTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATGAATCGTCAATTTAAGCGAAAGAGTTACACCGTCACCCCTTCCGTGCACTGTATAGAAATCCTATACTGTACGTCTCGACAGGGTCGTGAGttttc     >  1:701216/3‑149 (MQ=255)
                                                                                                                                                  cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAGAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTa     >  1:1627054/1‑150 (MQ=255)
                                                                                                                                                  cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTa     <  1:1070928/150‑1 (MQ=255)
                                                                                                                                                  cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTa     <  2:895526/150‑1 (MQ=255)
                                                                                                                                                     tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAATATAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  >  1:2911945/1‑150 (MQ=255)
                                                                                                                                                     tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  >  2:1029872/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCA  >  NC_000913/1809913‑1810211

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