Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,352 C→T 5.0% R151W (CGG→TGG)  ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,3520CT 5.0% 391.8 / 7.3 127R151W (CGG→TGG) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base C (69/46);  new base T (3/3);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

Only 100 of 146 total aligned reads displayed.

GAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATT  >  NC_000913/569204‑569501
                                                                                                                                                    |                                                                                                                                                     
gAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCg                                                                                                                                                      >  1:1183894/1‑150 (MQ=255)
gAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCg                                                                                                                                                      >  1:2288546/1‑150 (MQ=255)
gAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCg                                                                                                                                                      <  2:2594333/150‑1 (MQ=255)
   aTTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGt                                                                                                                                                   >  1:1790742/1‑150 (MQ=255)
    ttGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGtt                                                                                                                                                  >  1:469794/1‑150 (MQ=255)
         tttGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAAtt                                                                                                                                             >  2:1480100/1‑150 (MQ=255)
         tttGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAAtt                                                                                                                                             >  1:2681398/1‑150 (MQ=255)
            gAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATc                                                                                                                                          >  1:1775462/1‑150 (MQ=255)
             aaTCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCa                                                                                                                                         <  1:561751/150‑1 (MQ=255)
             aaTCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCa                                                                                                                                         <  1:400071/150‑1 (MQ=255)
              aTCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCAt                                                                                                                                        >  1:1650265/1‑150 (MQ=255)
              aTCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCAt                                                                                                                                        >  1:1438687/1‑150 (MQ=255)
               tCACGGTATTTATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATg                                                                                                                                       <  2:2695238/150‑1 (MQ=255)
               tCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATg                                                                                                                                       <  2:1025177/150‑1 (MQ=255)
               tCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATg                                                                                                                                       <  2:2695205/150‑1 (MQ=255)
                 aCGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGga                                                                                                                                     <  1:165861/150‑1 (MQ=255)
                 aCGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGga                                                                                                                                     <  1:484545/150‑1 (MQ=255)
                  cGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTTGGGTTAAATTATCATGgaa                                                                                                                                    >  2:583698/1‑150 (MQ=255)
                  cGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATCAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTTGGGTTAAATTATCATGgaa                                                                                                                                    >  2:583722/1‑150 (MQ=255)
                         gATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGaaaaaa                                                                                                                             <  1:1557702/150‑1 (MQ=255)
                         gATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGaaaaaa                                                                                                                             <  2:1374249/150‑1 (MQ=255)
                          aTGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGTAAGAAAAAAc                                                                                                                            >  1:1633882/1‑150 (MQ=255)
                                 aaCTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATTTTCATTAATAAAAGCCATACTTATAGCACAACGGGCAAATGAAGAAAGCGAGATAACGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGa                                                                                                                     <  2:463328/150‑1 (MQ=255)
                                 aaCTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGa                                                                                                                     >  1:1546258/1‑150 (MQ=255)
                                 aaCTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATAATGGAAGAAAAAACGGCAGGa                                                                                                                     >  2:1015566/1‑150 (MQ=255)
                                    tctTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGc                                                                                                                  >  1:955003/1‑150 (MQ=255)
                                    tctTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGc                                                                                                                  >  1:533505/1‑150 (MQ=255)
                                    tctTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGc                                                                                                                  >  1:2299130/1‑150 (MQ=255)
                                      tttGCTACAATAAAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTNGGGNTAAATTATCATGGANGAAANAACGNCAGGATGCAc                                                                                                                <  2:2097169/150‑1 (MQ=255)
                                      tttGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCAc                                                                                                                <  2:533505/150‑1 (MQ=255)
                                           gACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGaa                                                                                                           >  1:680271/1‑150 (MQ=255)
                                              aaTACAGTCTATAATATTGACTCTTTGAATTAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAATAAAAAACGGCAGGATGCACTGGAATCa                                                                                                        <  2:1103607/150‑1 (MQ=255)
                                               aTACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAg                                                                                                       <  2:227072/150‑1 (MQ=255)
                                                          aataTTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATg                                                                                            >  2:1594102/1‑150 (MQ=255)
                                                            tttttaCTCTTTGAATGAGCCAAATTAATTAAAAAACCCCAAAAATAACCCAAAACGGCAAAATCAAGAAAGCGATATAAAGACATGTCGGTTAAATTTATCATGGAATAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGAc                                                                                          <  1:2041547/145‑1 (MQ=255)
                                                                 aCTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGt                                                                                     >  2:2181633/1‑150 (MQ=255)
                                                                  ctctTTGAATGAGCCATTTTAATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCTAGATAAAGTCAAGTCGGGTTAATTTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    <  2:2131738/150‑1 (MQ=255)
                                                                  ctctTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAATGGAAGAAAGCGAGATAACGCCAAGTCGGGTTAAATTATCTTGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    <  2:2131728/150‑1 (MQ=255)
                                                                  ctctTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    <  2:197541/150‑1 (MQ=255)
                                                                  ctctTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    >  2:2203083/1‑150 (MQ=255)
                                                                  ctctTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    >  1:1078378/1‑150 (MQ=255)
                                                                  ctctTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTc                                                                                    >  1:1078364/1‑150 (MQ=255)
                                                                     tttGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTg                                                                                 >  2:1437390/1‑150 (MQ=255)
                                                                      ttGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGt                                                                                >  2:951696/1‑150 (MQ=255)
                                                                      ttGAATGAGCAATATAAATAAATAAAACCCATACTTATAGCACAATGGTCAATTGAAGAAAGCGAGATAACGACAAGTGGGTTTAAATTATCATGGAACAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGt                                                                                <  1:2522434/150‑1 (MQ=255)
                                                                         aaTGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCg                                                                             >  1:2730894/1‑150 (MQ=255)
                                                                         aaTGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCg                                                                             >  1:2691018/1‑150 (MQ=255)
                                                                            gAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCgaga                                                                          >  2:1398437/1‑150 (MQ=255)
                                                                            cacaaatatTAATAAAAAAAACAAATACTTATAGCAAAATGGCAAAATCAATAAAGCTTTATAAATTAAAGTCGGTTTAATTTATAATGGAATAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCgaga                                                                          <  1:2438710/145‑1 (MQ=255)
                                                                             aGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGAt                                                                         >  1:2314550/1‑150 (MQ=255)
                                                                              gCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATg                                                                        >  2:2376917/1‑150 (MQ=255)
                                                                               ccATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATgg                                                                       >  1:978799/1‑150 (MQ=255)
                                                                                     taaTTAAAAAAAGCAATATATATACCACAATGGCAAATTGAATAAAGCGATATAATGTAATGTGGGTTTAATTTATCATGCAATAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTcc                                                                 <  1:2131658/148‑1 (MQ=255)
                                                                                      aaTTAATAAAAGCCATACTTATCGCACACAGGGAAAATTAAGAAAGCTAGATAAATTCATGTCGTGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCt                                                                <  1:2344329/149‑1 (MQ=255)
                                                                                      aaTTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCATGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCt                                                                <  2:1722890/149‑1 (MQ=255)
                                                                                      aaTTAATAAAAGCAATACTTATAGCACAACGGGCATATGAATAAAGCGAGATAAAGTCATGTCGGGTTAATTTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCt                                                                <  1:1842977/149‑1 (MQ=255)
                                                                                      aaTTAAAAAAACCAATACTTATAGCACAACGGCAAAATTAAGAACGCTATATAAAGTCAAGTGGGGTTAAATTATCATGGAACAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCt                                                                <  1:2344295/149‑1 (MQ=255)
                                                                                       aTTAATAAAAGCCATACTTATCGCACAATGGGCAAATGAAGAAAGCGAGATAAAGTCATGTCGGGTTATATTATCATGTAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCtt                                                               <  1:2457180/150‑1 (MQ=255)
                                                                                       aTTAATAAAAGCCATACTTATAGCACAACGGGCAACTGACGAACGCGAGATACAGTCACGGGGGGTTAAATTAGCCTGGCACAACCAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCtt                                                               <  2:955003/150‑1 (MQ=255)
                                                                                       aTTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCtt                                                               >  2:2174216/1‑150 (MQ=255)
                                                                                        ttAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTa                                                              <  2:1666476/150‑1 (MQ=255)
                                                                                        ttAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTa                                                              <  2:2288546/150‑1 (MQ=255)
                                                                                          aataaAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGa                                                            <  2:2282872/150‑1 (MQ=255)
                                                                                          aataaAAGAAATACTTATAGAACAAAGGTAAAATTAATAAAGCGATATAATGTAATTTTGTTTTAAATTATCTTTTAATAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGa                                                            <  1:583722/150‑1 (MQ=255)
                                                                                          aataaAACCATTACTTATAGAAAAAAGGGAAATTTAATAAAGCTATATAAAGTAATTTTGGGTTAAATTATCATGGAATAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGa                                                            <  1:583454/150‑1 (MQ=255)
                                                                                           ataaAAGCCATCCTCATAGCACACAGGGCAAATGAAGAACGCGAGATAACGTCATGTCGGGTTAACTTATCATGGTAGAAAATTCGGCAGGATGCACTGGATTCAGTTACGATTATGACGGCGTCTTGTGCGAGATGGCTCTCCTTAGAt                                                           <  2:2469955/150‑1 (MQ=255)
                                                                                           ataaAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGTAGAAAAATCGGCTGGTTGCACTGGAATCTGGTAGGATTATGTCGGCGTCTTGTGCGAGATGGCTCTCGtttgtt                                                           >  2:416879/1‑146 (MQ=255)
                                                                                              aaaGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGAc                                                        >  2:1210849/1‑150 (MQ=255)
                                                                                                    aTACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGa                                                  >  2:81826/1‑150 (MQ=255)
                                                                                                        ttATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACgg                                              >  1:402457/1‑150 (MQ=255)
                                                                                                        ttATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACgg                                              >  2:904596/1‑150 (MQ=255)
                                                                                                          ataGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCt                                            >  2:1764599/1‑150 (MQ=255)
                                                                                                          ataGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGACAAGACGGTTTAAATTATCAAGGAAGAAAAAAAGGCAGGATTCACTGCAATCAGGAACTATTATGACGGAGACTTTTCCGAGATGGCTCTCCTTAGATTACAAAATAACTGCt                                            >  1:2595927/1‑150 (MQ=255)
                                                                                                             gCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCtttt                                         <  2:915813/150‑1 (MQ=255)
                                                                                                              cacaAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTg                                        >  2:1794732/1‑150 (MQ=255)
                                                                                                              cacaAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTg                                        >  1:1124749/1‑150 (MQ=255)
                                                                                                              cacaAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGACTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTg                                        >  1:1504961/1‑150 (MQ=255)
                                                                                                                caAAGGGCAAATGAATAAAGCGAGATAACGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGtt                                      <  2:521689/150‑1 (MQ=255)
                                                                                                                caAAGGCCAAATGAAGAAAGCCTCATAAAGTAATTTGGGGTTAAATTATCATGTAACAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGtt                                      <  1:2403327/150‑1 (MQ=255)
                                                                                                                 aaaGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTc                                     >  2:2528046/1‑150 (MQ=255)
                                                                                                                  aaGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTcc                                    <  1:1594102/150‑1 (MQ=255)
                                                                                                                     ggCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGa                                 >  2:2649368/1‑150 (MQ=255)
                                                                                                                        aaaTGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGAccc                              <  2:2015767/150‑1 (MQ=255)
                                                                                                                          aTGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCg                            >  2:2277825/1‑150 (MQ=255)
                                                                                                                             aagaaAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACa                         >  1:1357334/1‑150 (MQ=255)
                                                                                                                              agaaAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAg                        >  2:1086627/1‑150 (MQ=255)
                                                                                                                              agaaAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAg                        <  2:435055/150‑1 (MQ=255)
                                                                                                                              agaaAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAg                        >  2:817145/1‑150 (MQ=255)
                                                                                                                                         tAAAGTCAAGTCGGGTTAAATTATCATGGAATAAAAAACGGCAGGATGCACTGGAATCAGGTACTATTATGACGGCGTCTTGTCCGAGATGGCTCTCCATAGATGACAAAATAAAGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTa             >  2:2265114/1‑150 (MQ=255)
                                                                                                                                          aaaGTCAAGTCGGGTTAAATTCTCAGGGAAGAAAACACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTAt            <  2:1546258/150‑1 (MQ=255)
                                                                                                                                          aaaGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTAt            <  2:1775462/150‑1 (MQ=255)
                                                                                                                                          aaaGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTAt            <  2:1576960/150‑1 (MQ=255)
                                                                                                                                           aagtcaagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGACTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTAtt           <  2:2299130/150‑1 (MQ=255)
                                                                                                                                           aagtcaagtcGGGTTAAATTAACATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTAtt           <  2:2156581/150‑1 (MQ=255)
                                                                                                                                             gtcaagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGAAAAAAGAACGGGTTTTGTTGCATACCCCGACAGGGTGAAAACTATTgt         >  1:2247516/1‑149 (MQ=255)
                                                                                                                                              tcaagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAg        <  2:1464132/150‑1 (MQ=255)
                                                                                                                                              tcaagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAg        <  2:128953/150‑1 (MQ=255)
                                                                                                                                              tcaagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAg        >  1:2529312/1‑150 (MQ=255)
                                                                                                                                               caagtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGc       <  2:1786989/150‑1 (MQ=255)
                                                                                                                                                 agtcGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTa     >  2:2087000/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGATGGCTCTCCTTAGATGACAAAAGAACGGCTTTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATT  >  NC_000913/569204‑569501

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: