Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 16.0% coding (855/897 nt) ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,3190C.16.0% 372.8 / 68.8 125coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (49/56);  new base . (12/8);  total (61/64)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTC  >  NC_000913/3066221‑3066418
                                                                                                    |                                                                                                   
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCa                                                                                                    >  2:900860/1‑100 (MQ=255)
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCa                                                                                                    >  1:2114828/1‑100 (MQ=255)
 ttAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCAt                                                                                                   >  2:695096/1‑100 (MQ=255)
    gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGCAGCATCGt                                                                                                >  2:122376/1‑100 (MQ=255)
      cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                              <  1:135058/100‑1 (MQ=255)
        ctgcTTCGACGTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  1:230169/97‑1 (MQ=18)
        cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                            <  1:815309/100‑1 (MQ=21)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                           <  2:1000018/100‑1 (MQ=21)
          tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCggn                                                                                          <  2:859666/100‑2 (MQ=255)
          tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                          >  1:1937563/1‑100 (MQ=255)
          tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                           <  2:582758/99‑1 (MQ=21)
          tCTTCAACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGGCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACggg                                                                                          <  1:2400724/100‑2 (MQ=18)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         >  2:542964/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  1:2697425/100‑1 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         >  2:2359525/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  2:2434124/100‑1 (MQ=255)
             tCGACTTACCCCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtata                                                                                       <  1:998419/100‑1 (MQ=255)
             tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtata                                                                                       >  1:32699/1‑100 (MQ=255)
             tCGACTGACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtata                                                                                       <  1:1659223/100‑1 (MQ=255)
                  ttACACCGTAATTTAATAGCAACCAGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgat                                                                                  <  2:33536/100‑1 (MQ=255)
                   tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatg                                                                                 <  1:214287/100‑1 (MQ=255)
                   tACACCGTAATTTAAGTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgat                                                                                  <  1:1094681/99‑1 (MQ=255)
                      tccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAtt                                                                              <  1:1099723/99‑1 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAtta                                                                             >  2:1099725/1‑99 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGGCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                             >  1:2577433/1‑100 (MQ=255)
                          taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                          <  2:2352159/100‑1 (MQ=255)
                           aatttaattGCAACCGGCCTGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                         <  1:2668714/100‑1 (MQ=255)
                            atttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                         >  1:2264020/1‑99 (MQ=255)
                            atttaattGCAACCGGCTAGATGATGGGTTAATGACAAGAACCTTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                          >  1:452688/1‑98 (MQ=255)
                             tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                       >  2:1582517/1‑100 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                      <  2:1115514/100‑1 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                      <  1:884971/100‑1 (MQ=255)
                               taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                     <  2:799540/100‑1 (MQ=255)
                                aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATACCGATGATTTCTGATTGaa                                                                    <  2:1458891/100‑1 (MQ=255)
                                aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                    >  2:720064/1‑100 (MQ=255)
                                 attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                   <  2:1468979/100‑1 (MQ=255)
                                    gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGAGGTATAACGATGATTTCTGATTGAAt                                                                   >  1:191692/1‑97 (MQ=255)
                                      aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                              <  2:2172034/100‑1 (MQ=255)
                                      aaCCGGCTCGATGCTGGGTTAATGACTCGGCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                              <  2:2028561/100‑1 (MQ=255)
                                       aCCGGCTCGATGCTGGGTTAATGACTCGTCCGGTCAAGGC‑CTTCAGCAATA‑TGCTGAGTCATCGTGGCGGTATAACGATGAGTTCTGATTGAATaaa                                                                >  1:1300983/1‑97 (MQ=255)
                                        ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCCGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat                                                             >  1:2598335/1‑99 (MQ=255)
                                            cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                        >  1:2014270/1‑100 (MQ=255)
                                             tCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                       >  2:2748331/1‑100 (MQ=255)
                                              cGATGCTGGGTTACTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                       >  1:995586/1‑99 (MQ=255)
                                              cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                      <  2:1982824/100‑1 (MQ=255)
                                              cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                      >  2:1610554/1‑100 (MQ=255)
                                                aTGGTGGCTTAATCACTCGTCCCGTCACGCT‑CTTCAGTAATA‑TGCTGAAGCATCGTGGCTCTATAACGAGGACTTCTCACTGAATAAACATATTCCTGCa                                                    <  1:1252111/100‑1 (MQ=255)
                                                aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                    <  2:1636385/100‑1 (MQ=255)
                                                aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                    <  1:2397984/100‑1 (MQ=255)
                                                 tGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCat                                                   >  2:745921/1‑100 (MQ=255)
                                                    tGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                                <  1:1220602/100‑1 (MQ=255)
                                                     gggTTAATGACTCGTGCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                        <  1:1105357/91‑1 (MQ=255)
                                                     gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                               >  2:1641882/1‑100 (MQ=255)
                                                       gTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAt                                             >  2:1280315/1‑100 (MQ=255)
                                                         tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                           >  1:1055227/1‑100 (MQ=255)
                                                          aaTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                          >  2:553203/1‑100 (MQ=255)
                                                            tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                        <  1:2362540/100‑1 (MQ=255)
                                                            tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                        <  2:2013992/100‑1 (MQ=255)
                                                              aCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                      <  2:2327967/100‑1 (MQ=255)
                                                               cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                     <  2:2085914/100‑1 (MQ=255)
                                                               cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                     <  2:2180654/100‑1 (MQ=255)
                                                                 cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                   >  1:2549563/1‑100 (MQ=255)
                                                                 cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                   >  2:1549547/1‑100 (MQ=255)
                                                                 cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                   <  2:1368722/100‑1 (MQ=255)
                                                                    cccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                <  1:1665859/100‑1 (MQ=255)
                                                                    cccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                >  1:74097/1‑100 (MQ=255)
                                                                     ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               <  1:229941/100‑1 (MQ=255)
                                                                     ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  2:1820642/1‑100 (MQ=255)
                                                                     ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               <  2:853129/100‑1 (MQ=255)
                                                                      cgccACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACGGTATAACGATGATTTCGGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              >  2:2365299/4‑100 (MQ=21)
                                                                      cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGCCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                 >  1:2069661/4‑97 (MQ=21)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  1:1820635/100‑1 (MQ=255)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  1:1245529/100‑1 (MQ=255)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  1:1868813/100‑1 (MQ=255)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  2:2114835/100‑1 (MQ=255)
                                                                       gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                             <  1:985307/100‑1 (MQ=255)
                                                                         cACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                           >  2:473244/1‑100 (MQ=37)
                                                                           cgatgcTTCAGCATAT‑TGCTGAAG‑ATAGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATACACATAGGCCTGTGCGGTGCATCaa                         >  1:378879/6‑100 (MQ=25)
                                                                           cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGCCGGTATACCGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         >  2:1238676/6‑100 (MQ=25)
                                                                            gatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                        >  1:2798128/5‑100 (MQ=25)
                                                                              tgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                      >  1:360636/3‑100 (MQ=37)
                                                                              t‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                      >  1:2216611/1‑100 (MQ=255)
                                                                                cTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTACTGCATATACATCGGCCTGTGCGGTGCATCAAAata                     >  2:2043987/1‑100 (MQ=255)
                                                                                  tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAGCGGCCTGTGCGGTGCATCAAAATAtt                   <  2:2291519/100‑1 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                   >  1:2546309/1‑99 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                 <  2:2521620/100‑1 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                 <  2:2606709/100‑1 (MQ=255)
                                                                                      cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  2:1463576/100‑1 (MQ=255)
                                                                                      cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               >  2:333278/1‑100 (MQ=255)
                                                                                      cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  1:122375/100‑1 (MQ=255)
                                                                                          a‑tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGtt           >  2:199305/1‑100 (MQ=255)
                                                                                            tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg          >  1:59252/1‑100 (MQ=255)
                                                                                            tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg          <  2:1937571/100‑1 (MQ=255)
                                                                                              cTGAAGCATCGTGGCGGTATAACGATGAGTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc        <  1:2536803/100‑1 (MQ=255)
                                                                                               tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg       <  2:1895897/100‑1 (MQ=255)
                                                                                                gcaGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      >  2:2267107/3‑100 (MQ=255)
                                                                                                gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      >  2:1306206/1‑100 (MQ=255)
                                                                                                gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      >  2:2735412/1‑100 (MQ=255)
                                                                                                gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      >  2:2547220/1‑100 (MQ=255)
                                                                                                  aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATGTCAGCGTTGCCGGAt    >  2:2863631/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTC  >  NC_000913/3066221‑3066418

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