Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,948 A→C 5.3% K465Q (AAG→CAG)  ynfF → S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,661,9480AC 5.3% 307.1 / 4.2 103K465Q (AAG→CAG) ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (51/46);  new base C (4/2);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATA  >  NC_000913/1661801‑1662095
                                                                                                                                                   |                                                                                                                                                   
cTGGGGGTAGAATGGTTCCCGATGCTCGAGAATACTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAg                                                                                                                                                   >  1:296095/1‑150 (MQ=255)
 tGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGt                                                                                                                                                  >  2:209987/1‑150 (MQ=255)
    gtgTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtttt                                                                                                                                               <  2:1079544/148‑1 (MQ=255)
    gggTAGAATGGTTCCCGATGCTCGAGATTCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtttt                                                                                                                                               <  1:1980154/150‑1 (MQ=255)
    gggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtttt                                                                                                                                               >  2:107069/1‑150 (MQ=255)
    gggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCGTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGAGACCCGTGATGGTGTTAGTGGATAAGATAAACTGGATGTCCCCATCAAGtttt                                                                                                                                               >  1:227419/1‑150 (MQ=255)
      gTAGAATGGGTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTa                                                                                                                                             <  2:492361/150‑1 (MQ=255)
         gAATGGTTCCCGATTCTCGAGAATCCTGTTAAAACGCATATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                          <  2:1340643/150‑1 (MQ=255)
         gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                          >  1:1780143/1‑150 (MQ=255)
         gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                          >  1:523846/1‑150 (MQ=255)
         gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                          >  1:1081906/1‑150 (MQ=255)
          aaTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGt                                                                                                                                         <  1:1993636/150‑1 (MQ=255)
          aaTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGt                                                                                                                                         <  2:850020/150‑1 (MQ=255)
           aTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTg                                                                                                                                        >  2:352643/1‑150 (MQ=255)
            ttgTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATTGTACGGAAATGTCCGCTACCCGTGATGGTTTTCGTGGAAATGAAAAACTGGATGTCCCCATCATTTTTTTATGGTGc                                                                                                                                       <  2:1368486/148‑1 (MQ=255)
            tGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTTGATGGGGACATCCAGtttttcttttcc                                                                                                                                       <  1:1617800/150‑8 (MQ=255)
            tGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGGCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTTGATGGGGACATCCAGtttttcttttcc                                                                                                                                       <  1:1076177/150‑8 (MQ=255)
                tCCCGATGCTCGAGAATCCTGATAAAACGCAGATTTCCTTCTTTACATGTACAGATGCCATCGACCATGGTACGGAAATGACCGCGACACGTGATGGTGTGCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACg                                                                                                                                   <  2:296095/150‑1 (MQ=255)
                   cGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCa                                                                                                                                >  1:665151/1‑150 (MQ=255)
                    gATGCTCGAGAATCCTGTTAAACCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAg                                                                                                                               <  2:756815/150‑1 (MQ=255)
                      tGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTa                                                                                                                             >  1:586030/1‑150 (MQ=255)
                      tGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTa                                                                                                                             >  1:876627/1‑150 (MQ=255)
                          cGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCTCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAacac                                                                                                                         >  1:54538/1‑150 (MQ=255)
                           gagaATCCTGTTAAAACGCAGATTTCCGTCTTTCCATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAAcaca                                                                                                                        <  1:1347742/150‑1 (MQ=255)
                           gagaATCCTGTTAAAACGCAGATTTCCGTCTTTACATCGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGGTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAAcaca                                                                                                                        <  2:347456/150‑1 (MQ=255)
                              aaTCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTg                                                                                                                     <  2:1912789/150‑1 (MQ=255)
                                   tGTTAAAACGCAGAATTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatcaa                                                                                                                <  1:1315231/150‑1 (MQ=255)
                                     ttAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGCCCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTATTATGGTGCTACGCCAGTAACACATTGatcaatc                                                                                                              >  2:1861770/1‑150 (MQ=255)
                                     ttAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatcaatc                                                                                                              <  1:887539/150‑1 (MQ=255)
                                           cGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATg                                                                                                        >  1:1860230/1‑150 (MQ=255)
                                           cGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATg                                                                                                        >  1:524564/1‑150 (MQ=255)
                                             cAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGc                                                                                                      >  1:1135422/1‑150 (MQ=255)
                                               gATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGa                                                                                                    >  2:1102575/1‑150 (MQ=255)
                                                aTTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGAc                                                                                                   >  1:1000581/1‑150 (MQ=255)
                                                   accGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGAGATGGTGTTCGTGGAAAAGAAAACCTGCATGTCCCCAACAAGTATATATGGAGCTACGCCAGTCACACATTGATCCATCAACATGGCGCCatc                                                                                                >  2:289552/2‑150 (MQ=255)
                                                     cGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaa                                                                                              >  1:1468100/1‑150 (MQ=255)
                                                       tCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatc                                                                                            <  1:1134030/150‑1 (MQ=255)
                                                       tCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatc                                                                                            <  1:641160/150‑1 (MQ=255)
                                                       tCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaata                                                                                            >  1:938807/1‑149 (MQ=255)
                                                       tCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCATGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatc                                                                                            >  1:1493368/1‑150 (MQ=255)
                                                        cTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatca                                                                                           <  1:2009767/150‑1 (MQ=255)
                                                        cTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatca                                                                                           >  1:1841931/1‑150 (MQ=255)
                                                         tttACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATcac                                                                                          >  1:1955594/1‑150 (MQ=255)
                                                            aCATGGACAGATGCCATCGACCATGGTACGGAAATGACCTCGACCCTTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACAcc                                                                                       <  1:1800133/150‑1 (MQ=255)
                                                             cATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGAGGAAAAGAAAAACTGGATGTCACCATCAAGTTATTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACAccc                                                                                      >  1:51980/1‑150 (MQ=255)
                                                               tGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTATTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCCCACCCAt                                                                                    >  2:1718501/1‑150 (MQ=255)
                                                                ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                                   >  2:1528599/1‑150 (MQ=255)
                                                                ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                                   >  2:1660471/1‑150 (MQ=255)
                                                                ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                                   <  2:1421904/150‑1 (MQ=255)
                                                                   cAGATGCCATCGACCATTGTACGGAAATGACCGCGACCCGTGATGGTGTGAGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTATTTATGGTGCTACGCCAGTATCACATTTATCACTCAACATGGCGACATCAATCACACCCATTAgg                                                                                <  2:965259/150‑1 (MQ=255)
                                                                   cAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAgg                                                                                <  2:1188313/150‑1 (MQ=255)
                                                                      aTGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACTCCAGTAACACATTCATCCATCACCATGGCGACATCAAACACACCCAGGCGGTGc                                                                             >  1:365696/1‑150 (MQ=255)
                                                                       tGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCt                                                                            <  1:107069/150‑1 (MQ=255)
                                                                               cccATGGTACTGAAATGACCGCTACCCGTGATGCGGTTCGTGGAAAAGAAAAACTGGATGTCGCCATCAAGTTTTTATGGTGCTACGCCAGTAAGACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                                    <  2:564104/149‑1 (MQ=255)
                                                                               aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGATGTGCTTCAGGATg                                                                    <  1:438527/150‑1 (MQ=255)
                                                                               aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                                    >  2:1010475/1‑150 (MQ=255)
                                                                                ccATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGa                                                                   >  2:850899/1‑150 (MQ=255)
                                                                                acATTGTACGGAAATGACCGCTACCCGTGATGGTTTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGGTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGAc                                                                  <  1:1234690/149‑1 (MQ=255)
                                                                                 cATGGTACGGAAATGACCGCGATCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGAc                                                                  <  2:1750118/150‑1 (MQ=255)
                                                                                        cGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCTAGt                                                           >  1:1079149/1‑150 (MQ=255)
                                                                                        cGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGt                                                           >  1:1242465/1‑150 (MQ=255)
                                                                                             aTGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCgag                                                      <  2:665151/150‑1 (MQ=255)
                                                                                              tgcccGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCgaga                                                     >  2:15744/4‑150 (MQ=255)
                                                                                              tGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCgaga                                                     >  1:1756864/1‑150 (MQ=255)
                                                                                                aCCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAgatg                                                   >  2:953910/1‑150 (MQ=255)
                                                                                                       cccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGAttgttg                                            >  2:625317/1‑150 (MQ=255)
                                                                                                        gcGTGATGGTGTTCGTTGAAAAGAAAACCTTGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTgg                                           <  1:1393501/149‑1 (MQ=255)
                                                                                                        ccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTgg                                           <  1:1471667/150‑1 (MQ=255)
                                                                                                        ccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTgg                                           <  2:1955594/150‑1 (MQ=255)
                                                                                                          gTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCa                                         <  1:53130/150‑1 (MQ=255)
                                                                                                             aTGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTg                                      <  2:523846/150‑1 (MQ=255)
                                                                                                              tGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGa                                     <  1:127859/150‑1 (MQ=255)
                                                                                                               ggTGTTCGTGGAAAAGAAAAACTTGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTTATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATTATTGTTGGCATTGAc                                    <  1:1504187/150‑1 (MQ=255)
                                                                                                               ggTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGTTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGAc                                    <  1:728240/150‑1 (MQ=255)
                                                                                                                    tCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACtt                               >  1:517410/1‑150 (MQ=255)
                                                                                                                      gTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCa                             >  1:657517/1‑150 (MQ=255)
                                                                                                                      gTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCGACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCa                             <  1:281952/150‑1 (MQ=255)
                                                                                                                       tGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCAt                            <  2:703427/150‑1 (MQ=255)
                                                                                                                       tGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCAt                            <  2:1607000/150‑1 (MQ=255)
                                                                                                                        ggAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGAGCTTCAGGATGACAGCAAGAGCGAGATGATTGTTGGCATTGACCACTTCaag                           >  2:1207233/1‑148 (MQ=255)
                                                                                                                         gaaaagaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa                          >  1:490780/1‑150 (MQ=255)
                                                                                                                         caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa                          >  1:1869319/7‑150 (MQ=255)
                                                                                                                         caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa                          >  1:1676050/7‑150 (MQ=255)
                                                                                                                         caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa                          >  1:1715966/7‑150 (MQ=255)
                                                                                                                          accataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGAc                         >  1:475069/6‑150 (MQ=255)
                                                                                                                          aaaagaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGAc                         >  2:1921643/1‑150 (MQ=255)
                                                                                                                           aaagaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACg                        <  2:54538/150‑1 (MQ=255)
                                                                                                                             agaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGc                      >  2:1094069/1‑150 (MQ=255)
                                                                                                                                 aaaCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCg                  <  1:1102575/150‑1 (MQ=255)
                                                                                                                                   aCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGc                >  1:726619/1‑150 (MQ=255)
                                                                                                                                   aCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGc                <  2:1500322/150‑1 (MQ=255)
                                                                                                                                      ggATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTaa             >  2:828596/1‑150 (MQ=255)
                                                                                                                                            ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTGAGTATTg       >  2:173031/1‑150 (MQ=255)
                                                                                                                                            ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg       >  1:707765/1‑150 (MQ=255)
                                                                                                                                            ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg       >  1:456148/1‑150 (MQ=255)
                                                                                                                                            ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg       >  2:1167394/1‑150 (MQ=255)
                                                                                                                                             cccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGc      <  2:1583012/150‑1 (MQ=255)
                                                                                                                                             cccAACAAGTATTTATGGTGCTACTCCAGTAACAGATTGATCAATCAACATGGCGAGATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGTGATGATTGGTGGCATTGACCACTTCTTGACGGCCTCGGCTAAGTATTGc      <  2:158670/150‑1 (MQ=255)
                                                                                                                                              ccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCg     >  2:578867/1‑150 (MQ=255)
                                                                                                                                                 tCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGATATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGata  <  1:709280/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATA  >  NC_000913/1661801‑1662095

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