Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,031,753 A→C 8.6% T203P (ACC→CCC) ‡ dppA‑IV → dipeptide transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,031,7530AC 8.6% 200.2 / 9.7 70T203P (ACC→CCC) ‡dppA‑IVdipeptide transport protein
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base C (3/3);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

CACCCTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAA  >  NC_002947/6031620‑6031892
                                                                                                                                     |                                                                                                                                           
ctcccTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGAACGGCCCGGTGCCGATCggtt                                                                                                                              <  7:281689/147‑4 (MQ=255)
 aCCCTCAGCCACCCTGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAATAACCAACCGATCGGCACCGGGCCGTTCgtgt                                                                                                                              <  3:190499/148‑1 (MQ=255)
 aCCCTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCgtgt                                                                                                                              <  6:41014/148‑1 (MQ=255)
   ccTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTc                                                                                                                            >  8:120701/1‑148 (MQ=255)
   ccTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGAACGGCCCGGTGCCGATCggtt                                                                                                                              >  8:281689/1‑143 (MQ=255)
    cTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCa                                                                                                                          <  4:161759/149‑1 (MQ=255)
    cTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGAACGGCCCGGTGCCGATCggttggtt                                                                                                                          <  1:137046/149‑8 (MQ=255)
       aGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCg                                                                                                                       >  1:140794/1‑149 (MQ=255)
           aCCCGGACGCCACCTTCCTGGCCACTCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTc                                                                                                                   <  4:272750/149‑1 (MQ=255)
             ccGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTc                                                                                                                                  >  7:175542/1‑132 (MQ=255)
             ccGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCa                                                                                                                 <  8:24509/149‑1 (MQ=255)
             ccGGAAGCCACCTTACTGGCCACCCTTAGCATGGGGTTCGCCTGGATCTACTCCTCAGAATACGCTGAGAATCTTCTCAAGGCCGGCACACCGGATAAGCTCAACAACCAACCGATCGGCACCTGGCCGTTc                                                                                                                                  <  8:175542/132‑1 (MQ=255)
                   gCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGa                                                                                                           >  3:174479/1‑149 (MQ=255)
                      aCCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGc                                                                                                        >  5:298068/1‑149 (MQ=255)
                        cTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCg                                                                                                      <  7:120701/149‑1 (MQ=255)
                                  aCCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGaa                                                                                                              <  3:93678/131‑1 (MQ=255)
                                  aCCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGATACAGCGTTTCAAGaa                                                                                                              >  4:93678/1‑131 (MQ=255)
                                          cATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGAACACCGGAAAAGCTCAACAACCAACCGATCGGCACAAGTCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGACGTGGAGCGATACCGAGcc                                                                                     >  4:252485/1‑148 (MQ=255)
                                             gggCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAAcc                                                                                 >  3:200200/1‑149 (MQ=255)
                                             gggCTTCGCCTCGATCTACTCCGAAGAATACGATGACAAGCTGCTAAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGCTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAAcc                                                                                 >  7:201414/1‑149 (MQ=255)
                                              ggCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAAccc                                                                                <  3:68043/149‑1 (MQ=255)
                                                cTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCg                                                                               >  3:151463/1‑148 (MQ=255)
                                                 ttCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCg                                                                               <  2:292317/147‑1 (MQ=255)
                                                 ttCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCg                                                                               >  1:292317/1‑147 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCACCAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        >  3:283547/1‑134 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        >  3:278886/1‑134 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        <  4:283547/134‑1 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        <  4:278886/134‑1 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGactac                                                                         <  7:201196/149‑1 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGacta                                                                          >  2:202757/1‑148 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGacta                                                                          <  1:202757/148‑1 (MQ=255)
                                                     ccTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACACGGactac                                                                         >  7:273383/1‑149 (MQ=255)
                                                             tACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        >  3:179262/1‑126 (MQ=255)
                                                             tACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACcgc                                                                                        <  4:179262/126‑1 (MQ=255)
                                                                tCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCaa                                                              >  8:83983/1‑149 (MQ=255)
                                                                  cGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAg                                                                                                                <  7:295120/97‑1 (MQ=255)
                                                                  cGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAg                                                                                                                >  8:295120/1‑97 (MQ=255)
                                                                     aGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcgg                                                         <  4:140393/149‑1 (MQ=255)
                                                                       aaTACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGc                                                                    <  6:157182/136‑1 (MQ=255)
                                                                       aaTACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGc                                                                    >  5:157182/1‑136 (MQ=255)
                                                                       aaTACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCcggc                                                        <  7:263251/148‑1 (MQ=255)
                                                                       aaTACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCCCGCCAACCCCGACTACTTAGCCGGCAAGCccgcg                                                       >  8:263251/1‑149 (MQ=255)
                                                                           cGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaag                                                    <  6:298068/148‑4 (MQ=255)
                                                                            gCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCgg                                                                              <  4:158123/121‑1 (MQ=255)
                                                                            gCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCgg                                                                              >  3:158123/1‑121 (MQ=255)
                                                                            gCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGc                                                                    <  3:270283/131‑1 (MQ=255)
                                                                            gCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGc                                                                    >  4:270283/1‑131 (MQ=255)
                                                                                 cAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGc                                             <  4:174479/149‑1 (MQ=255)
                                                                                           aaGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTc                                                                      <  8:78056/114‑1 (MQ=255)
                                                                                           aaGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTc                                                                      >  7:78056/1‑114 (MQ=255)
                                                                                           aaGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGc                                   >  8:51998/1‑149 (MQ=255)
                                                                                           aaGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGc                                   >  2:181784/1‑149 (MQ=255)
                                                                                           aaGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGc                                   >  4:108901/1‑149 (MQ=255)
                                                                                                  gCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAAc                                                                                  >  8:53956/1‑95 (MQ=255)
                                                                                                  gCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAAc                                                                                  <  7:53956/95‑1 (MQ=255)
                                                                                                  gCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGAGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATcac                             >  3:185086/1‑148 (MQ=255)
                                                                                                    acacCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCaa                                                                                   <  3:294108/92‑1 (MQ=255)
                                                                                                    acacCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCaa                                                                                   >  4:294108/1‑92 (MQ=255)
                                                                                                            aaaGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCGGATCTTCGCCATCCCCCCCGATGcc                   >  8:45691/1‑148 (MQ=255)
                                                                                                                cTCAACAACCGACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGt              >  2:132940/1‑149 (MQ=255)
                                                                                                                 tCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTg             >  6:298873/1‑149 (MQ=255)
                                                                                                                  caacacCCAACCGATCGTCACCGGGCCTTTGGTGTTCCCGCTTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCAAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTgc            <  2:140794/149‑1 (MQ=255)
                                                                                                                   aacaacCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGCCTACTTCGCCGGCAAGCCGGCGGTCGACCCGCGGATCTTCGCCATCACCACCGATGCCACCGGgcg           >  2:274399/1‑149 (MQ=255)
                                                                                                                      tggaacacGAACGGCCCGGTGCCGATCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCt        >  1:7453/8‑149 (MQ=255)
                                                                                                                      aaccaaccGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCt        >  3:303695/1‑149 (MQ=255)
                                                                                                                        ccaaccGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGc      <  5:298873/149‑1 (MQ=255)
                                                                                                                         aacacGAACGGCCCGGTGCCGATCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATg                     >  4:257657/5‑133 (MQ=255)
                                                                                                                         aacacGAACGGCCCGGTGCCGATCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATg                     <  3:257657/129‑1 (MQ=255)
                                                                                                                           accGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAg    >  5:103884/1‑148 (MQ=255)
                                                                                                                            ccGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGaa  <  1:132940/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                           
CACCCTCAGCCACCCGGACGCCACCTTCCTGGCCACCCTGAGCATGGGCTTCGCCTCGATCTACTCCGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAA  >  NC_002947/6031620‑6031892

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: