Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,669 G→T L133I (CTA→ATA)  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,6690GT98.3% 215.9 / NA 60L133I (CTA→ATA) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (27/32);  minor base C (1/0);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATC  >  NC_000913/4616423‑4616909
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                                                                                                                                                                                              <  1:347758/251‑1 (MQ=255)
      tGCTACATTCACAGATGTTCTTCGCCACTGCCAGCAGCAGGTCATCAGGGGTGATTTCAGGAGCGTAGATAAAGGTCAGGGGCGGTGAAACCTGCTTCAACGCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAc                                                                                                                                                                                                                                         <  2:271069/250‑1 (MQ=255)
          aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTGTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                                                                                                                                                                                                      <  1:53549/249‑1 (MQ=255)
                   gATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                                                                                                                                                                            >  2:397969/1‑250 (MQ=255)
                     tGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                                                                                                                                                                                         >  1:163074/1‑251 (MQ=255)
                         cttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATATGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGATCTGTCGCGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTGGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGcct                                                                                                                                                                                                                     >  1:259400/1‑249 (MQ=255)
                            cGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                                                                                                                                                                                                  <  1:397969/251‑1 (MQ=255)
                                   tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATGTTGGCGTGCATAAAGGCGTCTG‑CAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                                                                                                                                                                           <  2:319620/250‑1 (MQ=255)
                                   tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTG‑CAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                                                                                                                                                                           <  2:319617/250‑1 (MQ=255)
                                      agcagcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                                                                                                                                                                        <  2:112110/251‑1 (MQ=255)
                                             ggTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                                                                                                                                                                                                 <  1:276408/251‑1 (MQ=255)
                                                catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTg                                                                                                                                                                                                                                >  1:120525/1‑217 (MQ=255)
                                                catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTg                                                                                                                                                                                                                                <  2:120525/217‑1 (MQ=255)
                                                    aGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                                                                                                                                                                          >  1:231061/1‑251 (MQ=255)
                                                     ggggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                                                                                                                                                                          >  1:349566/1‑250 (MQ=255)
                                                          gATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg                                                                                                                                                                                     <  2:179858/250‑1 (MQ=255)
                                                                 ggATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGAACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACgg                                                                                                                                                                             >  2:455992/1‑251 (MQ=255)
                                                                         gATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt                                                                                                                                                                     <  1:331636/251‑1 (MQ=255)
                                                                         gATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATAGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAAGTGCCCGACGGTACgtggt                                                                                                                                                                     >  2:398270/1‑251 (MQ=255)
                                                                                   aGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCCTCGATGTGACCGACGGTACGTGGTGGAAAAGAAt                                                                                                                                                           >  2:177049/1‑251 (MQ=255)
                                                                                    ggTTCGGTGAAACCTGCTTCAACTCTGCAGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg                                                                                                                                                          <  2:76442/251‑1 (MQ=255)
                                                                                                         aCTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                                                                                                                                                      >  2:160666/1‑186 (MQ=255)
                                                                                                         aCTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                                                                                                                                                      <  1:160666/186‑1 (MQ=255)
                                                                                                              gcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGtcatca                                                                                                                                  <  2:231061/249‑1 (MQ=255)
                                                                                                               catctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGa                                                                                                                               <  1:354458/251‑1 (MQ=255)
                                                                                                                          aaGATCGCGCGGGTAATGGGCGAATCAGCCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTACAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGcc                                                                                                                    >  2:314573/1‑251 (MQ=255)
                                                                                                                               cgcgcgGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTtcca                                                                                                               <  2:361118/251‑1 (MQ=255)
                                                                                                                                         aTGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGcc                                                                                                      >  1:285888/1‑250 (MQ=255)
                                                                                                                                         aTGGGCGAATCAGACGGGCCGATATTGGCCTGCATAAAGGCGTCTGGCAGGGTTCTGACGAGGTAACGCCAGACACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACACCGCATCCAGATGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGCGGTGGTAAAGAATGGTCAGCAGAGGGAGTGCGTCATCACGATCTTTCGCGCTTTCCAGCTCAAGccc                                                                                                     >  2:384780/1‑250 (MQ=255)
                                                                                                                                         aTGGGCGAACCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCa                                                                                                     <  2:343754/251‑1 (MQ=255)
                                                                                                                                                aaTCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGa                                                                                               >  1:282123/1‑250 (MQ=255)
                                                                                                                                                   cAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAAcc                                                                                            <  1:247830/250‑1 (MQ=255)
                                                                                                                                                     gACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGtt                                                                                         <  1:368816/251‑1 (MQ=255)
                                                                                                                                                             gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTACTGTCGAGGTAACGCCAGAAAAGTTTTATTCGAAAATAGCTCTCGTCTTATGTTATACTTCTAACATACGGTTGCCACAACGCCTCCCGTTGCCCCAGGTAGCCCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCCGCCGAGCGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGa                                                                                  >  1:260722/1‑250 (MQ=255)
                                                                                                                                                                   ggCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCAAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGGCAGCAGAGAGAGTGCGGCATCAAGATCTTTCTCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACggt                                                                           <  1:445162/251‑2 (MQ=255)
                                                                                                                                                                            aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCgg                                                                  >  1:61419/1‑251 (MQ=255)
                                                                                                                                                                             aaGGCGTATGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGGCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGAGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCggg                                                                 >  2:412747/1‑251 (MQ=255)
                                                                                                                                                                              aGGCGTCTGGCAGGGTTCTGTCGATGTAACGCCAGAAACGTTTTATTCGAACATCGATCGCTTCGTGTGTTATAATTCTAACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGATAGACCGGCATAGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGGGAGGGCGTCATCAAGATCTTTCGCGCCTTCTAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCTGGt                                                                <  2:176493/251‑1 (MQ=255)
                                                                                                                                                                                     tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACAt                                                         >  2:170428/1‑251 (MQ=255)
                                                                                                                                                                                        cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAgc                                                      >  1:227034/1‑251 (MQ=255)
                                                                                                                                                                                           ggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGcgcg                                                   >  1:335078/1‑251 (MQ=255)
                                                                                                                                                                                                   cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTAc                                           <  1:163028/251‑1 (MQ=255)
                                                                                                                                                                                                               aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCt                               <  2:365458/251‑1 (MQ=255)
                                                                                                                                                                                                                gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCtt                              <  1:104072/251‑1 (MQ=255)
                                                                                                                                                                                                                gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGGACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCtt                              <  2:37348/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGGTGCAACAACGCATCCAGTTGCCCCAGGGAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCaa                           <  2:349566/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCGTTTTATTCGAACATAGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGATCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCAGCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCaa                           <  2:340152/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAAc                         <  1:131978/251‑1 (MQ=255)
                                                                                                                                                                                                                        ttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACat                       <  1:249274/250‑1 (MQ=255)
                                                                                                                                                                                                                             cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGc                  >  2:174342/1‑250 (MQ=255)
                                                                                                                                                                                                                                aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGc                  >  1:137944/1‑247 (MQ=255)
                                                                                                                                                                                                                                 cATCGATCTCGTCTTGTGTTATAATTCTAAAATACGGTTGCAACAAAGCATCCAGTTGCCCCAGGTAGACCGGCCTCGATGTGACAGACGGTACGTGGTGGAAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGAtt             >  1:301093/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTac           >  2:130836/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gatatcGTCTTGTGTTCTAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGAAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATAGCAGATTacac         >  1:239395/5‑251 (MQ=255)
                                                                                                                                                                                                                                     gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTacac         >  1:239385/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCtt      >  1:163389/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTacac         <  2:163074/246‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACAcc        <  1:85595/247‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCt       <  1:284566/248‑1 (MQ=255)
                                                                                                                                                                                                                                            tCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATc  >  1:135069/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATC  >  NC_000913/4616423‑4616909

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