Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,940 T→G 6.2% V462G (GTC→GGC) ‡ 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,9400TG 6.2% 288.9 / 8.6 98V462G (GTC→GGC) ‡ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base G (4/2);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

TAGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAA  >  NC_000913/1661791‑1662084
                                                                                                                                                     |                                                                                                                                                
tAGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGt                                                                                                                                                  >  1:1449591/1‑150 (MQ=255)
  gCTGGGATCTGGGGGTAGAATGGTTCCCGCTGCTCGAGACTCCTGTGTCACCGCTGATTTCCGTCTTTACATGGACAGATGCGATCGACCATGGTACGGAAATGCCCGCGACCCGTGATGGTGGTCGTGGAAAAGAAAAACTGGATGTcc                                                                                                                                                >  1:327781/1‑150 (MQ=255)
   cTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAAATGCCATCGACCATGGTTCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTccc                                                                                                                                               <  2:390585/150‑1 (MQ=255)
   cTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAAATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTccc                                                                                                                                               <  2:1415428/150‑1 (MQ=255)
          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:1081906/1‑150 (MQ=255)
                   gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                               >  1:523846/1‑150 (MQ=255)
                   gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                               >  1:1780143/1‑150 (MQ=255)
                    aaTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGt                                                                                                                              <  2:850020/150‑1 (MQ=255)
                    aaTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGt                                                                                                                              <  1:1993636/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)
                                                       cAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGc                                                                                           <  1:344151/150‑1 (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:1660471/1‑150 (MQ=255)
                                                                          ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                        >  2:1528599/1‑150 (MQ=255)
                                                                          ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                        <  2:1421904/150‑1 (MQ=255)
                                                                           gACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGa                                                                       <  2:1465864/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)
                                                                                                  cGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACCCCCATGAGGTGCTTCAGGATGACAGCAAGt                                                >  2:97543/1‑150 (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                                <  2:1955594/150‑1 (MQ=255)
                                                                                                                  ccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTgg                                <  1:1471667/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:1607000/150‑1 (MQ=255)
                                                                                                                                 tGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCAt                 <  2:703427/150‑1 (MQ=255)
                                                                                                                                  ggAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGAGCTTCAGGATGACAGCAAGAGCGAGATGATTGTTGGCATTGACCACTTCaag                >  2:1207233/1‑148 (MQ=255)
                                                                                                                                   gaaaagaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa               >  1:490780/1‑150 (MQ=255)
                                                                                                                                   caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa               >  1:1676050/7‑150 (MQ=255)
                                                                                                                                   caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa               >  1:1715966/7‑150 (MQ=255)
                                                                                                                                   caccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGa               >  1:1869319/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     <  2:1500322/150‑1 (MQ=255)
                                                                                                                                             aCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGc     >  1:726619/1‑150 (MQ=255)
                                                                                                                                                ggATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTaa  >  2:828596/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TAGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAA  >  NC_000913/1661791‑1662084

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: