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,6690GT97.1% 250.5 / 0.6 70L133I (CTA→ATA) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (34/34);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

GCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCA  >  NC_000913/4616423‑4616907
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
gCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCCTCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGCGGTAACGCCAGAAACGTTTTATTCGCACATCGATCTCGTCTTGTGTTATAAtt                                                                                                                                                                                                                                            >  2:277960/1‑251 (MQ=255)
  ggTTTGCTACATTCACAGATGTTCTTCGCCACTGCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                                                                                                                                                                                          <  2:344604/251‑1 (MQ=255)
    tttGCTTCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                                                                                                                                                                                         >  1:120195/1‑250 (MQ=255)
      tGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAc                                                                                                                                                                                                                                       >  1:94563/1‑250 (MQ=255)
        caaCATTCACAGATGTTCTTCGCCACGTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGGAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                                                                                                                                                                                                    <  1:375352/249‑1 (MQ=255)
        cTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGGTCTCGTCTTGTGTTATAATTCTAACAt                                                                                                                                                                                                                                     >  1:222428/1‑250 (MQ=255)
          aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                                                                                                                                                                                                    <  1:426312/249‑1 (MQ=255)
          aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACGCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                                                                                                                                                                                                    <  2:106055/249‑1 (MQ=255)
            aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGTTTCGGTGAAACCTGCTTCAACTCTGCATCTGGACGTAAGATGTCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGTATATAGGCGTCTGGGAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCTATCTCGTCTTGTGTTATAATTCTAACATACGGt                                                                                                                                                                                                                                <  1:172966/251‑1 (MQ=255)
                  aGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                                                                                                                                                                          <  1:37826/251‑1 (MQ=255)
                          ttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATc                                                                                                                                                                                                                  >  1:154421/1‑251 (MQ=255)
                             gCCACTTCCAGCAGCAGGTCATCAGGGGTTATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACCCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAATGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAg                                                                                                                                                                                                               <  1:201989/251‑1 (MQ=255)
                              ccACTTCCACCAGCAGGTCCTCAGGGGTGATTTCAGGATCGTACATAAAGGACAGGATCGGTGCAACCTGCCTCAACTCTGCATCTGCACGCAAGCTCGCGCGGGTAATGGGCGACTCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGCGGTAACGCCAGAAACGTTTTATTCGAACATCGCTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                                                                                                                                                                              >  2:362640/1‑251 (MQ=255)
                                     cagcagcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCCGTTGcccc                                                                                                                                                                                                        >  2:362187/1‑250 (MQ=255)
                                          gcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                                                                                                                                                    >  2:193850/1‑249 (MQ=255)
                                            agGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                                                                                                                                                                                <  2:92169/251‑1 (MQ=255)
                                                          gATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg                                                                                                                                                                                   >  2:87035/1‑250 (MQ=255)
                                                             ttCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc                                                                                                                                                                                <  1:390426/250‑1 (MQ=255)
                                                               cAGGATCGTAGATAAAGGTCAGTTTCGGTGAAACCTGCTACAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGTCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGGGACCGAc                                                                                                                                                                             <  1:369305/251‑1 (MQ=255)
                                                               cAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc                                                                                                                                                                             >  1:357200/1‑251 (MQ=255)
                                                                      gTAGATAACGATCACGATCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgt                                                                                                                                                                      >  2:159179/1‑251 (MQ=255)
                                                                        aGATAAGGGTCAGGTTCGGTAAATCCTGCTTCAACTCTGCTTCTGCACGTAAGATCGCGCGGGTAATGGGCGAAGCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACAGCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAGCGCATCCAGTTGCCCCAGGTAGACCGGCATGGATGTGACCGACGGTACgtg                                                                                                                                                                     <  2:359347/250‑1 (MQ=255)
                                                                         gATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGCTGTGACCGACGGTACgtggt                                                                                                                                                                   >  1:250339/1‑251 (MQ=255)
                                                                           tAAAGGCCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                                                                                                                                                                 >  1:114251/1‑251 (MQ=255)
                                                                            aaaGGTCAGTTTCGTTGAAACCTGCTTCAACTCTGCAGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTggtggt                                                                                                                                                                <  2:138833/251‑1 (MQ=255)
                                                                                  ccggTTCGGTGAAACCTGCTTCAATTCTGCATCTGAACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGGCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa                                                                                                                                                          <  2:311607/249‑1 (MQ=255)
                                                                                  cAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa                                                                                                                                                           >  1:241290/1‑250 (MQ=255)
                                                                                    ggTTCGGTGACACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGGAAAGAATg                                                                                                                                                        >  2:292069/1‑251 (MQ=255)
                                                                                            gAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGAGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc                                                                                                                                                  >  1:63643/1‑249 (MQ=255)
                                                                                              aaCCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag                                                                                                                                                >  1:301181/1‑249 (MQ=255)
                                                                                               aCCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag                                                                                                                                                >  1:375423/1‑248 (MQ=255)
                                                                                                   gCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGt                                                                                                                                         >  1:48218/1‑251 (MQ=255)
                                                                                                    cTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGCATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTg                                                                                                                                        >  2:204311/1‑251 (MQ=255)
                                                                                                           tctgcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGtcatc                                                                                                                                 >  1:155911/1‑251 (MQ=255)
                                                                                                              gcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGtcatca                                                                                                                                >  1:159499/1‑249 (MQ=255)
                                                                                                                   ggcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCtt                                                                                                                         <  2:154421/250‑1 (MQ=255)
                                                                                                                       cGTAAGATCGCGCGGGTAATGGGCGAATCAGCCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTcgc                                                                                                                     >  1:417042/1‑251 (MQ=255)
                                                                                                                        gTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGCAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg                                                                                                                                                        >  2:63398/1‑215 (MQ=255)
                                                                                                                        gTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg                                                                                                                                                        <  1:63398/215‑1 (MQ=255)
                                                                                                                            gATCGCGCGGGAACTGGGCGAATCAACAGGGCCAATATTGGCGTGCATACCGGCGTCTGGCACGGTTCTATAGGGGTAACGCCAGAAACGATTCATTCGAACAACGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCCTCGATGTGACCGACGGTACGTCGTGGAAAAGAATGGTAAGCCGAGAGATAGCGGCATCAAGATCTTTcacgacgt                                                                                                                >  2:275300/1‑244 (MQ=255)
                                                                                                                                gcgcgGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTtccag                                                                                                            <  1:204311/251‑1 (MQ=255)
                                                                                                                                gcgcgGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCGtccag                                                                                                            <  2:345031/251‑1 (MQ=255)
                                                                                                                                  gcgGGTAATGTGCTAAACAGCCGGGCCGTTATTGGCTTGCATAAAGGCGTCTAGCAGTGTTCTGGCGAGGTAACGCCAGAAACATTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCATTTGCCCCAGGTATACCGGCATCGTTGTGACCGACGGTACGTGTTGGTAAAGAATGGTCTGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCTTtccagc                                                                                                           <  1:300970/250‑1 (MQ=255)
                                                                                                                                               gAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCgg                                                                                              >  2:63831/1‑250 (MQ=255)
                                                                                                                                                       cGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCg                                                                                     <  2:114563/251‑1 (MQ=255)
                                                                                                                                                       cGGGCCGATATTGGCGGGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACAGCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAATCGTTCg                                                                                     <  2:372827/251‑1 (MQ=255)
                                                                                                                                                         gggCCGATATTGGCGTGCATAAAGGCGTCTGGCCGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGc                                                                                    >  1:385723/1‑251 (MQ=255)
                                                                                                                                                              gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGGGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTtccagct                                                                                                          >  2:196280/1‑224 (MQ=255)
                                                                                                                                                              gATATTGGCGTGCATAAAGGCGGCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTtccagct                                                                                                          <  1:196280/224‑1 (MQ=255)
                                                                                                                                                               atatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGaaa                                                                              <  2:94563/251‑1 (MQ=255)
                                                                                                                                                                             aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAAAGCCAGAAACGTTTTATTCGACCATCCATCTCGTCTTGTGTTATAATTCTAACATCCGGTTGAAACAACGCATCCAGTTGCCCCAGGAAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGCACCGTTCGCCAGAAAACGGGCGTAATCgg                                                                >  1:411369/1‑251 (MQ=255)
                                                                                                                                                                                ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGt                                                              >  2:2736/1‑250 (MQ=255)
                                                                                                                                                                                 gCGTCTGGCAGGGTTCTGTCGAGTTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAGCGATGTGACCGACGGTACTTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTGTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTaa                                                            <  2:112490/251‑1 (MQ=255)
                                                                                                                                                                                   gTCTGGCAGAGTTATGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTCACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGCGGTAAAGAATCGTCAGCAGAGCGAGTGCGTCATCAAGATCTCTCGCGCCTTCCAGCTCCGGCCATTCGGTACAGTTCGCCAGAAATCGGGCGTACTCGGGTa                                                             >  2:88277/1‑248 (MQ=255)
                                                                                                                                                                                             gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCgg                                                <  1:335616/251‑1 (MQ=255)
                                                                                                                                                                                                        gTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGgcg                                     >  1:327728/1‑251 (MQ=255)
                                                                                                                                                                                                        gTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGgcg                                     <  2:327728/251‑1 (MQ=255)
                                                                                                                                                                                                               cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGAcc                              <  1:206060/251‑1 (MQ=255)
                                                                                                                                                                                                                 gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCtt                            <  1:38931/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTc                           <  2:375423/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCa                          <  1:63831/251‑1 (MQ=255)
                                                                                                                                                                                                                           aTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATc                  >  1:367244/1‑251 (MQ=255)
                                                                                                                                                                                                                             tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGGTTGTACGGCGCATGACCTTCAAACATATCg                 <  2:385723/250‑1 (MQ=255)
                                                                                                                                                                                                                             tCGAACATCGACCTTGTCTTGTGTTTTAAATCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGCGGACCGCCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCACCAAGATATTTCGCGCCATCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATct                 <  2:222428/250‑2 (MQ=255)
                                                                                                                                                                                                                              cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCAATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCCAACATATCGc                >  1:383661/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGAt            <  1:249517/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTacac       >  1:122325/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTCTTGTGTTATAATTCTAACATACGGTTGCAACCACGCATCCAGTTGCCCCAGGTAGTCCGGCATCGATGTGTCCGTCGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCa  <  2:301181/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCtt    <  1:193850/249‑1 (MQ=255)
                                                                                                                                                                                                                                           cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCa  <  1:137021/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
GCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCA  >  NC_000913/4616423‑4616907

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