Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,337,390 C→T K180K (AAG→AAA)  adiY ← DNA‑binding transcriptional activator AdiY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,337,3900CT92.5% 126.5 / 2.4 40K180K (AAG→AAA) adiYDNA‑binding transcriptional activator AdiY
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (24/13);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-02

AGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGG  >  NC_000913/4337146‑4337628
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
aGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGGCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGTTGATTATTTGGCTATAACTTGTGTTTTCGCTTTTCAACTTTTTTTTTAAc                                                                                                                                                                                                                                          <  1:137611/251‑1 (MQ=255)
   tctTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGGGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAgc                                                                                                                                                                                                                                       >  2:59201/1‑250 (MQ=255)
     ttCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACt                                                                                                                                                                                                                                     >  2:507473/1‑251 (MQ=255)
         ttAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTggg                                                                                                                                                                                                                                 >  1:555836/1‑251 (MQ=255)
               ttAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTa                                                                                                                                                                                                                              >  2:501289/1‑248 (MQ=255)
                            gACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTgg                                                                                                                                                                                                              <  1:218021/251‑1 (MQ=255)
                                                       aTAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTtgt                                                                                                                                                                                   >  1:427154/1‑251 (MQ=255)
                                                          aTGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATCTTTTTACCGTCCATCCTTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGa                                                                                                                                                                                 >  1:242520/1‑250 (MQ=255)
                                                              aGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATc                                                                                                                                                                            >  2:516260/1‑251 (MQ=255)
                                                                cggcgTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGc                                                                                                                                                                          >  2:216280/1‑251 (MQ=255)
                                                                 ggcgTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCt                                                                                                                                                                         >  1:502292/1‑251 (MQ=255)
                                                                       aTACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATACCGCTTTCaat                                                                                                                                                                   <  1:251078/251‑1 (MQ=255)
                                                                          ccGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCaataat                                                                                                                                                                <  1:216280/251‑1 (MQ=255)
                                                                           cGTAGAAGTCTTTAAAGACAGACATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAAtt                                                                                                                                                               >  2:546026/1‑251 (MQ=255)
                                                                            gTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCTTACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAAttt                                                                                                                                                              >  1:564175/1‑251 (MQ=255)
                                                                                 aGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATaa                                                                                                                                                         <  2:502292/251‑1 (MQ=255)
                                                                                  gTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATaaa                                                                                                                                                        <  1:511946/251‑1 (MQ=255)
                                                                                      ttAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTGTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAAc                                                                                                                                                       <  1:65941/248‑1 (MQ=255)
                                                                                         aaGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTc                                                                                                                                                 >  1:180726/1‑251 (MQ=255)
                                                                                                         acgacgTACTGTTTGAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCTTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTTTCACTTACACAATTACGt                                                                                                                                 <  2:271749/251‑1 (MQ=255)
                                                                                                           gacgTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCCACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTAAATAATTTGATAAACGCTGTCACTTACACAATTACGTa                                                                                                                                >  2:208730/1‑250 (MQ=255)
                                                                                                                               ggACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAAGTACGTAACATATACATCATTAATGaa                                                                                                            >  2:390377/1‑250 (MQ=255)
                                                                                                                                gACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAAc                                                                                                          <  2:180726/251‑1 (MQ=255)
                                                                                                                                           cTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTAc                                                                                               >  1:3597/1‑251 (MQ=255)
                                                                                                                                                agaTATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTct                                                                                          >  2:8529/1‑250 (MQ=255)
                                                                                                                                                agaTATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTAACATACTAAGAATCCAGTCTTTGTGAATATCGATTTCAATAATTTGATAAAAGCTGTCACTTACACAATTACGTAACATATACCTCATTAATGAAAAGAATTTTTTACTGTCg                                                                                          >  1:581214/1‑251 (MQ=255)
                                                                                                                                                    atTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGaa                                                                                      <  1:59201/251‑1 (MQ=255)
                                                                                                                                                                 atcatTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCAc                                                                         >  1:408990/1‑251 (MQ=255)
                                                                                                                                                                      tAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTg                                                                     <  1:131521/250‑1 (MQ=255)
                                                                                                                                                                       aaTTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGaa                                                                   <  2:555836/251‑1 (MQ=255)
                                                                                                                                                                                      tAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCgtgt                                                    >  2:508278/1‑251 (MQ=255)
                                                                                                                                                                                                          gATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGATAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAACCGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTc                                  >  2:145750/1‑249 (MQ=255)
                                                                                                                                                                                                                 tGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTg                          <  1:507473/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCg                        >  1:110406/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggCTATAACTGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCg                        <  1:468009/251‑1 (MQ=255)
                                                                                                                                                                                                                           tGGTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGACTTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGcatactg               >  2:74180/1‑247 (MQ=255)
                                                                                                                                                                                                                            ggTGTTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATg               <  1:516260/250‑1 (MQ=255)
                                                                                                                                                                                                                              tgtTTTCGCTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCtg            >  2:567086/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTTTTCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTgcgc    >  1:3818/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tttCAACTTTTTTTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCgg  >  2:359294/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
AGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGG  >  NC_000913/4337146‑4337628

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