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.

TAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATT  >  NC_000913/3066223‑3066417
                                                                                                  |                                                                                                  
 aGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGC‑ATATTGCTGAAG‑ATCGt                                                                                               <  2:3026598/100‑1 (MQ=255)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGt                                                                                               >  2:4317714/1‑100 (MQ=255)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGt                                                                                               >  1:2827120/1‑100 (MQ=25)
  gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCAGCTCTATGCTGGGTTAATGACTCATCCCGTCACAAT‑CTTCAGCAATA‑TGCTGAAGCAt                                                                                                  >  1:366032/1‑97 (MQ=255)
     aCGTCTTCGACTTACACCGTAATTTACTTGCAACCGGCTCGATGCTGCGTTAATGACTCGTCCCGTCACGAT‑CTTCCGCAATA‑TGCTGAAGCa                                                                                                   >  1:1457915/1‑93 (MQ=255)
      cGTCTTCGACTTACACCGTAATTTAATTGCAATCGGCTCGATGCTGGGTTAATTACTCGTCCCGTCACGATGCTTCATCATAT‑TGCTGAAG‑ATCGTGACg                                                                                           <  1:1973236/100‑1 (MQ=21)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCg                                                                                           >  2:3702737/1‑100 (MQ=255)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAG‑CTTCAGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                          <  1:3045819/100‑1 (MQ=37)
      cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                           <  2:48300/100‑1 (MQ=21)
      cGTCTTCGACTTACACCGTAATTTAAGTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGc                                                                                            <  1:2600064/99‑1 (MQ=255)
       gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                          <  2:2827136/100‑1 (MQ=21)
           tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtata                                                                                      >  2:1511123/1‑100 (MQ=255)
             gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGTT‑CTTCAGCAACA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                    <  2:409402/100‑1 (MQ=255)
               cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TCCTGAAGCATCGTGGCGGTATAACga                                                                                  >  2:1170438/1‑100 (MQ=255)
                ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgat                                                                                 <  2:3873734/100‑1 (MQ=255)
                 tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatg                                                                                >  1:3510414/1‑100 (MQ=255)
                 tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGGTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTCTAACgatg                                                                                >  1:155680/1‑100 (MQ=255)
                 tACACCGTAAGTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACga                                                                                  <  1:1334675/98‑1 (MQ=255)
                   cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatgat                                                                              >  2:3445131/1‑100 (MQ=255)
                     ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                            >  2:3635142/1‑100 (MQ=255)
                     ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                            >  1:4166749/1‑100 (MQ=255)
                      cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                           <  1:1881863/100‑1 (MQ=255)
                         aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                        <  2:1445222/100‑1 (MQ=255)
                           tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                      <  2:11044/100‑1 (MQ=255)
                           tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTCGCGGTATAACGATGATTTCTGAtt                                                                      >  1:221568/1‑100 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     >  1:4091941/1‑100 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  1:504461/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  1:2698985/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                     <  2:2694976/100‑1 (MQ=255)
                            ttaattGCAACCGGCTCGATGCTGGGTTAACGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                       <  1:1013017/98‑1 (MQ=255)
                             taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                    >  1:187095/1‑100 (MQ=255)
                               attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                  <  2:814645/100‑1 (MQ=255)
                               attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                  <  1:3906737/100‑1 (MQ=255)
                                  gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                               <  1:4110616/100‑1 (MQ=255)
                                  gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACCAT‑CTTCAGCAATA‑TGCTGAAGCCTCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                               >  1:956718/1‑100 (MQ=255)
                                   cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                              >  2:1866081/1‑100 (MQ=255)
                                   cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                              <  2:540310/100‑1 (MQ=255)
                                      ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGata                                                           <  2:62403/100‑1 (MQ=255)
                                          cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TTCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                       <  2:1626606/100‑1 (MQ=255)
                                          cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                       <  2:973026/100‑1 (MQ=255)
                                            cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                     >  2:4030735/1‑100 (MQ=255)
                                             gATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                    <  2:2037890/100‑1 (MQ=255)
                                              aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                   <  2:1602465/100‑1 (MQ=255)
                                                gCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                 <  1:124322/100‑1 (MQ=255)
                                                gCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                 >  2:2163516/1‑100 (MQ=255)
                                                 cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCatat                                                >  2:290802/1‑100 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGTT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCatat                                                <  1:3304603/99‑1 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                               <  2:4342435/100‑1 (MQ=255)
                                                  tGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                               >  2:3709664/1‑100 (MQ=255)
                                                   gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                              >  2:2595171/1‑100 (MQ=255)
                                                   gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                              >  1:774102/1‑100 (MQ=255)
                                                    ggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                             >  1:3595722/1‑100 (MQ=255)
                                                        aaTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                         >  2:1357995/1‑100 (MQ=255)
                                                        aaTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                         <  1:3813566/100‑1 (MQ=255)
                                                         aTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCCTCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                        >  1:4183553/1‑100 (MQ=255)
                                                         aTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                        >  1:1627527/1‑100 (MQ=255)
                                                          tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  1:3774354/100‑1 (MQ=255)
                                                          tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  1:1667581/100‑1 (MQ=255)
                                                          tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                       <  2:3781159/100‑1 (MQ=255)
                                                           gACTCGTCCCGTCACGAG‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCTTGCATATAACTCGGctt                                      >  1:1902690/1‑98 (MQ=255)
                                                             cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                    >  1:3165410/1‑100 (MQ=255)
                                                             cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                    >  1:1729574/1‑100 (MQ=255)
                                                             cTCGTCCCGTCACGAT‑CGTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                     <  1:4025943/99‑1 (MQ=255)
                                                               cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  2:1875463/1‑100 (MQ=255)
                                                               cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  1:1561711/1‑100 (MQ=255)
                                                               cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  <  1:4317041/100‑1 (MQ=255)
                                                                gTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                 >  2:1547731/1‑100 (MQ=255)
                                                                 taccgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                >  2:2474527/7‑100 (MQ=18)
                                                                  cccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                               >  2:3816751/1‑100 (MQ=255)
                                                                  accGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                  >  1:3199630/2‑97 (MQ=18)
                                                                  accGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                               >  2:924644/2‑100 (MQ=18)
                                                                   ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                              <  1:1791609/100‑1 (MQ=255)
                                                                   ccGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGACTGTGCGGTg                              >  1:2520488/1‑100 (MQ=18)
                                                                    cgccACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                              >  1:2082894/4‑99 (MQ=25)
                                                                     gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            >  1:4234146/1‑100 (MQ=255)
                                                                     gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            <  2:3419837/100‑1 (MQ=255)
                                                                      tCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                           <  2:3746536/100‑1 (MQ=255)
                                                                      ccACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTACTGCATATACATCGGCCTGTGCGGTGCAt                           >  1:888161/2‑100 (MQ=25)
                                                                       cACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                          >  2:4211841/1‑100 (MQ=37)
                                                                       cACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGAGGTGCAt                           >  2:2257032/1‑99 (MQ=25)
                                                                       cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                            >  1:2567881/1‑98 (MQ=255)
                                                                         cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  1:1930706/6‑100 (MQ=25)
                                                                         cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  2:4234247/6‑100 (MQ=25)
                                                                         cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  1:1185072/1‑100 (MQ=255)
                                                                         cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                        >  2:1321632/1‑100 (MQ=255)
                                                                           aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                      >  2:1391851/1‑100 (MQ=255)
                                                                           aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                      <  1:1068295/100‑1 (MQ=255)
                                                                              cTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                    >  1:3822618/1‑100 (MQ=255)
                                                                               ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCAGCATATACATCGGTCTGTTCGGTGCGTCAAAgta                    >  1:2913345/1‑99 (MQ=255)
                                                                                  aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                <  2:3490892/100‑1 (MQ=255)
                                                                                  aGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                <  2:329932/100‑1 (MQ=255)
                                                                                    cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg              <  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)
                                                                                                  |                                                                                                  
TAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATT  >  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: ‑