Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,062 T→C 5.3% 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.3% 1340.7 / 41.1 395intergenic (+14/‑149)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (180/193);  new base C (8/13);  total (189/206)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 402 total aligned reads displayed.

GATCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCA  >  NC_000913/1809913‑1810211
                                                                                                                                                     |                                                                                                                                                     
gATCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGt                                                                                                                                                       <  1:2028155/150‑1 (MQ=255)
  tCGCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGa                                                                                                                                                     <  2:4137933/150‑1 (MQ=255)
      ctctAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtctc                                                                                                                                                 <  2:3374083/150‑1 (MQ=255)
       tctAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCg                                                                                                                                                >  2:2434798/1‑150 (MQ=255)
          aaaTCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAATCCCAAAAAGTTCCACGTTTTGTATTATCAGACGTCAAAATTTCAAAGCTGACAGAACTCACCGCAAAAGACCATATCGGTTAATGAACAGGGGCCGTGATATCGcag                                                                                                                                             >  1:6763896/1‑148 (MQ=255)
          aaaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATATCGCTg                                                                                                                                             >  1:190382/1‑150 (MQ=255)
          aaaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTAATCGCTGACAGAACTCACCGCAAAAGACCATCTCGGTTAATGACCAGGGGCAGTGATCTCGCTg                                                                                                                                             >  1:3399461/1‑150 (MQ=255)
          aaaGCAGCCCGCATGCGTTCCATCGTCGTTCCAGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCGTCTCGGGTAATGACCAGGGGCAGTGATCTAGCTg                                                                                                                                             >  1:2955788/1‑150 (MQ=255)
            agcagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcc                                                                                                                                           >  2:3507583/1‑150 (MQ=255)
             gcagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                                                                          >  2:672048/1‑150 (MQ=255)
             gcagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                                                                          <  1:3143731/150‑1 (MQ=255)
              cagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                         <  1:91181/150‑1 (MQ=255)
              cagcCCGCATGCGTTCCATCGTCGGTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                         <  1:1105228/150‑1 (MQ=255)
               agcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                                                                        >  2:6526671/1‑150 (MQ=255)
                 cccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTAATCGCTGACAGAACACACCGCAAAAGACCATCTCGGTTAATTACCAGGGGCAGTGATCTCGCTGCCCCAgg                                                                                                                                      >  1:5392199/1‑150 (MQ=255)
                  ccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGt                                                                                                                                     >  1:5534295/1‑150 (MQ=255)
                  ccGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGt                                                                                                                                     <  2:3863891/150‑1 (MQ=255)
                    gCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTc                                                                                                                                   >  2:5096591/1‑150 (MQ=255)
                       tGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                                                                <  2:2612455/150‑1 (MQ=255)
                         cGTTCCCTCGTCGTTCCTGCGCCAGAAGCGCAACATGATCCACGTTTTGTATTAGCAGACGTCAAACTGTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                                                                              <  2:2235397/150‑1 (MQ=255)
                         cGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                                                                              <  1:4178005/150‑1 (MQ=255)
                           ttCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCt                                                                                                                            <  1:4031068/150‑1 (MQ=255)
                           gtCCATCGGCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACGTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCt                                                                                                                            <  1:5815803/149‑1 (MQ=255)
                                 cgtcgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                                                                      <  1:4735787/150‑1 (MQ=255)
                                 cgtcgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                                                                      <  2:2986982/150‑1 (MQ=255)
                                  gtcgtTCCAGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTATCAGACGTCAAACTTTCAACGCTGACAGAACACACCCCAAAATACCTTATCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTATTTATATGAAAttt                                                                                                                     >  1:2808482/1‑149 (MQ=255)
                                    cgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                                                                                                   <  2:6624915/150‑1 (MQ=255)
                                         cTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                                                                                              >  1:4101253/1‑150 (MQ=255)
                                          tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                                                                             >  2:6426719/1‑150 (MQ=255)
                                          tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                                                                             >  2:2124844/1‑150 (MQ=255)
                                              ccAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                                                                         >  1:3039481/1‑150 (MQ=255)
                                               cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        <  1:6639647/150‑1 (MQ=255)
                                               cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        >  2:5633185/1‑150 (MQ=255)
                                               cAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                                                                        <  2:590656/150‑1 (MQ=255)
                                                aGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                                                                                       >  2:4433116/1‑150 (MQ=255)
                                                 gAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                                                                                      <  1:1667008/150‑1 (MQ=255)
                                                  aaGCGCAAACTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                                                                     <  2:3326937/150‑1 (MQ=255)
                                                    gcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                                                                                                   >  1:743799/1‑150 (MQ=255)
                                                    gcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                                                                                                   >  1:1268257/1‑150 (MQ=255)
                                                      gcAAAATGATCCGCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                                                                                                 >  1:5305342/1‑150 (MQ=255)
                                                      gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                                                                                                 <  2:2038231/150‑1 (MQ=255)
                                                       cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                                <  1:4452504/150‑1 (MQ=255)
                                                          aaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                                                                             >  1:13646/1‑150 (MQ=255)
                                                           aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                                                                            <  1:5633185/150‑1 (MQ=255)
                                                            tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                                                                                           >  1:2325818/1‑150 (MQ=255)
                                                              aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                                                                         >  1:3631460/1‑150 (MQ=255)
                                                               tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                                                                        >  2:781277/1‑150 (MQ=255)
                                                               tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                                                                        >  1:3336888/1‑150 (MQ=255)
                                                                  aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                                                                     <  2:4294255/150‑1 (MQ=255)
                                                                   cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                                                                    <  1:3670226/150‑1 (MQ=255)
                                                                   cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                                                                    <  2:5753619/150‑1 (MQ=255)
                                                                       ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa                                                                                <  1:5458320/150‑1 (MQ=255)
                                                                         gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                              >  2:5535897/1‑150 (MQ=255)
                                                                         gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                              >  2:3180178/1‑150 (MQ=255)
                                                                          tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                                                                             <  1:4550278/150‑1 (MQ=255)
                                                                             tAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                                                                          <  1:1866681/150‑1 (MQ=255)
                                                                              aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTac                                                                         >  2:4704506/1‑150 (MQ=255)
                                                                                   aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGt                                                                    >  1:5705187/1‑150 (MQ=255)
                                                                                     gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATc                                                                                                                             <  1:2410874/91‑1 (MQ=255)
                                                                                     gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                                                                  >  2:109634/1‑150 (MQ=255)
                                                                                      tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                                                                 <  1:5035458/150‑1 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               <  1:4336081/150‑1 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  1:4878524/1‑150 (MQ=255)
                                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                                                               >  1:2535382/1‑150 (MQ=255)
                                                                                               cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                        <  1:6199663/150‑1 (MQ=255)
                                                                                                  cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa                                                     <  1:6426719/150‑1 (MQ=255)
                                                                                                    cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                                                                             >  1:5242971/1‑124 (MQ=255)
                                                                                                    cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACt                                                   <  1:3180178/150‑1 (MQ=255)
                                                                                                       aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCGCTGta                                                <  2:1204766/150‑1 (MQ=255)
                                                                                                       aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGta                                                >  2:5988550/1‑150 (MQ=255)
                                                                                                       aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGta                                                >  1:6016456/1‑150 (MQ=255)
                                                                                                         aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGtata                                              >  2:1326075/1‑150 (MQ=255)
                                                                                                           aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACAGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaaa                                            >  2:1782211/1‑150 (MQ=255)
                                                                                                              tCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                                >  2:2726997/1‑95 (MQ=255)
                                                                                                              tCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAAt                                         >  2:5021359/1‑150 (MQ=255)
                                                                                                              tCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAAt                                         >  1:6118891/1‑150 (MQ=255)
                                                                                                               cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATc                                        <  2:368006/150‑1 (MQ=255)
                                                                                                                aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATcc                                       >  2:4830326/1‑150 (MQ=255)
                                                                                                                  cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCta                                     >  2:1489099/1‑150 (MQ=255)
                                                                                                                    cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtata                                   <  1:3910377/150‑1 (MQ=255)
                                                                                                                         gACCTTCTCGGTTAATGACCAGGCGCAGTGATCTCGCTGCCCCTGGAACATTATCTGAAATTGCATTCAACTGACGGATTAATCGACAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                              <  1:461667/150‑1 (MQ=255)
                                                                                                                          aCCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTAc                             <  2:3771243/150‑1 (MQ=255)
                                                                                                                            cTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGt                           <  1:327403/150‑1 (MQ=255)
                                                                                                                            cTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGt                           >  1:5988409/1‑150 (MQ=255)
                                                                                                                             ttCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTa                          >  2:2689805/1‑150 (MQ=255)
                                                                                                                                  ggTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGAc                     >  1:3771009/1‑150 (MQ=255)
                                                                                                                                      aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGt                               >  1:2716857/1‑136 (MQ=255)
                                                                                                                                      aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGtt                 <  2:5348921/150‑1 (MQ=255)
                                                                                                                                      aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGtt                 <  2:1635752/150‑1 (MQ=255)
                                                                                                                                      aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGtt                 >  1:2237469/1‑150 (MQ=255)
                                                                                                                                       aTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGttt                <  1:172201/150‑1 (MQ=255)
                                                                                                                                          aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGt             <  2:496998/150‑1 (MQ=255)
                                                                                                                                           ccAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTg            >  1:100010/1‑150 (MQ=255)
                                                                                                                                             aGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAg          >  1:1626041/1‑150 (MQ=255)
                                                                                                                                             aGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTa               <  2:4879489/145‑1 (MQ=255)
                                                                                                                                                 gCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtata                                   <  1:5226260/121‑1 (MQ=255)
                                                                                                                                                   aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGt                                                                                              >  2:5242473/1‑60 (MQ=255)
                                                                                                                                                   aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTAt    >  2:2507347/1‑150 (MQ=255)
                                                                                                                                                     tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  >  1:6495832/1‑150 (MQ=255)
                                                                                                                                                     tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCATGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATCa  >  2:6107350/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 < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: