Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,706 T→A 7.8% intergenic (+30/‑210) rho → / → wecA transcription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,7060TA 7.8% 416.2 / 11.3 144intergenic (+30/‑210)rho/wecAtranscription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base T (70/62);  new base A (5/7);  total (75/69)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

TCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCT  >  NC_000913/3967561‑3967855
                                                                                                                                                 |                                                                                                                                                     
tCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACg                                                                                                                                                   <  1:1489627/150‑1 (MQ=255)
  cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTg                                                                                                                                                 >  2:707201/1‑150 (MQ=255)
  cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTg                                                                                                                                                 >  1:927263/1‑150 (MQ=255)
   tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTgg                                                                                                                                                <  2:982195/150‑1 (MQ=255)
    gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGCCCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGc                                                                                                                                               >  1:448318/1‑150 (MQ=255)
     cgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCg                                                                                                                                              <  2:193564/150‑1 (MQ=255)
       cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGtt                                                                                                                                            >  2:401055/1‑150 (MQ=255)
       cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtt                                                                                                                                            <  1:696416/150‑1 (MQ=255)
          aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTg                                                                                                                                         <  1:908626/150‑1 (MQ=255)
          aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGttttt                                                                                                                                         <  2:466524/150‑2 (MQ=255)
           aTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGc                                                                                                                                        <  2:601783/150‑1 (MQ=255)
            tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtttttgg                                                                                                                                       <  2:358729/150‑4 (MQ=255)
              aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCttt                                                                                                                                     <  2:961214/150‑1 (MQ=255)
                 cACCCGTTGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGAACAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTtat                                                                                                                                  <  2:1254197/150‑1 (MQ=255)
                  aCCCGATGGGCGAAATCGATGCTATGGAATTCCTCATTAATAAACTGGAAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTtata                                                                                                                                 <  2:769523/150‑1 (MQ=255)
                         gggCGAAATCGATGCAATGGATTTCCTCATTAATAAACTGTAAATGACCAAGACCAATGACTATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCAACGTGTTTACGTGTCGTATTGCTTTTATATCTGTaa                                                                                                                          >  1:218705/1‑150 (MQ=255)
                            cGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCt                                                                                                                       >  2:1020059/1‑150 (MQ=255)
                               aaTCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTaa                                                                                                                    <  2:592451/150‑1 (MQ=255)
                                 tCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATg                                                                                                                  <  2:523438/150‑1 (MQ=255)
                                  cGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGc                                                                                                                 >  1:147444/1‑150 (MQ=255)
                                  cGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGc                                                                                                                 >  2:861951/1‑150 (MQ=255)
                                   gATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGcc                                                                                                                >  2:102795/1‑150 (MQ=255)
                                   gATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGcc                                                                                                                >  1:1253337/1‑150 (MQ=255)
                                   gATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCAAAAAACGCCACGTGTATACGTGGCGTGATGCTTTAATATCTGAAATCTTAATTcc                                                                                                                >  2:97621/1‑150 (MQ=255)
                                    aTGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCcg                                                                                                               <  1:1021451/150‑1 (MQ=255)
                                       cAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtttt                                                                                                                                          <  1:641641/120‑1 (MQ=255)
                                         aTGGAATTCCTCATTAATTAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTg                                                                                                          <  1:863381/150‑1 (MQ=255)
                                          tGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTgg                                                                                                         <  2:1470731/150‑1 (MQ=255)
                                            gAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCg                                                                                                       <  1:357664/150‑1 (MQ=255)
                                            gAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCg                                                                                                       <  1:526938/150‑1 (MQ=255)
                                             aaTTCCTCATCAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAATACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGa                                                                                                      <  1:1535976/150‑1 (MQ=255)
                                                tCCTCATTAATAAACTTGCAATGACCAAGACCAATGACGATTTCTTAGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGt                                                                                                   <  2:1092786/150‑1 (MQ=255)
                                                tCCTCATTAATAAACTGGCAATGTCCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGt                                                                                                   <  2:1098795/150‑1 (MQ=255)
                                                 ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGtt                                                                                                  >  1:629182/1‑150 (MQ=255)
                                                   tCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAGAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAg                                                                                                <  1:1277378/150‑1 (MQ=255)
                                                   tCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAg                                                                                                >  2:1425193/1‑150 (MQ=255)
                                                    cATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAgg                                                                                               >  1:696145/1‑150 (MQ=255)
                                                        aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGaaaa                                                                                           >  2:1250383/1‑150 (MQ=255)
                                                        aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAAGTGTTAACGTGGCGTTTTGCTTTTATATCTGTAATATTAATGCCGAGCTGGCGATGTTAGGaaaa                                                                                           >  1:244567/1‑150 (MQ=255)
                                                        aataaACTGGCAATGACCAAGACCAATGACGATTACTTCGAAATGAAGAAACGCTCCAAAATTTGTCTTATGCCAAAAACGGCCCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGaaaa                                                                                           <  2:279914/150‑1 (MQ=255)
                                                            aaCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTcc                                                                                       >  2:1486036/1‑150 (MQ=255)
                                                             aCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCt                                                                                      <  2:520771/150‑1 (MQ=255)
                                                             aCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCt                                                                                      <  2:1225611/150‑1 (MQ=255)
                                                             aCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCt                                                                                      <  1:1425193/150‑1 (MQ=255)
                                                              cTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTg                                                                                     >  1:362902/1‑150 (MQ=255)
                                                              cTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTg                                                                                     <  1:1020059/150‑1 (MQ=255)
                                                                ggCAATGACCAAGACCAATGACGATTTCTTCGAAAAGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGa                                                                                   <  2:257954/150‑1 (MQ=255)
                                                                  cAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAAt                                                                                 >  2:794215/1‑150 (MQ=255)
                                                                       accaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCt                                                                            <  1:401055/150‑1 (MQ=255)
                                                                        ccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTg                                                                           >  1:1378839/1‑150 (MQ=255)
                                                                         caagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTgg                                                                          <  2:1123183/150‑1 (MQ=255)
                                                                         caagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTgg                                                                          <  2:147444/150‑1 (MQ=255)
                                                                         caagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTgg                                                                          <  1:508925/150‑1 (MQ=255)
                                                                             accaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATg                                                                      <  2:1507418/150‑1 (MQ=255)
                                                                             accaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATg                                                                      >  2:1590844/1‑150 (MQ=255)
                                                                              ccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGt                                                                     >  1:182937/1‑150 (MQ=255)
                                                                              ccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGt                                                                     >  2:1029937/1‑150 (MQ=255)
                                                                                aaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTATTATTCCAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTAt                                                                  <  1:989327/150‑1 (MQ=255)
                                                                                 aTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAgg                                                                                               >  1:1457266/1‑121 (MQ=255)
                                                                                   gACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGc                                                                <  1:154402/150‑1 (MQ=255)
                                                                                    aCGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCa                                                               >  1:1110617/1‑150 (MQ=255)
                                                                                     cGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCaa                                                              >  1:948552/1‑150 (MQ=255)
                                                                                       aTTTCTTCGAAATGATGAATCGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAAAGCCGCGCTGGCGATGTTAGGAAATTTCCTGTAATTTGCTGTCATGTTATGCAAtt                                                            >  1:180937/1‑150 (MQ=255)
                                                                                          tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                         <  2:1253337/150‑1 (MQ=255)
                                                                                          tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                         <  2:1150226/150‑1 (MQ=255)
                                                                                          tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                         <  2:1427295/150‑1 (MQ=255)
                                                                                            ttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCat                                                       <  1:31009/150‑1 (MQ=255)
                                                                                               gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATc                                                    >  2:1224240/1‑150 (MQ=255)
                                                                                               gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATc                                                    >  1:1112044/1‑150 (MQ=255)
                                                                                                 aaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAATTGCATAACat                                                  <  2:948552/150‑2 (MQ=255)
                                                                                                   tgatgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAAt                                                >  1:1025389/1‑150 (MQ=255)
                                                                                                   tgatgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAAt                                                >  1:121558/1‑150 (MQ=255)
                                                                                                    gatgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATg                                               >  1:1512547/1‑150 (MQ=255)
                                                                                                     atgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATgg                                              <  2:383110/150‑1 (MQ=255)
                                                                                                        aaaCGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAttt                                                           <  1:435790/134‑1 (MQ=255)
                                                                                                        aaaCGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAttt                                                           >  2:435790/1‑134 (MQ=255)
                                                                                                          aCGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAt                                         >  2:1034851/1‑150 (MQ=255)
                                                                                                           cGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtt                                        >  2:1564105/1‑150 (MQ=255)
                                                                                                             cTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtttt                                      >  2:426644/1‑150 (MQ=255)
                                                                                                               cATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTATTGCTTTAATATCTGTAATCTTAATGACGCGCTGGCGATGTTATGAAAATTCCTGTAATTTGCTGTCATGTTATGAAATTTTCATATCAAATGGTTAATTTTTg                                    >  1:1267052/1‑150 (MQ=255)
                                                                                                                 tAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGca                                  <  2:696145/150‑1 (MQ=255)
                                                                                                                   aaTTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGcaca                                >  1:759312/1‑150 (MQ=255)
                                                                                                                     tttGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAgg                              <  1:1188107/150‑1 (MQ=255)
                                                                                                                     tttGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAgg                              <  2:121558/150‑1 (MQ=255)
                                                                                                                     tttGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAgg                              >  1:1244037/1‑150 (MQ=255)
                                                                                                                               gCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                    <  1:812478/150‑1 (MQ=255)
                                                                                                                               gCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                    <  2:927263/150‑1 (MQ=255)
                                                                                                                               gCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                    >  1:404592/1‑150 (MQ=255)
                                                                                                                                  aaaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTg                 <  2:448318/150‑1 (MQ=255)
                                                                                                                                  aaaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCCTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTg                 <  1:710530/150‑1 (MQ=255)
                                                                                                                                   aaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTgg                <  1:1564105/150‑1 (MQ=255)
                                                                                                                                    aaaCGCCACGTAAACACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGa               >  1:663576/1‑150 (MQ=255)
                                                                                                                                    aaaCGCCACGTAAACACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAGTTTTGCAGAGGACTGGTGGGTTTggc               >  2:284433/1‑149 (MQ=255)
                                                                                                                                      aCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAAc             >  2:366181/1‑150 (MQ=255)
                                                                                                                                      aCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAAc             >  2:1592669/1‑150 (MQ=255)
                                                                                                                                        accACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACgg           >  1:931249/2‑150 (MQ=255)
                                                                                                                                         ccACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGa          >  1:1670935/1‑150 (MQ=255)
                                                                                                                                           aCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGaa              >  1:1175902/1‑144 (MQ=255)
                                                                                                                                           aCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGaa              <  2:1175902/144‑1 (MQ=255)
                                                                                                                                            cGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACtt       <  2:558685/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCT  >  NC_000913/3967561‑3967855

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