Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,967 Δ1 bp 35.6% intergenic (+205/‑341) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,9670A.35.6% 118.9 / 75.2 59intergenic (+205/‑341)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (14/23);  new base . (12/9);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAA  >  NC_002947/5988826‑5989116
                                                                                                                                             |                                                                                                                                                     
cTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGaaacaaac                                                                                                                                              <  2:687201/150‑1 (MQ=2)
 ttCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                                              <  2:198542/150‑1 (MQ=255)
 ttCGCGGGTAAACCCACTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAATAAGGCGGTCGCGCGCGCCTGGCACAGGAGTGATTGGCCAGaaacaaac                                                                                                                                              <  2:487315/150‑1 (MQ=255)
   cgcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCGGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGAAGCTGGAGGGATTGGCCAGAAACAAACGt                                                                                                                                            <  1:29928/150‑1 (MQ=37)
       ggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGa                                                                                                                                        >  1:1043128/1‑150 (MQ=255)
       ggTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAg                                                                                                                                       >  1:271914/1‑150 (MQ=2)
             cccGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGtt                                                                                                                                 >  2:1075815/1‑150 (MQ=2)
                gCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTAGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                                                                                                          <  2:355433/139‑1 (MQ=255)
                gCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTAGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                                                                                                          >  1:355433/1‑139 (MQ=255)
                  tCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAAGAGCGGGTTTAcc                                                                                                                             <  2:214173‑M2/150‑3 (MQ=255)
                   ccTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCcg                                                                                                                           <  2:388972/150‑1 (MQ=2)
                    cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccg                                                                                                                           <  2:729796‑M2/150‑5 (MQ=255)
                     tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgc                                                                                                                          <  1:664255‑M2/150‑6 (MQ=255)
                       cGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcga                                                                                                                        <  2:962665‑M2/150‑8 (MQ=255)
                        gtttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaa                                                                                                                       >  2:683149‑M2/1‑142 (MQ=255)
                         tttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaag                                                                                                                      <  2:34171‑M2/150‑10 (MQ=255)
                         tttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaag                                                                                                                      <  1:319429‑M2/150‑10 (MQ=255)
                         tttgtttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTCTACCCGCGAAgc                                                                                                                     >  2:14525/1‑150 (MQ=2)
                             tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCAACATGAGGCATGTATAGATATCGCGTCGTACAAACAAGGCGGACGCGCGCGACTGGTGCAGGAGGGACTGGCCGCAACAAAAGTCGGAGCGGGTTTACCCGCGAAGccccg                                                                                                                 >  2:38882/1‑150 (MQ=255)
                                    ccAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcggg                                                                                                           <  1:1083987‑M2/150‑21 (MQ=255)
                                     cAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgagcgccgcccgcgcg                                                                                                          <  2:523963‑M2/150‑22 (MQ=255)
                                      aattATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgggcgg                                                                                                        >  2:649750/3‑150 (MQ=14)
                                          aTTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgc                                                                                                     >  1:729796‑M2/1‑124 (MQ=255)
                                              cTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcga                                                                                                 <  2:805514‑M2/150‑31 (MQ=255)
                                                gtgGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACTTGGCTGCATGATGCCCCTTTCAACGGAGGGATTTGCCAGAAACCAACGTAGGAGTTTGTTTATCCGCGAGCGCCGCCCGCGCGGCGCttcgcgg                                                                                              <  1:523963/150‑7 (MQ=255)
                                                           aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCCGAAACAAACGTAGGAGCGGGAGTACCCGCAGAGAGCCGCGCGGGCGGCGCTCGCTCCCATGGGCGcc                                                                                    >  1:284571/1‑150 (MQ=21)
                                                           aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGcc                                                                                    >  1:916236/1‑150 (MQ=255)
                                                             gcgAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCa                                                                                  >  2:728139/1‑150 (MQ=255)
                                                                  cacccTGGAGAACGCCTCGATATAACGTCGTACGAACAAGGCTGTGGCGCGCGCCTGGCACCGTAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCCCCGCTCGGGCGCCGCTCGATCTCATGGGCGCCCAAAAAc                                                                             >  2:68029/3‑150 (MQ=25)
                                                                     ccTGTCTCACGCCTCGATATCGCGGCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTAt                                                                          <  2:284571/150‑1 (MQ=255)
                                                                               gCCTCGATATCGCGTCGTACAAACAAGGCGGGCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCAc                                                                <  1:34171/150‑1 (MQ=255)
                                                                                  tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcg                                                                                                         >  1:707172‑M2/1‑83 (MQ=255)
                                                                                  tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcg                                                                                                         <  2:707172‑M2/105‑23 (MQ=255)
                                                                                    gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaa                                                          <  2:1002070‑M2/150‑70 (MQ=255)
                                                                                     atatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGaa                                                          <  2:916236/150‑1 (MQ=255)
                                                                                     atatCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaa                                                         >  1:221498‑M2/1‑80 (MQ=255)
                                                                                      tatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCCCCAGAGGTATTGGCCAGAAACAAACGTAGGCGCCCGATTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAg                                                         <  1:805514/150‑1 (MQ=25)
                                                                                      tatCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaat                                                        >  2:143507‑M2/1‑79 (MQ=255)
                                                                                                tACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGc                                                <  1:593483/149‑1 (MQ=255)
                                                                                                tACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGc                                                >  2:593483/1‑149 (MQ=255)
                                                                                                   aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAANAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAAtt                                            >  1:8526/1‑150 (MQ=255)
                                                                                                   aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAAtt                                            >  2:716494/1‑150 (MQ=255)
                                                                                                   aaacaaGGCGGTCGCGCGCGACTGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccg                                                                        <  2:8695‑M2/120‑56 (MQ=255)
                                                                                                   aaacaaGGCGGTCGCGCGCGACTGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccg                                                                        >  1:8695‑M2/1‑65 (MQ=255)
                                                                                                     acaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTc                                          <  2:1043128/150‑1 (MQ=255)
                                                                                                     acaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTAGCCGAAACAAACGTAGGAGCGGGTTTAcccactaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattcccacgccgtccac                                         <  2:266792‑M2/150‑87 (MQ=255)
                                                                                                      caaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTcc                                         <  2:319429/150‑1 (MQ=255)
                                                                                                      caaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTcc                                         <  2:1083987/150‑1 (MQ=255)
                                                                                                       aaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTg                                       >  2:648749/1‑150 (MQ=255)
                                                                                                       aaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGCGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactc                                       <  2:407631‑M2/150‑89 (MQ=255)
                                                                                                         ggCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGa                                      >  1:638155/1‑150 (MQ=255)
                                                                                                         ggCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaa                                     <  1:951863‑M2/150‑91 (MQ=255)
                                                                                                                    gcgaCTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgc                          >  1:1075815‑M2/1‑49 (MQ=255)
                                                                                                                          ggtgcaggagGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgcccaccacccaaaattccaacgccgcccc                    >  2:271914‑M2/5‑43 (MQ=255)
                                                                                                                                ggaggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaa              <  2:337002‑M2/150‑114 (MQ=255)
                                                                                                                                   gggACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctc                                                                          >  1:147304‑M2/1‑34 (MQ=255)
                                                                                                                                   gggACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctc                                                                          <  2:147304‑M2/87‑54 (MQ=255)
                                                                                                                                           ccAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACaaa    >  1:909542/1‑150 (MQ=255)
                                                                                                                                            cAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTg                                                             <  2:862283/92‑1 (MQ=255)
                                                                                                                                            cAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGAGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTg                                                             >  1:862283/1‑92 (MQ=255)
                                                                                                                                             aGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACaaaaa  >  1:170271/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAA  >  NC_002947/5988826‑5989116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: