Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 16.1% 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.1% 348.8 / 66.1 118coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (56/43);  new base . (11/8);  total (67/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 118 total aligned reads displayed.

CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAT  >  NC_000913/3066221‑3066416
                                                                                                    |                                                                                                 
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCa                                                                                                  <  2:1575552/100‑1 (MQ=255)
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGACTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCa                                                                                                  >  2:2826199/1‑100 (MQ=255)
    gCCACGTCTTCGACTTACACCGTAATTTAATTTCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGGCACGAG‑CGTCAGCATAT‑TGCTGAAG‑ATCGTg                                                                                             <  1:283952/100‑1 (MQ=255)
    gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGt                                                                                              >  1:966962/1‑100 (MQ=255)
     ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTg                                                                                             <  2:468313/100‑1 (MQ=37)
     ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTg                                                                                             <  1:361578/100‑1 (MQ=255)
      cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                            >  2:1680560/1‑100 (MQ=255)
        cGTCTTCGACTTACACCGTAATTTAATTGCATCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                            <  1:475449/98‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                         <  1:1566675/100‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                         <  2:2626015/100‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                         >  1:1281626/1‑100 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCAGAT‑TGCTGAAG‑ATCGTGACgn                                                                                         <  1:1777227/100‑2 (MQ=21)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgn                                                                                         <  2:888123/100‑2 (MQ=21)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                         <  2:600457/100‑1 (MQ=21)
          tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCGGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACggg                                                                                        <  2:425103/100‑2 (MQ=18)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                       >  2:1316738/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                       >  2:1536231/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                       <  2:2860242/100‑1 (MQ=255)
            ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtat                                                                                      >  2:1023866/1‑100 (MQ=255)
              cGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATaa                                                                                    >  2:1290937/1‑100 (MQ=255)
                   tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatg                                                                               <  2:2799459/100‑1 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                           >  2:2004823/1‑100 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTACCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                           >  2:373616/1‑100 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTAGATGATGGGTTAATGTCTCATCCCGTCAGGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGTCGGTATaaaggtg                                                                               >  1:1608990/1‑91 (MQ=255)
                          taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                        <  1:2543113/100‑1 (MQ=255)
                            atttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                      >  2:2070194/1‑100 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                    <  2:1721671/100‑1 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                    <  2:1278883/100‑1 (MQ=255)
                                   tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                               <  1:626218/100‑1 (MQ=255)
                                   tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                               <  2:67073/100‑1 (MQ=255)
                                    gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                              >  2:164732/1‑100 (MQ=255)
                                    gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CCTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                              <  2:1308446/100‑1 (MQ=255)
                                     cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                             >  2:2821034/1‑100 (MQ=255)
                                      aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                            >  2:2708235/1‑100 (MQ=255)
                                       aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat                                                           >  1:2620055/1‑100 (MQ=255)
                                            cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                      <  2:2289973/100‑1 (MQ=255)
                                            cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                      <  2:356216/100‑1 (MQ=255)
                                              cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                     >  1:1979512/1‑99 (MQ=255)
                                                aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                  <  2:406071/100‑1 (MQ=255)
                                                    taggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:1971852/3‑100 (MQ=255)
                                                    tGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:798368/1‑100 (MQ=255)
                                                    tGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:2575219/1‑100 (MQ=255)
                                                     gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                             >  1:1040717/1‑100 (MQ=255)
                                                     gggTTAACGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAATCATCGTGGCGGTATAACGATGATTTCTGATTTAATAAAGATATTCCTGCAtata                                              <  1:2783737/99‑1 (MQ=255)
                                                       gTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAt                                           >  2:223151/1‑100 (MQ=255)
                                                         tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                         >  2:679317/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGAT‑CTTCAGCAATG‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:397633/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTTGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:2040070/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:1541886/1‑100 (MQ=255)
                                                            tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:2583769/100‑1 (MQ=255)
                                                            tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:1549132/100‑1 (MQ=255)
                                                            tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:1884939/1‑99 (MQ=255)
                                                             gACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:2399913/99‑1 (MQ=255)
                                                              aCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                    <  1:976915/100‑1 (MQ=255)
                                                               cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                   >  1:2226383/1‑100 (MQ=255)
                                                                tCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgtg                                  <  1:1769481/100‑1 (MQ=255)
                                                                 cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                 <  2:1575421/100‑1 (MQ=255)
                                                                  gTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAGTAAAGATATTCCTGCATATACATCGGCCtgtg                                  >  1:2501854/1‑98 (MQ=255)
                                                                    cccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                              <  2:1816925/100‑1 (MQ=255)
                                                                     ccGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGACTAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             >  1:1518138/1‑100 (MQ=18)
                                                                      cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            >  2:495612/4‑100 (MQ=21)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            >  2:2001354/1‑100 (MQ=255)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            <  2:966965/100‑1 (MQ=255)
                                                                      cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            <  2:1307277/100‑1 (MQ=255)
                                                                       gccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             >  1:2077702/3‑98 (MQ=25)
                                                                       gTCACGAT‑CTTCAGGAATA‑TGCTGAAGCATCGTGGCTGTATAACGATGATTTCTGATGGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             <  1:1291801/98‑1 (MQ=255)
                                                                       gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                           >  1:1908907/1‑100 (MQ=255)
                                                                        ccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                          >  2:2555156/2‑100 (MQ=21)
                                                                         cACGATGCTTCAGCATAT‑TGCTGAAGCATCGTGGCGGTATCAAGATGATTTCTGATTGAATAAAGATATTCTTGCTTATACATAGGCCTGTGCGGt                              >  1:756468/1‑96 (MQ=25)
                                                                         cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                         <  2:874644/100‑1 (MQ=255)
                                                                          acgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                        >  1:2288707/7‑100 (MQ=37)
                                                                           cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                       >  1:135743/6‑100 (MQ=25)
                                                                           cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                       >  2:1044014/6‑100 (MQ=25)
                                                                            gatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                      >  2:1726060/5‑100 (MQ=25)
                                                                              t‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                    <  2:2757186/100‑1 (MQ=255)
                                                                                 ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                  >  2:252749/1‑100 (MQ=255)
                                                                                  tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGTGGTGCATCAAAATAtt                 >  1:2478273/1‑100 (MQ=255)
                                                                                  tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                 >  2:969805/1‑100 (MQ=255)
                                                                                  tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                 <  1:264884/100‑1 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:113370/1‑100 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:1601418/1‑100 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:254679/1‑100 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:255929/1‑100 (MQ=255)
                                                                                   cAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  2:2272984/1‑100 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATTATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  1:2871706/100‑1 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATCTACATCGGCCTGTGCGGTGCATCAAAATATTTc               >  1:1365125/1‑100 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  2:1541892/100‑1 (MQ=255)
                                                                                    aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  2:1553306/100‑1 (MQ=255)
                                                                                     gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                    >  1:2483827/1‑94 (MQ=255)
                                                                                      cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg             >  2:2023969/1‑100 (MQ=255)
                                                                                        ata‑tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg           >  2:2384969/1‑100 (MQ=255)
                                                                                          tgtGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGt          <  2:397635/98‑1 (MQ=255)
                                                                                            tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg        >  1:224954/1‑100 (MQ=255)
                                                                                             gCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc       >  2:57961/1‑100 (MQ=255)
                                                                                              cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc      >  1:2351527/1‑100 (MQ=255)
                                                                                              cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc      >  1:153602/1‑100 (MQ=255)
                                                                                               tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg     <  1:924884/100‑1 (MQ=255)
                                                                                                gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg    >  2:2065482/1‑100 (MQ=255)
                                                                                                 aaGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGGGCGGTGCATCAAAATATTTCAGCGTTGCCGGa   <  1:524713/100‑1 (MQ=255)
                                                                                                  aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg    >  2:28801/1‑98 (MQ=255)
                                                                                                    |                                                                                                 
CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAT  >  NC_000913/3066221‑3066416

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: ‑