Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,336,411 2 bp→AT coding (83‑84/1260 nt) murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,4110GA97.3% 128.4 / ‑4.1 38I28N (ATC→AAT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/21);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
*NC_0009133,336,4120AT97.4% 120.8 / NA 38I28N (ATC→AAT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (16/21);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGC  >  NC_000913/3336169‑3336648
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                            
aCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGGCGCGGGCATCAATATGCACAGAGCCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTGCCGCCAGTAGTGCGGCAAAAAGATTAGGCag                                                                                                                                                                                                                                       <  1:57355/251‑1 (MQ=255)
          gACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCTGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTGTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCAttt                                                                                                                                                                                                                             <  2:443834/251‑1 (MQ=255)
                  cgcgcTACCAGCGGCCCCAGCGCCCAGAGAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGcc                                                                                                                                                                                                                     <  1:517496/251‑1 (MQ=255)
                   gcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCg                                                                                                                                                                                                                    <  2:118635/251‑1 (MQ=255)
                           aGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCATCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATtgt                                                                                                                                                                                                            >  2:138134/1‑251 (MQ=255)
                            gCGGCCCCAGCGCCCAGAAAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATtgtg                                                                                                                                                                                                           <  2:130271/251‑1 (MQ=255)
                                          cAGATAGAAGCACGCATGGTTTTAACCAGATCGTCAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAg                                                                                                                                                                                             <  2:325989/251‑1 (MQ=255)
                                             aTAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCtt                                                                                                                                                                                          <  1:164158/251‑1 (MQ=255)
                                               agaagccagCATGGTTTTAACCAGATTGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCTGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCg                                                                                                                                                                                        <  2:123901/243‑1 (MQ=255)
                                                  aGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATTCATTAACGTCGCGGGCAGCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTg                                                                                                                                                                                     <  1:528857/251‑1 (MQ=255)
                                                         aTGGTTTTAACCAGAGCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTCCGTTCTACTTTCGCCCCCAGCTGGCGTAGCAGCTTCATTGATGT‑ATCGACGTCTTGCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCACCTg                                                                                                                                                                              <  2:416167/251‑1 (MQ=255)
                                                              tttAACCATATCGTAAGGGTCGCAGAATACATTTACGTCCCAGGCACCAATATGCACAGACCCATTTCGTTCTTCTTTCGCACCCAGCTGTCTTAGCAGCTTCATTTA‑GTGATTGGCGTCTTTCAGTTTCGGGACTTTCTTGATCGCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGTATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGTAGCTTCGTTGGCCCCTGatca                                                                                                                                                                          <  2:530776/250‑4 (MQ=255)
                                                                tAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACg                                                                                                                                                                        >  2:424887/1‑250 (MQ=255)
                                                                      gATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAAttt                                                                                                                                                                  >  1:364545/1‑250 (MQ=255)
                                                                         cGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATcc                                                                                                                                                              >  2:302734/1‑251 (MQ=255)
                                                                              ggTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTa                                                                                                                                                         <  1:287610/251‑1 (MQ=255)
                                                                                     aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGtt                                                                                                                                                  >  2:198545/1‑251 (MQ=255)
                                                                                      gAATACATTAACGACGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTTTtc                                                                                                                                                 >  1:442091/1‑251 (MQ=255)
                                                                                         tACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtctc                                                                                                                                               >  1:256472/1‑250 (MQ=255)
                                                                                            aTTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGtt                                                                                                                                           >  1:37890/1‑251 (MQ=255)
                                                                                                     gcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTc                                                                                                                                   >  1:380113/1‑250 (MQ=255)
                                                                                                        ggCAGCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCatat                                                                                                                               <  1:47088/251‑1 (MQ=255)
                                                                                                         gCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATc                                                                                                                              >  1:428548/1‑251 (MQ=255)
                                                                                                               atatGCACAGAACCATTACGTTCTACTTTCGCTCCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTAc                                                                                                                        <  1:48102/251‑1 (MQ=255)
                                                                                                                   gCACAGAACCATTACGTTCTACTTTCGCCCCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGGTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGc                                                                                                                    <  1:360566/251‑1 (MQ=255)
                                                                                                                                aCGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCCGTTAACAATTCATATCCGCTACCGGCGAATCGCCCAtc                                                                                                        >  2:63246/1‑249 (MQ=255)
                                                                                                                                   ttCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTc                                                                                                    >  1:522324/1‑251 (MQ=255)
                                                                                                                                       accttCGCACCCAGCTGGCTTAGCAGCTTCATTGAGGT‑ATCGACGTCTTTCAGTTTCGTGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCAGTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTGGTCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCaaa                                                                                                 <  2:68796/247‑1 (MQ=255)
                                                                                                                                       aCTTGCGCCCCCAGCTGGCTTAGCAGCTTCATTGATGC‑ATCGACGGCTTTCAGGTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCa                                                                                                   <  2:133325/248‑1 (MQ=255)
                                                                                                                                            cGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTTCGGCAAAAAGCTTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTAATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGt                                                                                           >  2:519166/1‑251 (MQ=255)
                                                                                                                                                          ttAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcc                                                                             >  2:297397/1‑250 (MQ=255)
                                                                                                                                                           tAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGTGACTGTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcg                                                                             <  1:317430/250‑1 (MQ=255)
                                                                                                                                                                                           aGTTACGGGACGGTCTGGATCCCTTCCGGTTCTTCCGCCAGTAGTTCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCTCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTGCGTTCTCAGTTAACATGTCATATCCGCTTCCGGCGAATCGCCCATAGCTCAAAAGCCGGTCAGTTTGCGATCGCGTGCTCACTCCGCATGGGTATACgtgttgat                                             <  2:559239/251‑8 (MQ=255)
                                                                                                                                                                                                     cGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGc                                    >  2:429448/1‑250 (MQ=255)
                                                                                                                                                                                                     cGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCa                                   <  1:222661/251‑1 (MQ=255)
                                                                                                                                                                                                             acctctACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCg                           >  1:15254/3‑251 (MQ=255)
                                                                                                                                                                                                                               gCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTc         <  2:576370/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gtGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGAATTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGc  >  1:527813/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                            
ACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGT‑ATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGC  >  NC_000913/3336169‑3336648

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