Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,943 C→A 6.2% P463H (CCC→CAC) ‡ 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,9430CA 6.2% 278.1 / 5.3 95P463H (CCC→CAC) ‡ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base A (4/2);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 104 total aligned reads displayed.

CTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCG  >  NC_000913/1661794‑1662092
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                      <  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:876627/1‑150 (MQ=255)
                             tGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTa                                                                                                                          >  1:586030/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: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:641160/150‑1 (MQ=255)
                                                              tCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatcaatc                                                                                         <  1:1134030/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)
                                                                                               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: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: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             <  2:1500322/150‑1 (MQ=255)
                                                                                                                                             ggATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTaa          >  2:828596/1‑150 (MQ=255)
                                                                                                                                                   ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg    >  1:456148/1‑150 (MQ=255)
                                                                                                                                                   ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg    >  1:707765/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)
                                                                                                                                                     |                                                                                                                                                     
CTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCG  >  NC_000913/1661794‑1662092

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: