Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,073 T→C E331E (GAA→GAG)  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,0730TC97.7% 132.8 / ‑4.6 44E331E (GAA→GAG) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (19/24);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGG  >  NC_000913/4615828‑4616307
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
gCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacc                                                                                                                                                                                                                                        >  1:182453/1‑250 (MQ=255)
 cATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacac                                                                                                                                                                                                                                       >  2:6503/1‑251 (MQ=255)
  atgatgCGGTGCGACAGTTTGCAGATGGTGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacaca                                                                                                                                                                                                                                      <  1:450103/251‑1 (MQ=255)
  atgatgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACCTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCCATCCCTTCTTTCTcaccc                                                                                                                                                                                                                                       >  1:297033/1‑248 (MQ=255)
    gagtcGGGGCGAAAGTTGGCAGAGGGGGGATCGGATCGGGCTCACCGTCATACGGCAACCGCGCGCCCGGCGTGGGGCCGACATCGGAGCGGATCCCCGACTGTGCGTGTAACATGGCGCGTTCTTGCCAGCCATATTTCACGGGGGTAGTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTGGCGGCGGCTTCTTTACCGTAGCGCGCGTCAATCCCTTCTTTCTcacccat                                                                                                                                                                                                                                     <  2:155583/246‑5 (MQ=255)
         ggggCGACAGTTTGCAGATGGGTGATCGTATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGtt                                                                                                                                                                                                                               <  1:376384/251‑1 (MQ=255)
              gACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCg                                                                                                                                                                                                                          <  1:120382/251‑1 (MQ=255)
                    ttGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTcc                                                                                                                                                                                                                     <  1:6503/250‑1 (MQ=255)
                                       tCTGGCTCATCGCCATACGGCAAACGCGCGACCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGCCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTACGGAGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTAACACCGAAAGTTAACAGCTTCCGCCAGCCAATACCTGCCaa                                                                                                                                                                                                  >  1:329095/1‑250 (MQ=255)
                                          ggCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCTATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCaaa                                                                                                                                                                                                 <  1:86302/248‑1 (MQ=255)
                                                 cgccgcaCGGCAAGCGCGCGCCGGGGGGGGGGCCGATATCGGACCTGTACCCGGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATGTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCGTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGcaca                                                                                                                                                                                       <  2:428869/245‑1 (MQ=255)
                                                           caacccgtcgcCCGGCGTGGGCCCGTTACCGGACCTGCCCCCCTACTGTGCGTTGACCACGCCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATATCCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGTCACAAAACGTTCAg                                                                                                                                                                             <  2:60912/243‑1 (MQ=255)
                                                                             ggtgCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTGTTGCCAGCCATATTTAACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATACCCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttctt                                                                                                                                                           <  1:161965/251‑1 (MQ=255)
                                                                                gCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCa                                                                                                                                                        <  1:113317/251‑1 (MQ=255)
                                                                                  cGATATCGGTACTGAGCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCAcc                                                                                                                                                      <  1:221737/251‑1 (MQ=255)
                                                                                               gATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATtc                                                                                                                                          <  1:76866/250‑1 (MQ=255)
                                                                                                              gtgtAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATGTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAtt                                                                                                                          <  1:308978/251‑1 (MQ=255)
                                                                                                                tgaaCATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAttg                                                                                                                         <  1:4443/248‑1 (MQ=255)
                                                                                                                gtAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACATACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATGCTCAAAGAAGTGTGAttg                                                                                                                         <  1:109849/250‑1 (MQ=255)
                                                                                                                    cATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACCCAGCAAGTTAACCGCTTCCGCCAGCCCATCCATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTAAAAGACGGTTGAttgttgtt                                                                                                                    >  2:118933/1‑249 (MQ=255)
                                                                                                                         cgcgTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGGCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTGCTCACACAGCAAGTTAACCGCTGCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAg‑g                                                                                                              <  1:273861/251‑1 (MQ=255)
                                                                                                                                ttGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTcac                                                                                                       >  1:3699/1‑251 (MQ=255)
                                                                                                                                        caTATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGTTTGTTGATAGAG‑GAATTCACCCCgcgcc                                                                                               >  2:111132/1‑250 (MQ=255)
                                                                                                                                          tatTTAACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATCACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATACCTTATTTCTCACACAGAAAGTTAACCGCATCCGCCAGACCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCACTTCTTTCACACGGAAGCTATTCTCAAAGAAGTGTGATGGTTGATAGAG‑GAATTCacacaacgcat                                                                                              >  1:336017/1‑243 (MQ=255)
                                                                                                                                                           ggCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATctgct                                                                            <  1:116263/251‑1 (MQ=255)
                                                                                                                                                           ggCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATctgct                                                                            <  1:235954/251‑1 (MQ=255)
                                                                                                                                                               aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAAACGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCCGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAACGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                        >  1:372948/1‑251 (MQ=255)
                                                                                                                                                                  aaCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATCCATGCCACACCTTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCCAAGGAGTGTGATTGTTGATA‑AGAGAATTCACACCGCGCATCGATGCTGGAGATCGGCTGCTGGcc                                                                     >  2:440657/1‑250 (MQ=255)
                                                                                                                                                                      ccGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACAAAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCCGGGTTAACCCGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGGGGGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                                                                   >  2:318194/1‑249 (MQ=255)
                                                                                                                                                                         ccAGTTGCCCGCTGATGCGATAACCCACTTCATTTGCGGCGGCCTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGTAAGTTACCCGCTTCCGCCAGCCCCTACATGCCAAACATTGGCACAATACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCACAGAAGGGTGATTGTTGATAGAG‑GAATTCCCACCGCGCCTCGATGATGGCGATCTGCTGCTGGCAGTAGTGc                                                              >  2:226229/1‑251 (MQ=255)
                                                                                                                                                                               gcgcgcTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTag                                                         <  1:371945/250‑1 (MQ=255)
                                                                                                                                                                                cgcgcTGATGCAATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACCTGCCAACCATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATGCTCAAAGAAGTGTGATTGTTTATAGAG‑GCATTCACACCGCGCATGGAAGATGGCGATCTGCTGCTGGCGGTAGTGCGGTag                                                         >  2:264055/1‑249 (MQ=255)
                                                                                                                                                                                   gcTGATGCCATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGACAGTTAACCGCTTCCGCCAGCCCATACCTGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGACGAGTGATTGTTGATTGAG‑GACTTCCCACCGCGCATCGATGATGGCGATCTGCTGCTGGCGGTAGGGCGGTagatagc                                                    >  2:477893/1‑247 (MQ=255)
                                                                                                                                                                                           gATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCAAAAAACGTTTAGGGTTAATCCGCCCTTCTTTCACCAGGACGCTATTCTCAAAGAAGTGTGATTGATGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCCGTAGTGCGGTAGAGTGCGCGTaa                                              >  1:433803/1‑249 (MQ=255)
                                                                                                                                                                                             tAACCCACTTCATTTGCGGCGGCTTCTTTAACGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGCTGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGGGTGCGCGTAaaga                                           >  1:378575/1‑250 (MQ=255)
                                                                                                                                                                                                  caCTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCAt                                     >  1:471331/2‑251 (MQ=255)
                                                                                                                                                                                                               gcggcTTCGTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCt                        <  2:169013/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggcTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCACATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTCGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGcagc                      >  2:436669/1‑248 (MQ=255)
                                                                                                                                                                                                                            cggaGCGCGCGGCAATCCCGTCTTTCTCACACAGCAAGTTAACCGCGGCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATg           <  2:297033/248‑1 (MQ=255)
                                                                                                                                                                                                                             ggaGCGCGCGGCAAGCCCTTCGTTCGCGCACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAGTGGCACAAAACGGTCAGGGTGAATCAGCCCTGCTGTCACCAGGAAGCTATGCTCAAAGAAGTGTGAGTGTTGAGAGAG‑GAATTCACGCCGCGCAGCGATGATGGCGAGCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATgg          <  1:141519/249‑1 (MQ=255)
                                                                                                                                                                                                                             gTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCttt                                                                                                                                                          <  2:287997/108‑1 (MQ=255)
                                                                                                                                                                                                                             gTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCttt                                                                                                                                                          >  1:287997/1‑108 (MQ=255)
                                                                                                                                                                                                                                   cgaggCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGAGCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCACCGATTCGCTGCGCTCGGCAATGGCTTTCa    <  1:383330/248‑1 (MQ=255)
                                                                                                                                                                                                                                     cgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAgg  <  1:454814/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
GCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAG‑GAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGG  >  NC_000913/4615828‑4616307

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