Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,625 T→A 6.7% intergenic (+37/+36) bfr‑I → / ← uvrA bacterioferritin/uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947566,6250TA 6.7% 222.8 / 5.3 75intergenic (+37/+36)bfr‑I/uvrAbacterioferritin/uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGAAC  >  NC_002947/566481‑566773
                                                                                                                                                 |                                                                                                                                                     
cgtgcgGAGCTGGCCGACACCGAAGAAGTACACACCTACTCCTCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCgg                                                                                                                                                    <  1:272222/149‑1 (MQ=255)
    cgcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAc                                                                                                                                                 >  2:100220/1‑147 (MQ=255)
    cgcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCa                                                                                                                                               <  3:315137/149‑1 (MQ=255)
     gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAg                                                                                                                                              <  3:197213/149‑1 (MQ=255)
       aGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTACGGGATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAggg                                                                                                                                            >  1:110053/1‑149 (MQ=255)
       aGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAggg                                                                                                                                            <  4:181341/149‑1 (MQ=255)
           ggCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCtt                                                                                                                                        >  1:292000/1‑149 (MQ=255)
            gCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGCtt                                                                                                                                        <  2:268586/148‑1 (MQ=255)
                acacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                                   <  3:310204/149‑1 (MQ=255)
                  acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                                 <  2:211220/149‑1 (MQ=255)
                  acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                                 >  2:217349/1‑149 (MQ=255)
                  acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                                 >  6:246829/1‑149 (MQ=255)
                  acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                                 <  5:223047/149‑1 (MQ=255)
                    cGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                              <  4:15999/149‑1 (MQ=255)
                    cGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                              <  6:147155/149‑1 (MQ=255)
                    cGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                              >  6:120564/1‑149 (MQ=255)
                     gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgggg                                                                                                                                           >  6:244915/1‑136 (MQ=255)
                     gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgggg                                                                                                                                           <  5:244915/136‑1 (MQ=255)
                               acacCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATa                                                                                                                    <  2:292000/148‑1 (MQ=255)
                                cacCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACt                                                                                                                  <  1:87632/149‑1 (MQ=255)
                                 acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACt                                                                                                                  >  5:314657/1‑148 (MQ=255)
                                  ccTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                                <  2:216964/149‑1 (MQ=255)
                                   cTACTGCCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACGATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAg                                                                                                               <  7:81501/149‑1 (MQ=255)
                                     aCTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGCttttactt                                                                                                                                  <  3:225418/129‑5 (MQ=38)
                                     aCTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGCttttactt                                                                                                                                  >  4:225418/1‑125 (MQ=38)
                                      cTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGcc                                                                                                            <  1:217349/149‑1 (MQ=255)
                                      cTGCCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGcc                                                                                                            >  7:249405/1‑149 (MQ=255)
                                       tGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGAACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGccc                                                                                                           >  5:38603/1‑149 (MQ=255)
                                          gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGACCTa                                                                                                         >  2:265275/1‑149 (MQ=255)
                                              gagaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATcc                                                                                                      <  7:51093/148‑1 (MQ=255)
                                                       cagGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATAGGTACCAGCGGCTTTTTCATTGCGCCGGCAAAACTCAGGCCCTATCCCGCACCAGc                                                                                             >  8:19699/1‑148 (MQ=255)
                                                         gggCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAAt                                                                                          <  8:249405/149‑1 (MQ=255)
                                                               gATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGttt                                                                                     >  4:245577/1‑148 (MQ=255)
                                                                aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                                    <  6:38603/148‑1 (MQ=255)
                                                                 tCAAGGCCATCGTGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATGCACAGACTCAAGTAAAAGCCCCTGGCACCAGCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  6:274885/149‑1 (MQ=255)
                                                                 tCAAGGCCATCGGGCTGGAAAACTACCTTCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCa                                                                                   <  7:19699/148‑1 (MQ=255)
                                                                 tCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  3:245577/149‑1 (MQ=255)
                                                                       ccATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATg                                                                            >  4:225485/1‑149 (MQ=255)
                                                                         aTCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTATCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGc                                                                          >  1:313791/1‑149 (MQ=255)
                                                                                tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAACAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtata                                                                   >  3:218399/1‑149 (MQ=255)
                                                                                 ggAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                                  >  7:115605/1‑149 (MQ=255)
                                                                                    aaaCTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGAc                                                               <  5:202126/149‑1 (MQ=255)
                                                                                     aaCTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCGCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACt                                                              <  6:314657/149‑1 (MQ=255)
                                                                                          ccTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                        <  5:72837/134‑1 (MQ=255)
                                                                                          ccTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTAATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                        >  6:72837/1‑134 (MQ=255)
                                                                                                 tCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGc                                                  >  8:9399/1‑149 (MQ=255)
                                                                                                        tGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATTTGACTGCTTCATCTCGCTCAGctc                                           <  5:120564/149‑1 (MQ=255)
                                                                                                            aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttct                                       >  7:186300/1‑149 (MQ=255)
                                                                                                            aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttct                                       <  1:100220/149‑1 (MQ=255)
                                                                                                             gTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                      <  5:246829/149‑1 (MQ=255)
                                                                                                                aTCCACAGACTCCAGTAAACTCCCCTCGCACCCTCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACAGCTTCATCTCGCTCAGCTCTTCTGgt                                    <  8:236020/148‑1 (MQ=255)
                                                                                                                   cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTc                                                         <  5:62546/124‑1 (MQ=255)
                                                                                                                   cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTc                                                         >  6:62546/1‑124 (MQ=255)
                                                                                                                   cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTAc                                >  7:170432/1‑149 (MQ=255)
                                                                                                                   cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTAc                                <  8:170432/149‑1 (MQ=255)
                                                                                                                    acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTgtt                                    >  2:317110/1‑142 (MQ=255)
                                                                                                                    acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgt                                    <  1:317110/144‑1 (MQ=255)
                                                                                                                      aGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCg                              >  2:210502/1‑148 (MQ=255)
                                                                                                                       gACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCa                            <  4:218399/149‑1 (MQ=255)
                                                                                                                          tCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGGCGCCGGCAATAGTCAGGCCGTATCCCGCTCCAGCATTGGTTTCAGGTAATGGCCGGTATATGACTGGTTCATCTCGCTCTGCTCTTCTGGTGTACCGCAgg                          >  1:235784/1‑149 (MQ=255)
                                                                                                                          tCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAgg                          >  2:138457/1‑148 (MQ=255)
                                                                                                                            aaGTAAAAGCCCCTCGCACAATCGCTACCAGGGGCGTTTTCATTGCGCCGGCAAAACTCAGGCCCTATCCCGCTCCAGCAATGGATTCAGGTAATGCGCGGTATATGACTGCTTCATCTCGCAGAGCGATTCTGGTGTACCGAAGGCaa                       >  8:226936/1‑149 (MQ=255)
                                                                                                                                aaaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGa                    >  6:201548/1‑148 (MQ=255)
                                                                                                                                      cccTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATg                     >  3:164720/1‑141 (MQ=255)
                                                                                                                                      cccTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATg                     <  4:164720/141‑1 (MQ=255)
                                                                                                                                        cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtata                                                                   <  4:173790/93‑1 (MQ=255)
                                                                                                                                        cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtata                                                                   >  3:173790/1‑93 (MQ=255)
                                                                                                                                         tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGc               <  4:267770/144‑1 (MQ=255)
                                                                                                                                         tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGc               >  3:267770/1‑144 (MQ=255)
                                                                                                                                            cACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt       <  4:241992/149‑1 (MQ=255)
                                                                                                                                             aCCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTAATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTg      >  5:298901/1‑149 (MQ=255)
                                                                                                                                               cATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTgg     >  2:143967/1‑148 (MQ=255)
                                                                                                                                                aTCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGaa   <  8:186300/149‑1 (MQ=255)
                                                                                                                                                aTCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGaa   >  6:4503/1‑149 (MQ=255)
                                                                                                                                                aTCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGaa   <  1:138457/149‑1 (MQ=255)
                                                                                                                                                 tCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGAAc  <  7:9399/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGAAC  >  NC_002947/566481‑566773

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: