Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,893,2710CT31.6% 111.1 / 56.7 79A169A (GCC→GCTyidZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base T (9/16);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGCAGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACA  >  NC_000913/3893133‑3893399
                                                                                                                                          |                                                                                                                                
cccGCAGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGcc                                                                                                                                  <  1:1025855/139‑1 (MQ=255)
cccGCAGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGcc                                                                                                                                  <  1:1105120/139‑1 (MQ=255)
  cGCAGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCAt                                                                                                                                <  2:91022/139‑1 (MQ=255)
     aGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGa                                                                                                                             <  2:774580/139‑1 (MQ=255)
     aGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGa                                                                                                                             >  1:1157690/1‑139 (MQ=255)
          aCCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATg                                                                                                                        >  1:1247299/1‑139 (MQ=255)
            cATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGaa                                                                                                                      <  1:1013198/139‑1 (MQ=255)
             aTCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAg                                                                                                                     <  2:1170217/139‑1 (MQ=255)
              tCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGt                                                                                                                    <  2:1157690/139‑1 (MQ=255)
               cAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTg                                                                                                                   <  2:962198/139‑1 (MQ=255)
               cAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTg                                                                                                                   <  1:544703/139‑1 (MQ=255)
               cAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   <  2:1031858/139‑1 (MQ=255)
               cAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   <  2:1073915/139‑1 (MQ=255)
               cAAAGTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTg                                                                                                                   <  2:951269/139‑1 (MQ=255)
                               gATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGc                                                                                                   >  1:786506/1‑139 (MQ=255)
                                aTTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGcc                                                                                                  <  2:1247299/139‑1 (MQ=255)
                                  tACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCt                                                                                                >  2:324317/1‑139 (MQ=255)
                                       ttCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTc                                                                                           <  1:324317/139‑1 (MQ=255)
                                        tCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCt                                                                                          >  2:291323/1‑139 (MQ=255)
                                         ccTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   <  2:712316/113‑1 (MQ=255)
                                         ccTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   >  1:712316/1‑113 (MQ=255)
                                         ccTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCt                                                                                                >  2:66346/1‑132 (MQ=255)
                                         ccTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCt                                                                                                <  1:66346/132‑1 (MQ=255)
                                          cTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTgg                                                                                        <  2:727263/139‑1 (MQ=255)
                                          cTTAGATGCCATTACTCGTGGTGAAGCGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTgg                                                                                        >  1:138363/1‑139 (MQ=255)
                                           ttAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       <  1:929781/139‑1 (MQ=255)
                                              gATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTAc                                                                                    <  1:931200/139‑1 (MQ=255)
                                               aTGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACg                                                                                   <  1:417210/139‑1 (MQ=255)
                                                tGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGc                                                                                  <  1:838790/139‑1 (MQ=255)
                                                   cATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGaa                                                                                                                      <  2:329187/100‑1 (MQ=255)
                                                   cATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGaa                                                                                                                      >  1:329187/1‑100 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAg                                                                              >  1:874142/1‑139 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAg                                                                              >  2:1106219/1‑139 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTa                                                                                                            <  1:105219/109‑1 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTa                                                                                                            >  2:105219/1‑109 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       <  1:666897/130‑1 (MQ=255)
                                                    aTTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       >  2:666897/1‑130 (MQ=255)
                                                     ttACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTg                                                                                                                   <  2:1278003/101‑1 (MQ=255)
                                                     ttACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTg                                                                                                                   >  1:1278003/1‑101 (MQ=255)
                                                     ttACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGa                                                                             <  1:291323/139‑1 (MQ=255)
                                                         tCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGatcat                                                                         >  1:899900/1‑139 (MQ=255)
                                                          cGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGatcatc                                                                        <  2:801142/139‑1 (MQ=255)
                                                            tggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAg                                                                                                                     <  2:283119/92‑1 (MQ=255)
                                                            tggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAg                                                                                                                     >  1:283119/1‑92 (MQ=255)
                                                             ggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCAt                                                                                                                                <  2:622976/80‑1 (MQ=255)
                                                             ggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCAt                                                                                                                                >  1:622976/1‑80 (MQ=255)
                                                             ggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       >  1:1292970/1‑121 (MQ=255)
                                                             ggtgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       <  2:1292970/121‑1 (MQ=255)
                                                               tgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCa                                                                   <  1:1106219/139‑1 (MQ=255)
                                                               tgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCa                                                                   <  2:138363/139‑1 (MQ=255)
                                                               tgAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCa                                                                   <  1:238329/139‑1 (MQ=255)
                                                                 aaGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   >  1:572303/1‑89 (MQ=255)
                                                                 aaGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTg                                                                                                                   <  2:572303/89‑1 (MQ=255)
                                                                   gTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGc                                                               >  1:97709/1‑139 (MQ=255)
                                                                           cGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATg                                                                                                                        >  1:509737/1‑74 (MQ=255)
                                                                           cGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCTATTGATTATg                                                                                                                        <  2:509737/74‑1 (MQ=255)
                                                                           cGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACg                                                       >  1:607040/1‑139 (MQ=255)
                                                                             gTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATg                                                                                                                        >  2:1107252/1‑72 (MQ=255)
                                                                             gTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATg                                                                                                                        <  1:1107252/72‑1 (MQ=255)
                                                                              tttttCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCAc                                                                  <  1:1208986/125‑1 (MQ=255)
                                                                              tttttCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCAc                                                                  >  2:1208986/1‑125 (MQ=255)
                                                                                      gTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCAt                                                                                                                                <  2:151834/55‑1 (MQ=255)
                                                                                      gTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCAt                                                                                                                                >  1:151834/1‑55 (MQ=255)
                                                                                      gTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGa                                            >  1:976566/1‑139 (MQ=255)
                                                                                        cgcgAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGaca                                          <  1:505025/139‑1 (MQ=255)
                                                                                        cgcgAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGaca                                          <  2:906114/139‑1 (MQ=255)
                                                                                                          tcgcGGGAGCTGTTAAGCTCTCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGa                                            >  2:772320/1‑119 (MQ=255)
                                                                                                          tcgcGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGa                                            <  1:772320/119‑1 (MQ=255)
                                                                                                              gggAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCatat                    >  2:71253/1‑139 (MQ=255)
                                                                                                                    tGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCAt              <  2:976566/139‑1 (MQ=255)
                                                                                                                    tGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCAt              <  2:1050051/139‑1 (MQ=255)
                                                                                                                    tGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCAt              <  2:899900/139‑1 (MQ=255)
                                                                                                                      ttAAGCTCGCTACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCAttt            >  2:1253657/1‑139 (MQ=255)
                                                                                                                      ttAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       <  1:1113101/64‑1 (MQ=255)
                                                                                                                      ttAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGt                                                                                       >  2:1113101/1‑64 (MQ=255)
                                                                                                                              gCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGaa    <  2:749689/139‑1 (MQ=255)
                                                                                                                                tACCGTTAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACa  <  2:483013/139‑1 (MQ=255)
                                                                                                                                     ttAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGaca                                          <  2:772120/94‑1 (MQ=255)
                                                                                                                                     ttAGCTATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGaca                                          >  1:772120/1‑94 (MQ=255)
                                                                                                                                        gCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACctt                                      <  2:375623/95‑1 (MQ=255)
                                                                                                                                        gCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACctt                                      >  1:375623/1‑95 (MQ=255)
                                                                                                                                          |                                                                                                                                
CCCGCAGGCGACCATCAAATTACGTAACTGGGATTACGATTCCTTAGATGCCATTACTCGTGGTGAAGTGGATATCGGTTTTTCCGGTCGCGAAAGCCATCCTCGCTCGCGGGAGCTGTTAAGCTCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACA  >  NC_000913/3893133‑3893399

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