Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,185,708 T→C L120P (CTG→CCG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,7080TC98.7% 267.9 / ‑4.1 77L120P (CTG→CCG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (24/52);  total (25/52)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

CTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCT  >  NC_000913/4185560‑4185855
                                                                                                                                                    |                                                                                                                                                   
cTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATcgccgc                                                                                                                                                   >  2:541419/1‑151 (MQ=255)
  gtgCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCgccgcc                                                                                                                                                  <  1:1051651/150‑1 (MQ=255)
    gCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTc                                                                                                                                               <  2:734426/151‑1 (MQ=255)
      ggTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTccc                                                                                                                                             <  2:1181765/151‑1 (MQ=255)
                aGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTct                                                                                                                                   <  2:151160/151‑1 (MQ=255)
                  cgcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgc                                                                                                                                 <  2:174559/151‑1 (MQ=255)
                   gcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgct                                                                                                                                <  2:776745/151‑1 (MQ=255)
                       gAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTc                                                                                                                            <  2:1389091/151‑1 (MQ=255)
                         aaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTgctgct                                                                                                                             <  2:1344408/148‑1 (MQ=255)
                         aaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGa                                                                                                                          <  2:1170017/151‑1 (MQ=255)
                          acacCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGat                                                                                                                         <  1:258326/151‑1 (MQ=255)
                               gTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGcc                                                                                                                    <  2:356081/151‑1 (MQ=255)
                                  gCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCg                                                                                                                   <  2:180836/149‑1 (MQ=255)
                                     tCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcg                                                                                                              >  1:147282/1‑151 (MQ=255)
                                      cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcgc                                                                                                             >  2:1064247/1‑151 (MQ=255)
                                        tCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGcgcg                                                                                                            >  1:1260734/1‑150 (MQ=255)
                                           gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGatat                                                                                                        <  1:689568/151‑1 (MQ=255)
                                            tgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATc                                                                                                       >  2:1078964/1‑151 (MQ=255)
                                             gagaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATc                                                                                                       <  2:831756/150‑1 (MQ=255)
                                              agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAGCTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGa                                                                                                     <  2:935840/151‑1 (MQ=255)
                                              agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCg                                                                                                      >  2:981540/1‑150 (MQ=255)
                                              agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGa                                                                                                     <  1:620824/151‑1 (MQ=255)
                                               gaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGaa                                                                                                    <  1:400762/151‑1 (MQ=255)
                                                aaGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAc                                                                                                   <  1:529436/151‑1 (MQ=255)
                                                aaGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACACCGAACTGGCTTCCCCGACTGCGCAAATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAc                                                                                                   >  1:1125731/1‑151 (MQ=255)
                                                     cggcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATGGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgt                                                                                               >  1:1199880/1‑150 (MQ=255)
                                                      ggcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtac                                                                                             >  2:1006669/1‑151 (MQ=255)
                                                       gcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtact                                                                                            <  1:1311827/151‑1 (MQ=255)
                                                        cgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtactg                                                                                           <  2:527165/151‑1 (MQ=255)
                                                         gTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtactgt                                                                                          <  2:140111/151‑1 (MQ=255)
                                                           tGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtactgtac                                                                                        >  2:324012/1‑151 (MQ=255)
                                                           tGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtactgtac                                                                                        >  2:168453/1‑151 (MQ=255)
                                                             aaGTGACCCAGGCTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACtt                                                                                      <  1:196478/151‑1 (MQ=255)
                                                             aaGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACtt                                                                                      <  2:246401/151‑1 (MQ=255)
                                                             aaGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCAAAAGCGTACTCTAAtt                                                                                      >  2:774759/1‑151 (MQ=255)
                                                              aGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACttt                                                                                     <  2:259038/151‑1 (MQ=255)
                                                               gTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTg                                                                                    <  1:383365/151‑1 (MQ=255)
                                                                        aCTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTAtg                                                                           <  2:917623/151‑1 (MQ=255)
                                                                             aGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGtt                                                                      <  2:1284189/151‑1 (MQ=255)
                                                                                aCGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTAt                                                                    <  1:1064822/150‑1 (MQ=255)
                                                                                aCGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATc                                                                   >  1:1349915/1‑151 (MQ=255)
                                                                                aCGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATc                                                                   <  1:1166668/151‑1 (MQ=255)
                                                                                 cGCCGGGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCg                                                                  <  2:161516/151‑1 (MQ=255)
                                                                                  gCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCg                                                                  >  2:514566/1‑150 (MQ=255)
                                                                                    cGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGa                                                                 >  2:543149/1‑149 (MQ=255)
                                                                                    cGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAg                                                               <  2:182551/151‑1 (MQ=255)
                                                                                       gAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCg                                                            >  2:766722/1‑151 (MQ=255)
                                                                                       gAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCg                                                            <  2:658266/151‑1 (MQ=255)
                                                                                       gAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCCGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCg                                                            <  2:954087/151‑1 (MQ=255)
                                                                                        aGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCgg                                                           >  2:350728/1‑151 (MQ=255)
                                                                                        aGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCgg                                                           <  1:841443/151‑1 (MQ=255)
                                                                                         gCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCgg                                                           <  1:1423866/150‑1 (MQ=255)
                                                                                         gCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGt                                                          <  2:47250/151‑1 (MQ=255)
                                                                                          cGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTa                                                         <  1:49231/151‑1 (MQ=255)
                                                                                            tATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATg                                                       >  2:207645/1‑151 (MQ=255)
                                                                                             aTGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGa                                                      <  1:1269433/151‑1 (MQ=255)
                                                                                               gggCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGAc                                                     <  2:888261/150‑1 (MQ=255)
                                                                                                ggCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCa                                                   <  2:115962/151‑1 (MQ=255)
                                                                                                 gCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCa                                                   >  1:439459/1‑150 (MQ=255)
                                                                                                   cacaTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAAcc                                                >  1:108811/1‑151 (MQ=255)
                                                                                                    acaTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCt                                               <  2:1054051/151‑1 (MQ=255)
                                                                                                         gAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAAc                                          <  2:136100/151‑1 (MQ=255)
                                                                                                           aCTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGt                                        <  1:1296027/151‑1 (MQ=255)
                                                                                                              ggCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTca                                      <  2:395142/150‑1 (MQ=255)
                                                                                                               gCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                    <  1:1080789/151‑1 (MQ=255)
                                                                                                                cTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                    >  1:1121694/1‑150 (MQ=255)
                                                                                                                   ccccGTCTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGa                                 <  1:501832/150‑1 (MQ=255)
                                                                                                                   ccccGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                    >  2:121798/1‑147 (MQ=255)
                                                                                                                      cGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCct                             <  2:553412/151‑1 (MQ=255)
                                                                                                                           gcgcACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATTCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactg                        <  1:852219/151‑1 (MQ=255)
                                                                                                                            cgcACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactga                       <  1:244769/151‑1 (MQ=255)
                                                                                                                               acaTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactga                       >  2:1417500/1‑148 (MQ=255)
                                                                                                                               acaTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAaga                    <  2:138995/151‑1 (MQ=255)
                                                                                                                                caTCTGGTTCCTGAAATCGCCGCTGTCCCGTATCGGTCTGCTGCTCGAGATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAagag                   <  1:1330487/151‑1 (MQ=255)
                                                                                                                                 aTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGc                  >  2:1168816/1‑151 (MQ=255)
                                                                                                                                 aTCTGGTTCCTGAAAGCGCCGCCGTCCCGTATCGGTTTGCTGCGCGATATGGCGCTGCGCGATATCGATGGCGTACTGTACTTTGAAGCCTATGTGGTTATCGAAGGCGGTATGTCCAGCCTGGAACGTCAGCAGAGCCTGTCTGAagag                   <  1:1114788/150‑1 (MQ=255)
                                                                                                                                        tCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCt           >  1:696349/1‑151 (MQ=255)
                                                                                                                                              aaTCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAcg      >  1:1364857/1‑150 (MQ=255)
                                                                                                                                               aTCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAcgcg    <  2:900556/151‑1 (MQ=255)
                                                                                                                                                 cgccgcCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCt  >  1:365589/1‑151 (MQ=255)
                                                                                                                                                 cgccgcCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCt  <  2:729230/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCT  >  NC_000913/4185560‑4185855

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