Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,597 C→T 100% P750L (CCG→CTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5970CT100.0% 480.8 / NA 142P750L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (84/58);  total (84/58)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 142 total aligned reads displayed.

CAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGT  >  NC_000913/4187452‑4187742
                                                                                                                                                 |                                                                                                                                                  
cAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGa                                                                                                                                                 >  2:870757/1‑149 (MQ=255)
 aaaCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCGTAACGCAGCATCTACATGATGGCCGTCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGAt                                                                                                                                                <  1:1698747/149‑1 (MQ=255)
 aaaCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGGATGCGTGGTCTGATGGCGAAGCTGGAt                                                                                                                                                >  1:876621/1‑149 (MQ=255)
   aCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATgg                                                                                                                                              >  1:1265055/1‑149 (MQ=255)
     tGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCt                                                                                                                                            >  1:1086955/1‑149 (MQ=255)
        aaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCca                                                                                                                                         >  1:633645/1‑149 (MQ=255)
         aCCGTGATTACCCGGGACGGGCAGGAAGAGAAGCAGGTTTCCTTTGACAGCATCTACATGATGTCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcat                                                                                                                                        <  2:102428/149‑1 (MQ=255)
          ccGTGATTAACCGGGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatc                                                                                                                                       <  2:927605/149‑1 (MQ=255)
                 tAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCGGATGGCGAAGCTGGATGGCCCCATCATCGaaa                                                                                                                                >  2:489470/1‑149 (MQ=255)
                 tAACCGTGACGGGTAGGAAGAGAAGCAGGTTTCCGTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGGCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGaaa                                                                                                                                <  2:1495824/149‑1 (MQ=255)
                   ccaggggcGGTCAGGAAGAGACTTTTGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACg                                                                                                                              <  2:1815834/142‑1 (MQ=255)
                   aCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACg                                                                                                                              <  2:651816/149‑1 (MQ=255)
                    ccGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGc                                                                                                                             >  1:564456/1‑149 (MQ=255)
                     cGTGACGGTCAGGAAGAGAAACAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGGTGGCATAGCTGGGTGGCTCCAGCATCGAAACGcc                                                                                                                            >  2:1272948/1‑149 (MQ=255)
                       tGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACGTGATTGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCaa                                                                                                                          <  1:785415/149‑1 (MQ=255)
                          cGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGACACGCCAATCa                                                                                                                       >  2:1490073/1‑149 (MQ=255)
                           ggTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGGACAGATTCGTCAGCTTGCTGGTATGCGTGGGCTGATGTCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCAc                                                                                                                      >  2:555967/1‑149 (MQ=255)
                           ggTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCCCCATCATCGAAACGCCAATCAc                                                                                                                      >  1:301419/1‑149 (MQ=255)
                             tCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAGGCTGGATGGCTCCATCATCGAAACGCCAATCACcg                                                                                                                    >  2:1545511/1‑149 (MQ=255)
                              cAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCGTGATCGAAACGCCAATCACcgc                                                                                                                   >  1:589587/1‑149 (MQ=255)
                                   agagAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACt                                                                                                              >  2:1423189/1‑149 (MQ=255)
                                       aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTGCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGTTCCATCAGCGGAACGCGAATCACCGCGAACTTcgg                                                                                                          >  2:1404005/1‑147 (MQ=255)
                                        aGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGt                                                                                                         <  1:1564646/149‑1 (MQ=255)
                                        aGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCCCGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGt                                                                                                         <  1:459781/149‑1 (MQ=255)
                                         gCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCGACGCGAACTTCCggt                                                                                                        >  2:542053/1‑147 (MQ=255)
                                         gCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTACGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTg                                                                                                        >  2:1157400/1‑149 (MQ=255)
                                           aGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGaa                                                                                                      <  2:445251/149‑1 (MQ=255)
                                           aGGTTGCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGaa                                                                                                      <  2:1157294/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACCACATGGCCGTGTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                     <  1:1123831/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGCGGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                     >  1:504623/1‑149 (MQ=255)
                                               ttCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTc                                                                                                  >  1:1448699/1‑149 (MQ=255)
                                                  cTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTAGATGGCTCCATCATCGAAACGCCAATCACCGCGGACTTCCGTGACGGTCTGa                                                                                               >  1:1619608/1‑149 (MQ=255)
                                                   ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                              <  1:1395629/149‑1 (MQ=255)
                                                   ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGACAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                              >  2:896978/1‑149 (MQ=255)
                                                    tCAACAGCATCTACGTGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                                                                             <  2:1666272/149‑1 (MQ=255)
                                                    tCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGGCTGCAc                                                                                             >  1:693505/1‑149 (MQ=255)
                                                     cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCCTCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCGGAACg                                                                                            >  2:110310/1‑149 (MQ=255)
                                                     cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGGTACGCGAATCACCGCGAACTTCCGTGAAGGTCCGAACg                                                                                            >  2:375831/1‑149 (MQ=255)
                                                     cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACg                                                                                            >  2:1330379/1‑149 (MQ=255)
                                                      aaCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCCCCGCGGACTTCCGTGAAGGTCCGCACGt                                                                                           >  2:49094/1‑149 (MQ=255)
                                                       aCAGCATCGACAGGATGGCCGCCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTa                                                                                          <  2:1086955/149‑1 (MQ=255)
                                                        caccATCTACCTGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTAc                                                                                         <  2:483450/146‑1 (MQ=255)
                                                         aGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACt                                                                                        >  2:89508/1‑149 (MQ=255)
                                                          gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                       <  1:1172021/149‑1 (MQ=255)
                                                            gtCTACAGGATGGCCGTCGCCGGGGCGCGTGGTTCGGCGGCACAGGTTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCa                                                                                     <  2:1533533/148‑1 (MQ=255)
                                                                aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                                 >  1:368955/1‑149 (MQ=255)
                                                                aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACCGCAGTAc                                                                                 >  2:1225273/1‑149 (MQ=255)
                                                                 cATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTCGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCCTCATCGAAACGCCAATCACCACGAACTGCCGTGAAGGTCTGAACCTACTCAAGTAcg                                                                                >  2:1486818/1‑148 (MQ=255)
                                                                    gatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCa                                                                             >  1:91009/1‑149 (MQ=255)
                                                                      tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGAACTCCAtc                                                                           >  2:811952/1‑149 (MQ=255)
                                                                      tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCGGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGGACTGCCGGGAAGGTCTGAACGTACGCCAGTACTTCAtc                                                                           <  2:835875/149‑1 (MQ=255)
                                                                       ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtct                                                                          <  1:870757/149‑1 (MQ=255)
                                                                       ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtct                                                                          >  1:1751645/1‑149 (MQ=255)
                                                                         ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                                        >  1:323376/1‑149 (MQ=255)
                                                                         ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                                        >  1:1581362/1‑149 (MQ=255)
                                                                          cGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCGTCTcca                                                                       >  1:12758/1‑149 (MQ=255)
                                                                            aCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacc                                                                     >  2:1631648/1‑149 (MQ=255)
                                                                            aCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTGGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacc                                                                     <  1:612946/149‑1 (MQ=255)
                                                                             cttcGGTTCGCGTGGTTCGTCGGTACAGATTCGTTAGCCTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGGCTGAACGTCCTCAAGTACTTCATCTccaccc                                                                    <  2:1165219/146‑1 (MQ=255)
                                                                              cgcGGTGCGCGTGGTTCTGCGGCACAGATGCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGACACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccca                                                                   <  2:204088/147‑1 (MQ=255)
                                                                               ccGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacccac                                                                  <  1:1071360/149‑1 (MQ=255)
                                                                                 ggTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCAATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACCCCAGTACTTCATCTCCACCCACgg                                                                >  2:763446/1‑149 (MQ=255)
                                                                                   tGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAGCTCCACCCACGGTg                                                              >  1:1542572/1‑149 (MQ=255)
                                                                                    gcgcGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGc                                                             <  1:1651766/149‑1 (MQ=255)
                                                                                    gcgcGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGc                                                             <  1:382791/149‑1 (MQ=255)
                                                                                    gcgaGTGGTTCTGCGGCACGGATTCGTCAGCTTGCTGGTATGCGTGGTCTGACGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGc                                                             <  1:1418801/149‑1 (MQ=255)
                                                                                     cgcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTGTGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCt                                                            <  1:797934/149‑1 (MQ=255)
                                                                                      gcgTGGTTCTGCGGCACAGATTCGTCTGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                           <  1:687483/149‑1 (MQ=255)
                                                                                      gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                           >  1:1031346/1‑149 (MQ=255)
                                                                                      gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                           <  1:423896/149‑1 (MQ=255)
                                                                                      gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                           <  1:890684/149‑1 (MQ=255)
                                                                                       tgTGGTTCTGCGGCACAGATTCGTCAGCCTGCGGGTATGCATGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCGAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCg                                                          <  2:338009/148‑1 (MQ=255)
                                                                                            ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                                     >  1:1126648/1‑149 (MQ=255)
                                                                                            ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                                     >  2:1021340/1‑149 (MQ=255)
                                                                                            ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                                     >  2:1440432/1‑149 (MQ=255)
                                                                                              cTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt                                                   <  1:1540188/149‑1 (MQ=255)
                                                                                                   gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGACGGACTGAACGTACTCCAGTACTTCCTCTCCACCCACGGTGCTCGGCAAGGCCGGGc                                              >  2:1356806/1‑149 (MQ=255)
                                                                                                   gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTACATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                              <  1:1157400/149‑1 (MQ=255)
                                                                                                   gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCCACATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                              <  1:606879/149‑1 (MQ=255)
                                                                                                             gTCAGCTTGCTGGTATGCGTGGTCTGATAGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTGCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                   >  1:1181703/1‑149 (MQ=255)
                                                                                                             gTCAGCTGGCTGGTATGCGGGGTCGGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                   <  2:219337/149‑1 (MQ=255)
                                                                                                              tCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGACACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCAc                                  >  2:589565/1‑149 (MQ=255)
                                                                                                              tCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGccc                                  >  1:466434/1‑147 (MQ=255)
                                                                                                                  cTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCGTCATCGAAACGCCAATGCCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTCCTTGATCCGCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCACCCTGaa                              >  1:1373292/1‑149 (MQ=255)
                                                                                                                    tGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaaaa                            >  1:246975/1‑149 (MQ=255)
                                                                                                                          tATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACGTCCGCGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGa                      <  1:867652/149‑1 (MQ=255)
                                                                                                                            tGCGTGGTCTGATGGCGAAGCTGGATGGCGCCATCAGCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc                    >  2:616836/1‑149 (MQ=255)
                                                                                                                             gCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACt                   >  1:39414/1‑149 (MQ=255)
                                                                                                                             gCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATAGAAACGCCAATCACCGCGAACTTCCGGGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGGCTGGCGGGTACCGCACTGACACCTGCGACCt                   >  2:767391/1‑149 (MQ=255)
                                                                                                                              cGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAGTCGCCGCGAGCTTCCGTGAAGGTCTGAAAGTACTCCAGTACTTCATCTCCACCCACGGTGCTGGTCAAGGACTGGCGAATAACGCACTGAACAGTGCGacatc                  >  2:807008/1‑145 (MQ=255)
                                                                                                                               gTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGCTACCGCACTGACAACTGCGAACTcc                 >  1:1674544/1‑149 (MQ=255)
                                                                                                                                tGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCg                <  1:1365643/149‑1 (MQ=255)
                                                                                                                                tGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCCCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACCCCg                >  2:1378894/1‑149 (MQ=255)
                                                                                                                                  gTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGt              >  2:861361/1‑149 (MQ=255)
                                                                                                                                   tCTGATGGCGCAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGACACCGCACTGAAAACTGCGAACTCCGGtt             <  2:914618/149‑1 (MQ=255)
                                                                                                                                      gATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAcc          >  1:898007/1‑149 (MQ=255)
                                                                                                                                       aTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGCGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCCCGGTGCTCGTAAAGGGCTGACGGAGCCCAGCATGACAACTTGGAACGCCGGTTACCt         >  2:123179/1‑149 (MQ=255)
                                                                                                                                          gCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGAc      >  2:1007614/1‑149 (MQ=255)
                                                                                                                                            gAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAGCTCCGGTTACCTGACtc    >  2:402668/1‑149 (MQ=255)
                                                                                                                                              aGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcgt  >  1:936913/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGT  >  NC_000913/4187452‑4187742

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