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.

CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAT  >  NC_000913/3066221‑3066416
                                                                                                    |                                                                                                 
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                                                                  <  2:1575552/100‑1 (MQ=255)
cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGACTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                                                                  >  2:2826199/1‑100 (MQ=255)
    gCCACGTCTTCGACTTACACCGTAATTTAATTTCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGGCACGAGCGTCAGCATATTGCTGAAGATCGTg                                                                                             <  1:283952/100‑1 (MQ=255)
    gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGt                                                                                              >  1:966962/1‑100 (MQ=255)
     ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGCATATTGCTGAAGATCGTg                                                                                             <  2:468313/100‑1 (MQ=37)
     ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTg                                                                                             <  1:361578/100‑1 (MQ=255)
      cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                                                                            >  2:1680560/1‑100 (MQ=255)
        cGTCTTCGACTTACACCGTAATTTAATTGCATCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                                                                            <  1:475449/98‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCgg                                                                                         <  1:1566675/100‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCgg                                                                                         <  2:2626015/100‑1 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCgg                                                                                         >  1:1281626/1‑100 (MQ=255)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCAGATTGCTGAAGATCGTGACgn                                                                                         <  1:1777227/100‑2 (MQ=21)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgn                                                                                         <  2:888123/100‑2 (MQ=21)
         gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                         <  2:600457/100‑1 (MQ=21)
          tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCGGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACggg                                                                                        <  2:425103/100‑2 (MQ=18)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                                                                       >  2:1316738/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                                                                       >  2:1536231/1‑100 (MQ=255)
           cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                                                                       <  2:2860242/100‑1 (MQ=255)
            ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtat                                                                                      >  2:1023866/1‑100 (MQ=255)
              cGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                                                                    >  2:1290937/1‑100 (MQ=255)
                   tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatg                                                                               <  2:2799459/100‑1 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                           >  2:2004823/1‑100 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTACCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                           >  2:373616/1‑100 (MQ=255)
                       ccGTAATTTAATTGCAACCGGCTAGATGATGGGTTAATGTCTCATCCCGTCAGGATCTTCAGCAATATGCTGAAGCATCGTGTCGGTATaaaggtg                                                                               >  1:1608990/1‑91 (MQ=255)
                          taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                        <  1:2543113/100‑1 (MQ=255)
                            atttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                      >  2:2070194/1‑100 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                    <  2:1721671/100‑1 (MQ=255)
                              ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                    <  2:1278883/100‑1 (MQ=255)
                                   tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                               <  1:626218/100‑1 (MQ=255)
                                   tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                               <  2:67073/100‑1 (MQ=255)
                                    gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                              >  2:164732/1‑100 (MQ=255)
                                    gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCCTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                              <  2:1308446/100‑1 (MQ=255)
                                     cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                             >  2:2821034/1‑100 (MQ=255)
                                      aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                            >  2:2708235/1‑100 (MQ=255)
                                       aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat                                                           >  1:2620055/1‑100 (MQ=255)
                                            cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                      <  2:2289973/100‑1 (MQ=255)
                                            cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                      <  2:356216/100‑1 (MQ=255)
                                              cGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                     >  1:1979512/1‑99 (MQ=255)
                                                aTGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                  <  2:406071/100‑1 (MQ=255)
                                                    taggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:1971852/3‑100 (MQ=255)
                                                    tGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:798368/1‑100 (MQ=255)
                                                    tGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                              >  2:2575219/1‑100 (MQ=255)
                                                     gggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                             >  1:1040717/1‑100 (MQ=255)
                                                     gggTTAACGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAATCATCGTGGCGGTATAACGATGATTTCTGATTTAATAAAGATATTCCTGCAtata                                              <  1:2783737/99‑1 (MQ=255)
                                                       gTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAt                                           >  2:223151/1‑100 (MQ=255)
                                                         tAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                         >  2:679317/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGATCTTCAGCAATGTGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:397633/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTTGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:2040070/1‑100 (MQ=255)
                                                           aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:1541886/1‑100 (MQ=255)
                                                            tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:2583769/100‑1 (MQ=255)
                                                            tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:1549132/100‑1 (MQ=255)
                                                            tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                       >  1:1884939/1‑99 (MQ=255)
                                                             gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                      <  2:2399913/99‑1 (MQ=255)
                                                              aCTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                    <  1:976915/100‑1 (MQ=255)
                                                               cTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                   >  1:2226383/1‑100 (MQ=255)
                                                                tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgtg                                  <  1:1769481/100‑1 (MQ=255)
                                                                 cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                 <  2:1575421/100‑1 (MQ=255)
                                                                  gTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAGTAAAGATATTCCTGCATATACATCGGCCtgtg                                  >  1:2501854/1‑98 (MQ=255)
                                                                    cccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                              <  2:1816925/100‑1 (MQ=255)
                                                                     ccGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGACTAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             >  1:1518138/1‑100 (MQ=18)
                                                                      cgccACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            >  2:495612/4‑100 (MQ=21)
                                                                      cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            >  2:2001354/1‑100 (MQ=255)
                                                                      cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            <  2:966965/100‑1 (MQ=255)
                                                                      cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                            <  2:1307277/100‑1 (MQ=255)
                                                                       gccACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             >  1:2077702/3‑98 (MQ=25)
                                                                       gTCACGATCTTCAGGAATATGCTGAAGCATCGTGGCTGTATAACGATGATTTCTGATGGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                             <  1:1291801/98‑1 (MQ=255)
                                                                       gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                           >  1:1908907/1‑100 (MQ=255)
                                                                        ccACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                          >  2:2555156/2‑100 (MQ=21)
                                                                         cACGATGCTTCAGCATATTGCTGAAGCATCGTGGCGGTATCAAGATGATTTCTGATTGAATAAAGATATTCTTGCTTATACATAGGCCTGTGCGGt                              >  1:756468/1‑96 (MQ=25)
                                                                         cACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                         <  2:874644/100‑1 (MQ=255)
                                                                          acgatgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                        >  1:2288707/7‑100 (MQ=37)
                                                                           cgatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                       >  1:135743/6‑100 (MQ=25)
                                                                           cgatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                       >  2:1044014/6‑100 (MQ=25)
                                                                            gatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                      >  2:1726060/5‑100 (MQ=25)
                                                                              tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                    <  2:2757186/100‑1 (MQ=255)
                                                                                 ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                  >  2:252749/1‑100 (MQ=255)
                                                                                  tCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGTGGTGCATCAAAATAtt                 >  1:2478273/1‑100 (MQ=255)
                                                                                  tCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                 >  2:969805/1‑100 (MQ=255)
                                                                                  tCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                 <  1:264884/100‑1 (MQ=255)
                                                                                   cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:113370/1‑100 (MQ=255)
                                                                                   cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:1601418/1‑100 (MQ=255)
                                                                                   cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:254679/1‑100 (MQ=255)
                                                                                   cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  1:255929/1‑100 (MQ=255)
                                                                                   cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                >  2:2272984/1‑100 (MQ=255)
                                                                                    aGCAATATGCTGAAGCATCGTGGCGGTATAACGATTATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  1:2871706/100‑1 (MQ=255)
                                                                                    aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATCTACATCGGCCTGTGCGGTGCATCAAAATATTTc               >  1:1365125/1‑100 (MQ=255)
                                                                                    aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  2:1541892/100‑1 (MQ=255)
                                                                                    aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc               <  2:1553306/100‑1 (MQ=255)
                                                                                     gCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                    >  1:2483827/1‑94 (MQ=255)
                                                                                      cAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg             >  2:2023969/1‑100 (MQ=255)
                                                                                        atatGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg           >  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)
                                                                                                    |                                                                                                 
CTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAT  >  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: