Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,937,497 A→G 6.1% D164G (GAC→GGC)  PP_4345 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,937,4970AG 6.1% 190.4 / 2.2 65D164G (GAC→GGC) PP_4345GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (2/2);  total (36/29)
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.79e-01

GCGCAGATGTACGAACTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCACCTTGCCGAG  >  NC_002947/4937355‑4937642
                                                                                                                                               |                                                                                                                                                 
gcgcAGATGTACGAACTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTc                                                                                                                                                <  1:322429/146‑1 (MQ=255)
gcgcAGATGTACGAACTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGt                                                                                                                                             <  1:283047/149‑1 (MQ=255)
          aCGAACTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCAt                                                                                                                                   >  8:83155/1‑149 (MQ=255)
              aCTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTc                                                                                                                               >  8:191389/1‑149 (MQ=255)
              aCTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTc                                                                                                                               <  7:191389/149‑1 (MQ=255)
                tctTATCTACCATGGAGCATGCCACCGACCCGCACGACCTTAGCGAACTCAATGAAAAGTTTCATGACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGt                                                                                                                             >  1:177293/1‑149 (MQ=255)
                tctTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTATGGGGAATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCg                                                                                                                              >  1:272443/1‑149 (MQ=255)
                  ttATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGttt                                                                                                                           >  1:153138/1‑149 (MQ=255)
                         ccATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGtttttt                                                                                                                        <  1:254329/145‑1 (MQ=255)
                         ccATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTttgt                                                                                                                      <  8:52005/147‑1 (MQ=255)
                          cATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGa                                                                                                                   >  5:180112/1‑149 (MQ=255)
                                     ccACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGa                                                                                                         <  7:83155/148‑1 (MQ=255)
                                     ccACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGa                                                                                                         <  3:163346/148‑1 (MQ=255)
                                      cACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTttgt                                                                                                                      <  7:190083/134‑1 (MQ=255)
                                      cACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTttgt                                                                                                                      >  8:190083/1‑134 (MQ=255)
                                      cACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGAc                                                                                                        >  2:73383/1‑148 (MQ=255)
                                             ccGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATACAAAATCAAAGCGCCCCAGATTACTTGCATTGATTAATTCCCTTCATGACTCGGTGAATCGCCATCTTCGTTATTTGTTCACAAAACTCGAcagaaagc                                                                                                 >  2:113115/1‑141 (MQ=255)
                                                       ttAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCaa                                                                                      >  6:319492/1‑149 (MQ=255)
                                                           cGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAgg                                                                                  <  2:177293/149‑1 (MQ=255)
                                                                tCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCa                                                                              >  4:161615/1‑148 (MQ=255)
                                                                  aaTGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGGCTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAAc                                                                                     <  6:300915/139‑1 (MQ=255)
                                                                  aaTGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGGCTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAAc                                                                                     <  6:296056/139‑1 (MQ=255)
                                                                  aaTGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGGCTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAAc                                                                                     >  5:300915/1‑139 (MQ=255)
                                                                  aaTGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGGCTCGGTGGATCGACATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAAc                                                                                     >  5:296056/1‑139 (MQ=255)
                                                                  aaTGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATa                                                                            >  5:199228/1‑148 (MQ=255)
                                                                    tGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGa                                                                                                                   >  6:291248/1‑107 (MQ=255)
                                                                    tGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGa                                                                                                                   <  5:291248/107‑1 (MQ=255)
                                                                      aaaaGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAg                                                                        >  3:198541/1‑148 (MQ=255)
                                                                             tcatcaCTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTg                                                                 <  1:129023/148‑1 (MQ=255)
                                                                              catcaCTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAg                                                               >  8:296647/1‑149 (MQ=255)
                                                                               atcaCTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTAGACTATCACAAGAAATCGCAACAGGATCAAAAAAGCATTTTGAgc                                                              >  8:131038/1‑149 (MQ=255)
                                                                               atcaCTTTCTTTATTCAAAACCAAACCGCCCCAGATTACTAGCAATGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAgc                                                              <  8:277820/149‑1 (MQ=255)
                                                                                 caCTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAgcgc                                                            >  7:141612/1‑149 (MQ=255)
                                                                                    tttctttATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTg                                                         >  7:134042/1‑149 (MQ=255)
                                                                                    tttctttATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTg                                                         >  4:111248/1‑149 (MQ=255)
                                                                                            ttCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCTCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAgg                                                 >  4:202533/1‑149 (MQ=255)
                                                                                            ttCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAgg                                                 <  2:272443/149‑1 (MQ=255)
                                                                                            ttCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTACCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAgg                                                 >  3:127393/1‑149 (MQ=255)
                                                                                              cAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGt                                               >  8:141337/1‑149 (MQ=255)
                                                                                               aaaaGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAgg                                                 <  1:50275/146‑1 (MQ=255)
                                                                                               aaaaGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAgg                                                 >  2:50275/1‑146 (MQ=255)
                                                                                                 aagcaaagcGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGAt                                            >  5:169588/1‑149 (MQ=255)
                                                                                                   gcaaagcGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATcg                                          >  2:113750/1‑149 (MQ=255)
                                                                                                    caaagcGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATcgc                                         <  6:199228/149‑1 (MQ=255)
                                                                                                        gcgcCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAg                                     <  8:141612/149‑1 (MQ=255)
                                                                                                             ccAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAgccgcc                                <  2:153138/149‑1 (MQ=255)
                                                                                                               aGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTTACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCACCaa                              <  4:305848/149‑1 (MQ=255)
                                                                                                                         gCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAAtttt                      >  7:69419/1‑148 (MQ=255)
                                                                                                                          cATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTa                    <  6:180112/149‑1 (MQ=255)
                                                                                                                          cATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCATCGAAATTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGAATCATAAAAGCATTTTGAGCGCTTTCGCCAAGGGTGATCGCGAAACCTCCAAGCAAATTTTa                    >  2:283384/1‑149 (MQ=255)
                                                                                                                           aTTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAg                   <  5:319492/149‑1 (MQ=255)
                                                                                                                              gATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGcc                 <  4:198541/148‑1 (MQ=255)
                                                                                                                                ttAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAAttttt                     <  2:102143/142‑1 (MQ=255)
                                                                                                                                ttAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAAAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAAttttt                     >  1:102143/1‑142 (MQ=255)
                                                                                                                                   aGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCTAAGCCGCCAAGCAATTTTTAGCCGAGCAc           >  8:1298/1‑149 (MQ=255)
                                                                                                                                   aGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAc           <  7:296647/149‑1 (MQ=255)
                                                                                                                                   aGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAc           >  3:127954/1‑149 (MQ=255)
                                                                                                                                   aGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAc           >  4:310725/1‑149 (MQ=255)
                                                                                                                                   aGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAc           >  8:5998/1‑149 (MQ=255)
                                                                                                                                    gTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAcc          >  5:272893/1‑149 (MQ=255)
                                                                                                                                       ccTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCACCTTg       <  6:272893/149‑1 (MQ=255)
                                                                                                                                        cTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAcc          <  7:251834/145‑1 (MQ=255)
                                                                                                                                        cTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCAcc          >  8:251834/1‑145 (MQ=255)
                                                                                                                                           cATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCACCTTGCCg    <  1:113750/148‑1 (MQ=255)
                                                                                                                                            aTGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCACCTTGCCGAg  <  4:127954/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCGCAGATGTACGAACTCTTATCTACCATGGAGCATGCCACCGACCCGCACCACCTTAGCGAACTCAATGAAAAGTTTCATCACTTTCTTTATTCAAAAGCAAAGCGCCCCAGATTACTTGCATTGATTAAGTCCCTTCATGACTCGGTGGATCGCCATCTTCGTTTTTTGTTGACAAATCTCGACTATCACAAGAAATCGCAACAGGATCATAAAAGCATTTTGAGCGCTTGCGCCAAGGGTGATCGCGAAGCCGCCAAGCAATTTTTAGCCGAGCACCTTGCCGAG  >  NC_002947/4937355‑4937642

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