Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,231 G→A 19.1% W52* (TGG→TAG)  yqjK → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,2310GA19.1% 330.1 / 54.0 145W52* (TGG→TAG) yqjKuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (47/70);  new base A (13/15);  total (60/85)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

AGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTG  >  NC_000913/3250132‑3250330
                                                                                                   |                                                                                                   
aGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTg                                                                                                     <  2:2541729/100‑1 (MQ=255)
 gCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTgg                                                                                                    >  2:3128288/1‑100 (MQ=255)
 gCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTgg                                                                                                    >  1:648973/1‑100 (MQ=255)
  caacaaCGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTggg                                                                                                   <  2:1391735/100‑1 (MQ=255)
   aacaacGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGgc                                                                                                  >  1:3137557/1‑100 (MQ=255)
     caacGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGgcgc                                                                                                >  1:730291/1‑100 (MQ=255)
      aacGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCt                                                                                               <  1:1364283/100‑1 (MQ=255)
       acGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTg                                                                                              >  2:1741879/1‑100 (MQ=255)
       acGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTg                                                                                              <  1:1622090/100‑1 (MQ=255)
        cGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTgg                                                                                             <  2:1931649/100‑1 (MQ=255)
          gCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGtt                                                                                           >  2:1280211/1‑100 (MQ=255)
           cTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTg                                                                                          <  1:2400298/100‑1 (MQ=255)
           cTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTg                                                                                          <  2:1963598/100‑1 (MQ=255)
            tGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTgg                                                                                         >  1:1562545/1‑100 (MQ=255)
             ggATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGc                                                                                        <  1:3483634/100‑1 (MQ=255)
             ggATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGc                                                                                        >  1:1347153/1‑100 (MQ=255)
              gATCTTTCCGCCAGTCGTCGTGAATGGCTAGAGACAACAGGCGCTTACGTTCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCa                                                                                       <  2:840719/100‑1 (MQ=255)
               atgtttCCGCCAGGCGTCGTGAGTGGCTGGAGACGACAGGCGCTTACGATCGTGGCTGTATTATGCTGCTAAGTCTGCGCGCCTAGGCGCTGGTTGGCAg                                                                                      <  2:711773/97‑1 (MQ=255)
               aTCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAg                                                                                      <  2:2912844/100‑1 (MQ=255)
                   ttCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGc                                                                                  >  2:227880/1‑100 (MQ=255)
                     ccGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGt                                                                                <  1:355427/100‑1 (MQ=255)
                        ccAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGAt                                                                             >  2:3129677/1‑100 (MQ=255)
                        ccAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGAt                                                                             >  2:1423019/1‑100 (MQ=255)
                         cAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTTGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATg                                                                            <  2:3173380/100‑1 (MQ=255)
                         cAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATg                                                                            >  2:2254001/1‑100 (MQ=255)
                          aGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATgg                                                                           <  1:328060/100‑1 (MQ=255)
                          aGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATg                                                                            <  1:1084565/99‑1 (MQ=255)
                           gtcgtcGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTACTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGc                                                                          >  2:1792770/1‑100 (MQ=255)
                               tcgtGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATc                                                                      <  1:1666141/100‑1 (MQ=255)
                                cgtGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCt                                                                     <  1:3038532/100‑1 (MQ=255)
                                cgtGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCt                                                                     <  1:2232561/100‑1 (MQ=255)
                                 gtGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTg                                                                    >  1:1830115/1‑100 (MQ=255)
                                   gAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGa                                                                  >  1:2293007/1‑100 (MQ=255)
                                     aTGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACg                                                                <  2:117296/100‑1 (MQ=255)
                                     aTGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGAGATCTGGACg                                                                >  1:90526/1‑100 (MQ=255)
                                      tGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGa                                                               >  2:276531/1‑100 (MQ=255)
                                       ggCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGAt                                                              >  2:3768663/1‑100 (MQ=255)
                                         cTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTc                                                            >  1:3746798/1‑100 (MQ=255)
                                           ggAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGc                                                          >  1:1976230/1‑100 (MQ=255)
                                           ggAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGc                                                          >  2:1538566/1‑100 (MQ=255)
                                           ggAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGc                                                          <  2:3078134/100‑1 (MQ=255)
                                           ggAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGc                                                          >  1:3560141/1‑100 (MQ=255)
                                                 aacaGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCt                                                    >  1:2552757/1‑100 (MQ=255)
                                                 aacaGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCt                                                    <  1:2736974/100‑1 (MQ=255)
                                                     ggCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAata                                                <  1:3768691/100‑1 (MQ=255)
                                                     ggCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAata                                                <  1:1041304/100‑1 (MQ=255)
                                                      gcgcTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAatat                                               <  2:3458415/100‑1 (MQ=255)
                                                       cgcTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATg                                              >  2:38944/1‑100 (MQ=255)
                                                       cgcTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCATGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATg                                              <  2:3527962/100‑1 (MQ=255)
                                                        gcTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGc                                             <  1:23521/100‑1 (MQ=255)
                                                         cTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCt                                            <  1:3126855/100‑1 (MQ=255)
                                                         cTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCt                                            <  2:841001/100‑1 (MQ=255)
                                                         cTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCt                                            >  2:2460960/1‑100 (MQ=255)
                                                         cTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGATTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCt                                            <  1:702097/100‑1 (MQ=255)
                                                         cTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCt                                            <  2:3079216/100‑1 (MQ=255)
                                                           tACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTgg                                          >  1:3371301/1‑100 (MQ=255)
                                                            aCGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCCTCCTAATATggtggt                                         >  2:945917/1‑100 (MQ=255)
                                                             cGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTAGTc                                        >  1:46850/1‑100 (MQ=255)
                                                               aTCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCg                                      <  1:2162032/100‑1 (MQ=255)
                                                                 cgtcgCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCt                                    <  1:209750/100‑1 (MQ=255)
                                                                  gtcgCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTg                                   <  2:1661356/100‑1 (MQ=255)
                                                                  gtcgCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTg                                   <  2:648966/100‑1 (MQ=255)
                                                                   tcgGTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTgg                                  <  1:982638/100‑1 (MQ=255)
                                                                   tcgCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTgg                                  <  2:3676137/100‑1 (MQ=255)
                                                                   tcgCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCt                                    >  2:165143/1‑98 (MQ=255)
                                                                   tcgCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTgg                                  >  2:2553143/1‑100 (MQ=255)
                                                                    cgCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTggg                                 <  1:3306022/100‑1 (MQ=255)
                                                                     gCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGc                                <  1:2254230/100‑1 (MQ=255)
                                                                     gCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGc                                <  1:3290800/100‑1 (MQ=255)
                                                                     gCTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGc                                <  1:2875814/100‑1 (MQ=255)
                                                                      cTGGAATATGCTGCTAAGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGcc                               <  2:263845/100‑1 (MQ=255)
                                                                           atatGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGAcg                          >  1:2463133/1‑100 (MQ=255)
                                                                             atGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGAcgcg                        <  2:2194135/100‑1 (MQ=255)
                                                                               gctgctAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGt                      <  1:2698308/100‑1 (MQ=255)
                                                                                ctgctAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGtt                     >  1:1919031/1‑100 (MQ=255)
                                                                                 tgctAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGttt                    <  1:3414666/100‑1 (MQ=255)
                                                                                  gctAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGtttt                   <  2:74140/100‑1 (MQ=255)
                                                                                   ctAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTg                  <  2:3329785/100‑1 (MQ=255)
                                                                                    tAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCTGACGCGGTTTTgg                 <  2:2366780/100‑1 (MQ=255)
                                                                                    tAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTgg                 <  1:3075916/100‑1 (MQ=255)
                                                                                    tAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTgg                 <  1:472817/100‑1 (MQ=255)
                                                                                    tAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTgg                 <  1:2770253/100‑1 (MQ=255)
                                                                                     aaGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGc                <  2:1562535/100‑1 (MQ=255)
                                                                                     aaGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGc                >  1:2491096/1‑100 (MQ=255)
                                                                                      aGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGc                <  2:2164150/99‑1 (MQ=255)
                                                                                      aGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCg               <  1:1905712/100‑1 (MQ=255)
                                                                                      aGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCg               <  2:3326438/100‑1 (MQ=255)
                                                                                      aGTCTGCGCTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCg               >  1:1885338/1‑100 (MQ=255)
                                                                                       gTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt              >  1:551661/1‑100 (MQ=255)
                                                                                       gTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGAGGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt              <  2:3152716/100‑1 (MQ=255)
                                                                                         cTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTAt            >  1:2921026/1‑100 (MQ=255)
                                                                                         cTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTAt            >  1:1997131/1‑100 (MQ=255)
                                                                                         cTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTAt            >  2:2269619/1‑100 (MQ=255)
                                                                                         cTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTAt            >  1:1682559/1‑100 (MQ=255)
                                                                                             gcTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAg        <  2:730284/100‑1 (MQ=255)
                                                                                             gcTCCTAGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAg        >  2:2910833/1‑100 (MQ=255)
                                                                                                ccTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAgcgc     <  1:992446/100‑1 (MQ=255)
                                                                                                  tGGGCGCTGGTTGGCAGTAGTGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGACt   >  2:3179886/1‑100 (MQ=255)
                                                                                                   gggCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTg  <  1:3617342/100‑1 (MQ=255)
                                                                                                   gggCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTg  <  2:3629777/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGGGCGCTGGTTGGCAGTAGCGTGATGGCGATCTGGACGATTCGCCATCCTAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTG  >  NC_000913/3250132‑3250330

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