Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 15.7% 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.15.7% 593.9 / 115.1 204coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (88/81);  new base . (19/13);  total (109/95)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 204 total aligned reads displayed.

TAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATT  >  NC_000913/3066223‑3066417
                                                                                                  |                                                                                                  
 aGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCATATTGCTGAAGATCGt                                                                                               <  2:3026598/100‑1 (MQ=255)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGt                                                                                               >  2:4317714/1‑100 (MQ=255)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGt                                                                                               >  1:2827120/1‑100 (MQ=25)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCAGCTCTATGCTGGGTTAATGACTCATCCCGTCACAATCTTCAGCAATATGCTGAAGCAt                                                                                                  >  1:366032/1‑97 (MQ=255)
     aCGTCTTCGACTTACACCGTAATTTACTTGCAACCGGCTCGATGCTGCGTTAATGACTCGTCCCGTCACGATCTTCCGCAATATGCTGAAGCa                                                                                                   >  1:1457915/1‑93 (MQ=255)
      cGTCTTCGACTTACACCGTAATTTAATTGCAATCGGCTCGATGCTGGGTTAATTACTCGTCCCGTCACGATGCTTCATCATATTGCTGAAGATCGTGACg                                                                                           <  1:1973236/100‑1 (MQ=21)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCg                                                                                           >  2:3702737/1‑100 (MQ=255)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                          <  1:3045819/100‑1 (MQ=37)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACg                                                                                           <  2:48300/100‑1 (MQ=21)
      cGTCTTCGACTTACACCGTAATTTAAGTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGc                                                                                            <  1:2600064/99‑1 (MQ=255)
       gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                          <  2:2827136/100‑1 (MQ=21)
           tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtata                                                                                      >  2:1511123/1‑100 (MQ=255)
             gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGTTCTTCAGCAACATGCTGAAGCATCGTGGCGGTATAAc                                                                                    <  2:409402/100‑1 (MQ=255)
               cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATCCTGAAGCATCGTGGCGGTATAACga                                                                                  >  2:1170438/1‑100 (MQ=255)
                ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgat                                                                                 <  2:3873734/100‑1 (MQ=255)
                 tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatg                                                                                >  1:3510414/1‑100 (MQ=255)
                 tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGGTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTCTAACgatg                                                                                >  1:155680/1‑100 (MQ=255)
                 tACACCGTAAGTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                  <  1:1334675/98‑1 (MQ=255)
                   cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                                                                              >  2:3445131/1‑100 (MQ=255)
                     ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                            >  2:3635142/1‑100 (MQ=255)
                     ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                            >  1:4166749/1‑100 (MQ=255)
                      cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                           <  1:1881863/100‑1 (MQ=255)
                         aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                        <  2:1445222/100‑1 (MQ=255)
                           tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                      <  2:11044/100‑1 (MQ=255)
                           tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTCGCGGTATAACGATGATTTCTGAtt                                                                      >  1:221568/1‑100 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     >  1:4091941/1‑100 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  1:504461/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  1:2698985/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  2:2694976/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAACGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                       <  1:1013017/98‑1 (MQ=255)
                             taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                    >  1:187095/1‑100 (MQ=255)
                               attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                  <  2:814645/100‑1 (MQ=255)
                               attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                  <  1:3906737/100‑1 (MQ=255)
                                  gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                               <  1:4110616/100‑1 (MQ=255)
                                  gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACCATCTTCAGCAATATGCTGAAGCCTCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                               >  1:956718/1‑100 (MQ=255)
                                   cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                              >  2:1866081/1‑100 (MQ=255)
                                   cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                              <  2:540310/100‑1 (MQ=255)
                                      ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGata                                                           <  2:62403/100‑1 (MQ=255)
                                          cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATTCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                       <  2:1626606/100‑1 (MQ=255)
                                          cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                       <  2:973026/100‑1 (MQ=255)
                                            cGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                     >  2:4030735/1‑100 (MQ=255)
                                             gATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                    <  2:2037890/100‑1 (MQ=255)
                                              aTGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                   <  2:1602465/100‑1 (MQ=255)
                                                gCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                 <  1:124322/100‑1 (MQ=255)
                                                gCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                 >  2:2163516/1‑100 (MQ=255)
                                                 cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCatat                                                >  2:290802/1‑100 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGTTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCatat                                                <  1:3304603/99‑1 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                               <  2:4342435/100‑1 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                               >  2:3709664/1‑100 (MQ=255)
                                                   gggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                              >  2:2595171/1‑100 (MQ=255)
                                                   gggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                              >  1:774102/1‑100 (MQ=255)
                                                    ggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                             >  1:3595722/1‑100 (MQ=255)
                                                        aaTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                         >  2:1357995/1‑100 (MQ=255)
                                                        aaTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                         <  1:3813566/100‑1 (MQ=255)
                                                         aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCCTCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                        >  1:4183553/1‑100 (MQ=255)
                                                         aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                        >  1:1627527/1‑100 (MQ=255)
                                                          tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  1:3774354/100‑1 (MQ=255)
                                                          tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  1:1667581/100‑1 (MQ=255)
                                                          tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  2:3781159/100‑1 (MQ=255)
                                                           gACTCGTCCCGTCACGAGCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCTTGCATATAACTCGGctt                                      >  1:1902690/1‑98 (MQ=255)
                                                             cTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                    >  1:3165410/1‑100 (MQ=255)
                                                             cTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                    >  1:1729574/1‑100 (MQ=255)
                                                             cTCGTCCCGTCACGATCGTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                     <  1:4025943/99‑1 (MQ=255)
                                                               cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  2:1875463/1‑100 (MQ=255)
                                                               cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  1:1561711/1‑100 (MQ=255)
                                                               cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  <  1:4317041/100‑1 (MQ=255)
                                                                gTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                 >  2:1547731/1‑100 (MQ=255)
                                                                 taccgccACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                >  2:2474527/7‑100 (MQ=18)
                                                                  cccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                               >  2:3816751/1‑100 (MQ=255)
                                                                  accGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  1:3199630/2‑97 (MQ=18)
                                                                  accGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                               >  2:924644/2‑100 (MQ=18)
                                                                   ccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                              <  1:1791609/100‑1 (MQ=255)
                                                                   ccGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGACTGTGCGGTg                              >  1:2520488/1‑100 (MQ=18)
                                                                    cgccACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                              >  1:2082894/4‑99 (MQ=25)
                                                                     gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            >  1:4234146/1‑100 (MQ=255)
                                                                     gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            <  2:3419837/100‑1 (MQ=255)
                                                                      tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                           <  2:3746536/100‑1 (MQ=255)
                                                                      ccACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTACTGCATATACATCGGCCTGTGCGGTGCAt                           >  1:888161/2‑100 (MQ=25)
                                                                       cACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                          >  2:4211841/1‑100 (MQ=37)
                                                                       cACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGAGGTGCAt                           >  2:2257032/1‑99 (MQ=25)
                                                                       cACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            >  1:2567881/1‑98 (MQ=255)
                                                                         cgatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  1:1930706/6‑100 (MQ=25)
                                                                         cgatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  2:4234247/6‑100 (MQ=25)
                                                                         cGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  1:1185072/1‑100 (MQ=255)
                                                                         cGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  2:1321632/1‑100 (MQ=255)
                                                                           aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                      >  2:1391851/1‑100 (MQ=255)
                                                                           aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                      <  1:1068295/100‑1 (MQ=255)
                                                                              cTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                    >  1:3822618/1‑100 (MQ=255)
                                                                               ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCAGCATATACATCGGTCTGTTCGGTGCGTCAAAgta                    >  1:2913345/1‑99 (MQ=255)
                                                                                  aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                <  2:3490892/100‑1 (MQ=255)
                                                                                  aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                <  2:329932/100‑1 (MQ=255)
                                                                                    cAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg              <  2:3072050/100‑1 (MQ=255)
                                                                                          tGCTGAAGTATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg         >  2:3184061/1‑100 (MQ=255)
                                                                                          tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg         <  2:1562686/100‑1 (MQ=255)
                                                                                          tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg         >  2:2840481/1‑100 (MQ=255)
                                                                                           gCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc        <  1:3855011/100‑1 (MQ=255)
                                                                                            cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc       >  2:2507818/1‑100 (MQ=255)
                                                                                             tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGTATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg      <  1:2440189/100‑1 (MQ=255)
                                                                                             tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg      <  2:3976642/100‑1 (MQ=255)
                                                                                              gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg     >  1:1375191/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATT  >  NC_000913/3066223‑3066417

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: