Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 305,259 C→T 5.0% R520H (CGT→CAT)  ecpD ← polymerized tip adhesin of ECP fibers

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913305,2590CT 5.0% 259.3 / 2.1 79R520H (CGT→CAT) ecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base T (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

ACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCG  >  NC_000913/305117‑305393
                                                                                                                                              |                                                                                                                                      
accgaAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGc                                                                                                                                   >  8:57219/1‑148 (MQ=255)
accgaAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGc                                                                                                                                   >  2:283842/1‑148 (MQ=255)
    aaaTCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAAt                                                                                                                              >  1:83354/1‑149 (MQ=255)
     aaTCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCACACCAGCCGTTATCATCTACCGTACGCAGCGAAATg                                                                                                                             >  2:165855/1‑149 (MQ=255)
       tCGCCTGTGTAAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc                                                                                                                           <  1:345410/149‑1 (MQ=255)
         gCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCg                                                                                                                          <  1:299906/148‑1 (MQ=255)
          ccTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGGCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt                                                                                                                        <  6:338011/149‑1 (MQ=255)
          ccTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTAAGGAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt                                                                                                                        <  2:68492/149‑1 (MQ=255)
           cTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTa                                                                                                                       >  5:227162/1‑149 (MQ=255)
              tgtCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATc                                                                                                                     >  2:32137/1‑148 (MQ=255)
                tCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTATGCAGCGAAATGGCGTTATCCAt                                                                                                                  <  7:181113/149‑1 (MQ=255)
                tCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt                                                                                                                  >  5:330408/1‑149 (MQ=255)
                tCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt                                                                                                                  <  2:83354/149‑1 (MQ=255)
                  aaaGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg                                                                                                                <  6:253047/149‑1 (MQ=255)
                  aaaGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGACTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt                                                                                                                 <  1:252898/148‑1 (MQ=255)
                             cAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAAt                                                                                                     >  6:182110/1‑149 (MQ=255)
                                 aTACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa                                                                                                 <  5:218373/149‑1 (MQ=255)
                                 aTACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa                                                                                                 <  3:76172/149‑1 (MQ=255)
                                 aTACAGACCGCGGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa                                                                                                 <  1:32137/149‑1 (MQ=255)
                                  tACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGAc                                                                                                >  4:4089/1‑149 (MQ=255)
                                   aCAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCTCCTGAACATGAATCTCTCCTGAAGCGCTGCCTTCGCCAAACCAGCCGTGATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACt                                                                                               >  6:305177/1‑149 (MQ=255)
                                       acggCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTgg                                                                                           >  8:130413/4‑149 (MQ=255)
                                       aCCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCAAGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTgg                                                                                           >  3:243561/1‑149 (MQ=255)
                                              tctgctgggCCTTAGTTAATGTTACGGCACGTCGCCTGAAGATGAATCTCTCCTGAAGCGCTTACTTCGCCTAACCAGCCGTTATCTTCTACCGTACGCAGGGAAATGGCGTTATCCATTTGAATCGAAAATTTGACTGTGGTCTTAt                                                                                     <  1:222198/142‑1 (MQ=255)
                                              tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTATGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATc                                                                                    >  5:16038/1‑149 (MQ=255)
                                              tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATTGAAAATTTGACTGTGGTCTTAt                                                                                     <  7:64251/148‑1 (MQ=255)
                                              tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTAt                                                                                     <  3:149356/148‑1 (MQ=255)
                                              tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTAt                                                                                     <  7:57219/148‑1 (MQ=255)
                                                  cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATc                                                                                >  7:45400/1‑149 (MQ=255)
                                                      gCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGcc                                                                            >  6:330911/1‑149 (MQ=255)
                                                        cTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGccc                                                                           >  7:342435/1‑148 (MQ=255)
                                                        cTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGccc                                                                           >  3:26047/1‑148 (MQ=255)
                                                            gTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACtt                                                                      >  8:274703/1‑149 (MQ=255)
                                                               aaTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGc                                                                   >  3:138345/1‑149 (MQ=255)
                                                                 tGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTATGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCa                                                                 <  6:16038/149‑1 (MQ=255)
                                                                  gTTACGCCACGTCGCCTGAAGATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCCGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGGGGTCTTATCCCTCTGCCCCACTTCGCCag                                                                >  2:248979/1‑149 (MQ=255)
                                                                   ttACGCAACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaaa                                                                                                      <  5:65762/110‑1 (MQ=255)
                                                                   ttACGCAACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaaa                                                                                                      >  6:65762/1‑110 (MQ=255)
                                                                     aCGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCagag                                                              <  1:283842/148‑1 (MQ=255)
                                                                      cGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAgaga                                                             <  7:71864/148‑1 (MQ=255)
                                                                      cGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGat                                                            <  3:4089/149‑1 (MQ=255)
                                                                         cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa                                                                                                 >  2:46405/1‑109 (MQ=255)
                                                                         cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa                                                                                                 <  1:46405/109‑1 (MQ=255)
                                                                         cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt                                                                                       >  3:125350/1‑119 (MQ=255)
                                                                         cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt                                                                                       <  4:125350/119‑1 (MQ=255)
                                                                         cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATc                                                         <  5:182110/149‑1 (MQ=255)
                                                                          aCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTTCGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCa                                                        <  6:324069/149‑1 (MQ=255)
                                                                               gCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCAAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTccc                                                   >  1:286253/1‑148 (MQ=255)
                                                                                       aTGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTgg                                           <  6:330408/149‑1 (MQ=255)
                                                                                         gAATCTCTCCTGAAGCGCTGACTTCGCCAAACCCGACGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAGCGAAAATTTGACTGTGGTCTTATCCATCTGCCCAACTTCGCCAGAGATATCAGTCCACTGTGTGGTc                                         >  4:175011/1‑149 (MQ=255)
                                                                                          aaTCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTc                                         <  8:183547/148‑1 (MQ=255)
                                                                                          aaTCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTc                                         <  5:330911/148‑1 (MQ=255)
                                                                                           aTCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCaa                                       <  8:342435/149‑1 (MQ=255)
                                                                                                tcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAAc                                      <  4:329230/145‑1 (MQ=255)
                                                                                                  cTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAgc                                <  8:45400/149‑1 (MQ=255)
                                                                                                  cTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAgc                                >  3:184969/1‑149 (MQ=255)
                                                                                                     aaGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCa                             >  2:246068/1‑149 (MQ=255)
                                                                                                        cgcgGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCg                          <  4:138345/149‑1 (MQ=255)
                                                                                                             aCTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTc                       <  4:26047/147‑1 (MQ=255)
                                                                                                              cTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCata                    <  7:274703/149‑1 (MQ=255)
                                                                                                               ttCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCatat                   >  5:193143/1‑149 (MQ=255)
                                                                                                                 cGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCatat                   <  5:316907/147‑1 (MQ=255)
                                                                                                                 cGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATc                  <  7:130413/148‑1 (MQ=255)
                                                                                                                  gCCAAACCAGCCGTTATCATCTACCGTATGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATccc                >  2:95607/1‑149 (MQ=255)
                                                                                                                  gCCAAACCAGCCGTTATCATCTACCGTACGCATCGAAATGGCGTTATCCATTGGAATAGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCAAGAGATATCAGTGCACGCTGTGGTCAAGCACAGCGCATCGGTCATATcc                 <  1:46203/148‑1 (MQ=255)
                                                                                                                  gCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATcc                                                                                   <  1:115322/82‑1 (MQ=255)
                                                                                                                  gCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATcc                                                                                   >  2:115322/1‑82 (MQ=255)
                                                                                                                    cAAACCGGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTc                       <  4:91617/140‑1 (MQ=255)
                                                                                                                    cAAACCGGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTc                       >  3:91617/1‑140 (MQ=255)
                                                                                                                    cAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGc              >  1:24127/1‑149 (MQ=255)
                                                                                                                      aaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgtg                                              >  4:242857/1‑115 (MQ=255)
                                                                                                                      aaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgtg                                              >  4:238407/1‑115 (MQ=255)
                                                                                                                      aaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgtg                                              <  3:242857/115‑1 (MQ=255)
                                                                                                                      aaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgtg                                              <  3:238407/115‑1 (MQ=255)
                                                                                                                       aCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCg                          <  5:233522/134‑1 (MQ=255)
                                                                                                                         cAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGc          >  7:142164/1‑148 (MQ=255)
                                                                                                                                 atcatcTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCg  >  3:297481/1‑148 (MQ=255)
                                                                                                                                       tACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgt                                               >  7:90537/1‑97 (MQ=255)
                                                                                                                                       tACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGgt                                               <  8:90537/97‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                      
ACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCG  >  NC_000913/305117‑305393

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