Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,535 A→G D177D (GAT→GAC yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,5350AG97.9% 143.5 / ‑4.5 47D177D (GAT→GACyjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/26);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGC  >  NC_000913/4616294‑4616774
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
caagggCTTTCGGGTTAAGGGGTACCGGCGTGCTGCCAGCCCCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGt                                                                                                                                                                                                                                        <  1:394656/247‑1 (MQ=255)
 aaTGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTAACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGGGATTTCAGGATCGTAGATAACGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGt                                                                                                                                                                                                                                        >  2:417630/1‑250 (MQ=255)
 aaTGGATTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACACGCCAGCAGCAGTAAGTTGTAACAGCTCACACTCCCGAAGCCCCCTTTAGCGAAAATTTTATAATGCACCGGACCGTTGGCGATGAGCGAATTGATACATTCACAGATGATCTTAGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTAAGGTTCGGGGAAAAATGCTTCAACTCGGCGTCAGCACGt                                                                                                                                                                                                                                        >  1:280715/1‑250 (MQ=255)
    ggttttCAGGTTAAGGGGTCCCAGCGGGCTGCCACCCCCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGa                                                                                                                                                                                                                                    <  2:363033/248‑1 (MQ=255)
    ggCTTTCAGGTTAAGGCGTTCCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGa                                                                                                                                                                                                                                    <  1:287628/251‑1 (MQ=255)
     gttttCATGTTAAGGCGTCCCAGCGTGGTTCAACCACCCGCCAGCGGCAGGGAGTTGTAACAGCTCACAATCCCGGAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGACAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAt                                                                                                                                                                                                                                   <  2:325802/249‑1 (MQ=255)
             gTTAATGCGTCCCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCggg                                                                                                                                                                                                                           <  2:471331/251‑1 (MQ=255)
               tAAGGCGTACCAGCGTGATGCCACCACCCGCCAGCGGCAGTGCGTTGTAACCGCTCACAATCCCGTAACCCCCTTTTGAGAAAATTTTATCATGCACCGGAACGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGGGATTTCAGGATCGTAGATAAAGGTCAGGTTCGAGGAACCCTGCTTCAACGCTGCGTCTGCACGTAAGATCGCGCggga                                                                                                                                                                                                                          >  2:284047/1‑249 (MQ=255)
                             gTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                                                                                                                                                                                                              >  1:162402/1‑248 (MQ=255)
                                  gcccccccccGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACgg                                                                                                                                                                                                        <  1:177175/242‑1 (MQ=255)
                                       caccCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGata                                                                                                                                                                                                 <  1:120808/251‑1 (MQ=255)
                                                      gggatttttatcagctaccatccccGTAGCCCCCTTTGGTTAAAATTTTATCTGTCCCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTGCTTCGCCAATTCCAGCAGCAGGTCATCAGGGGGGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGGAAGATCGCGAGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAg                                                                                                                                                                                  <  2:182453/229‑1 (MQ=255)
                                                        gAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGTGCCGATATTGGCGTGCATAAAGGc                                                                                                                                                                                <  1:118933/251‑1 (MQ=255)
                                                                          aTCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGCATGGCAGGGTTCTGTc                                                                                                                                                              >  1:43911/1‑251 (MQ=255)
                                                                               gTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGt                                                                                                                                                         >  1:85307/1‑251 (MQ=255)
                                                                                  gCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAAc                                                                                                                                                      >  1:10880/1‑251 (MQ=255)
                                                                                                aaaTTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGtttt                                                                                                                                        <  1:368068/251‑1 (MQ=255)
                                                                                                      tATCATGCCCCGGACCGGTGGCGATGTGCGGTTTGCAACATTCACAGATGTTCGTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCg                                                                                                                                   <  1:264055/250‑1 (MQ=255)
                                                                                                              cACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGCTCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAATGCCAGACACGTTTTATTCGAATATCgg                                                                                                                           >  2:323418/1‑250 (MQ=255)
                                                                                                                ccGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCGGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGAtc                                                                                                                         <  1:477893/251‑1 (MQ=255)
                                                                                                                 cGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGCCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGt                                                                                                                                                         >  2:240530/1‑218 (MQ=255)
                                                                                                                 cGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGt                                                                                                                                                         <  1:240530/218‑1 (MQ=255)
                                                                                                                 cGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGAtct                                                                                                                        >  1:229152/1‑251 (MQ=255)
                                                                                                                   gACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCg                                                                                                                      >  2:154048/1‑251 (MQ=255)
                                                                                                                   gACCGTTGGCGAGTGGCTGTTTGTTAAATTCCCAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCg                                                                                                                      <  2:309852/251‑1 (MQ=255)
                                                                                                                                  gCGGTTTGCTACATTCACAGATGTTCTTCGCCCCTTCCAGCAGCAGGTCATCAGGGTTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGTTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGTCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAAtt                                                                                                       <  1:284353/251‑1 (MQ=255)
                                                                                                                                  gCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAAtt                                                                                                       <  1:111132/251‑1 (MQ=255)
                                                                                                                                                        ggtacttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTGCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAAGATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCaacaac                                                                                  <  1:226229/247‑1 (MQ=255)
                                                                                                                                                         gTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGATGTAACGCCAGAAACGGTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACg                                                                                 <  1:436669/251‑1 (MQ=255)
                                                                                                                                                              tcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATcc                                                                            >  2:65686/1‑251 (MQ=255)
                                                                                                                                                                 ccACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGt                                                                         <  2:390249/251‑1 (MQ=255)
                                                                                                                                                                      tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGc                                                                      <  1:217736/249‑1 (MQ=255)
                                                                                                                                                                                      caGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAACCGTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGACGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCAAGGAACGTTTTATTCGTATATCGGTCTAGTCTTGTGTTAGAATGCTAACATATGGATGCAAAAACGCATCCAGTTGCCCCAGGTAGAACGGCATc                                                    >  2:286334/1‑251 (MQ=255)
                                                                                                                                                                                          ggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg                                                >  1:436665/1‑251 (MQ=255)
                                                                                                                                                                                                   aGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                                       >  1:124378/1‑251 (MQ=255)
                                                                                                                                                                                                     gATCGTAGCTAAAGGTCAGGTTCGGGGAAACCTGGTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAACGGCGTATGGCAGGGTTCTGTCGAGGTAACGCCAGTACCGTTTTTTTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATCTGGTTGTCACAACGCATCCAGTTCACCTAGGTAGACCGTCTTCGTTGTGAACGCCGGt                                     >  2:210595/1‑251 (MQ=255)
                                                                                                                                                                                                     gATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt                                     <  1:154048/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCGTAGATAAAGGTCAGGCTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTACTGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAAACGCTTTATTCGAATATCGATCTCGTCTTGTGTTAGCATTCTAACATATGGTTGCAACAACGCCTCCAGTTGCCCCAGGTAGACAGGCATCGATGTGACCGACGGt                                     >  2:376784/1‑251 (MQ=255)
                                                                                                                                                                                                            gATAAAGTTCAGGTTTGGGGAAACCTGCTTCAACTTTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGTTAACGCCAGAAACGTTTTATTCGAATATCGATCTCTTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt                              <  1:349320/251‑1 (MQ=255)
                                                                                                                                                                                                               aaaGGTCAGGTTCGGGGAAACCTGCTGCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCTGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTGTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCaacaac                                                                                  <  1:248009/196‑1 (MQ=255)
                                                                                                                                                                                                               aaaGGTCAGGATCGGGGAAACCTGCTTCAGCTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCTTCTTGTGTTAGAATTCTAACATATGGTTGCaacaac                                                                                  >  2:248009/1‑196 (MQ=255)
                                                                                                                                                                                                                  ggTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa                        <  1:65686/251‑1 (MQ=255)
                                                                                                                                                                                                                    tCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaa                         >  2:309443/1‑248 (MQ=255)
                                                                                                                                                                                                                     cAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa                      >  2:198090/1‑250 (MQ=255)
                                                                                                                                                                                                                      gggTTCGGGGAAACCTGCTTCAACGCTGCGGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt                    <  2:10880/250‑1 (MQ=255)
                                                                                                                                                                                                                            ggggAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTc                >  1:270695/4‑249 (MQ=255)
                                                                                                                                                                                                                                         tCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGc  <  1:220264/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
CAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGC  >  NC_000913/4616294‑4616774

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