Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,391,407 A→C 12.5% F208V (TTT→GTT)  tldD ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,391,4070AC12.5% 1715.1 / 90.8 655F208V (TTT→GTT) tldDputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (401/159);  new base C (2/93);  total (403/252)
Fisher's exact test for biased strand distribution p-value = 1.89e-41
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 684 total aligned reads displayed.

CGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGACACCACTGAGGC  >  NC_000913/3391258‑3391556
                                                                                                                                                     |                                                                                                                                                     
 ggCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCaaaa                                                                                                                                                      >  1:8003525/1‑150 (MQ=255)
       ccTGCGCATAGTACTACCGGCTTGGTGCCCGCTGGTGAAGAAACGGCATAAAGTTTGACCCGCGCCATACGCACTGTTTCTTTTGCCCTTGCATCCGCACGGACTTCGCCGTCTTTATGGGAATGGATTATTTCATAACAAAAACGAccg                                                                                                                                                <  2:7629129/150‑1 (MQ=255)
        cagcgcCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgc                                                                                                                                               >  1:171212/3‑150 (MQ=255)
         tGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgcc                                                                                                                                              >  2:3294314/1‑150 (MQ=255)
          gcgcCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgccg                                                                                                                                             >  1:3538237/1‑150 (MQ=255)
          gcaccAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgccg                                                                                                                                             >  1:1104305/4‑150 (MQ=255)
              cAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACcgccgccgc                                                                                                                                         >  1:3299518/1‑150 (MQ=255)
              cAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACcgccgccgc                                                                                                                                         >  1:2736368/1‑150 (MQ=255)
                  tactacCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACt                                                                                                                                     >  2:6775717/1‑150 (MQ=255)
                     tacCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGc                                                                                                                                  >  1:8573790/1‑150 (MQ=255)
                        cGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCAcc                                                                                                                               >  2:4027788/1‑150 (MQ=255)
                           cATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcg                                                                                                                            >  2:2054498/1‑150 (MQ=255)
                              ggTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttc                                                                                                                         >  1:7330394/1‑150 (MQ=255)
                               gTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGACTTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcg                                                                                                                        >  2:7022352/1‑150 (MQ=255)
                               gTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcg                                                                                                                        >  2:1231349/1‑150 (MQ=255)
                               gTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcg                                                                                                                        >  1:3056316/1‑150 (MQ=255)
                               gTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcg                                                                                                                        >  2:7827833/1‑150 (MQ=255)
                                 gCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcgcg                                                                                                                      >  1:999055/1‑150 (MQ=255)
                                    cGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGttt                                                                                                                   >  2:5147179/1‑150 (MQ=255)
                                    cGCTGGTGAATCAACGGCTCAAATATTGACCAGCCCAATCCGCACTGCTTATTTTGCCCATGCATCCGAACGGACTTCGCCGTCGTGATCGGAAAGTAACAATTCATAACAAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGttt                                                                                                                   <  2:8030558/150‑1 (MQ=255)
                                            cagcaACGGCATAAAGTTGGACCAGCGCCATCCGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCACGGAAGATTTCATAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           <  2:4466052/150‑1 (MQ=255)
                                            cagcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           >  1:4019724/1‑150 (MQ=255)
                                            cagcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           >  1:6591441/1‑150 (MQ=255)
                                            cagcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           >  1:5815006/1‑150 (MQ=255)
                                            cagcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGAAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           <  2:1968187/150‑1 (MQ=255)
                                            cagcaACGGCAGAAAGATTCACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                                                                                                           >  1:7552462/1‑150 (MQ=255)
                                               ccaCGGCAGAACGATGGACCCGCGCCAGCCGCACTGCTTCTTTTGCCCAGGCATCCGCACGGACGTCGCCGTCGAGATCGGCAAGGAAGATTTCATCACCACAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATcttctt                                                                                                        <  2:6428396/148‑1 (MQ=255)
                                               caACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATcttctt                                                                                                        >  2:3173647/1‑150 (MQ=255)
                                                aaCGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCttcttc                                                                                                       >  1:2093594/1‑150 (MQ=255)
                                                  cGGCATAAAGATTGACCAGCGCCATCCGCACTGCTTCTTTTGCCCAGGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGCAATAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTcga                                                                                                     <  2:2647131/150‑1 (MQ=255)
                                                   ggCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTcgac                                                                                                    >  2:5525448/1‑150 (MQ=255)
                                                     cAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTcgacga                                                                                                  >  2:5438438/1‑150 (MQ=255)
                                                      aGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACgag                                                                                                 >  2:363054/1‑150 (MQ=255)
                                                        aaaGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGaa                                                                                               >  2:2403388/1‑150 (MQ=255)
                                                        aaaGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGaa                                                                                               >  1:7838629/1‑150 (MQ=255)
                                                            aTTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCt                                                                                           >  1:4412622/1‑150 (MQ=255)
                                                               gACCAGCGCCATACGCACTGCTTCTTTTGCCCAGGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCAc                                                                                        <  2:4635722/150‑1 (MQ=255)
                                                                aCCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACg                                                                                       >  1:5942666/1‑150 (MQ=255)
                                                                aCCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACg                                                                                       >  2:2964651/1‑150 (MQ=255)
                                                                aCCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACg                                                                                       >  2:2317979/1‑150 (MQ=255)
                                                                  cAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCaaaa                                                                                                                                                      >  1:7234566/1‑85 (MQ=255)
                                                                  cAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCaaaa                                                                                                                                                      <  2:7234566/85‑1 (MQ=255)
                                                                   cgcgcAATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCATGGAATAATTCATAACCACAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGaa                                                                                    <  1:5116513/149‑1 (MQ=255)
                                                                   aGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGaa                                                                                    >  1:1603852/1‑150 (MQ=255)
                                                                      cccATCCGCACTGCTTCTTTTGCCCAGGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAATAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGa                                                                                 <  1:7027066/149‑1 (MQ=255)
                                                                      cccATACGCACTGCTTCTTTTGCCCATGAACCCGCACGTAGTTCGCCGTCGAGTTCGGCATGTAATAATTAATAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGa                                                                                 <  1:6947389/149‑1 (MQ=255)
                                                                       ccATCCGCAGTGCTTCTTTTGCCAATGCTTCCGCAGGGACTTCGCCGTCGAGAGCGGCATGGAATAATTCATAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGAc                                                                                <  2:7510771/150‑1 (MQ=255)
                                                                           aCGCACTGCTTCTTTTGCCCATGCATCCGCACGTACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCAc                                                                            <  2:171212/150‑1 (MQ=255)
                                                                             gCACGGCTTCTTTGGCCCCTCCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCa                                                                          <  2:5885803/150‑1 (MQ=255)
                                                                              cACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAg                                                                         >  1:6195471/1‑150 (MQ=255)
                                                                                 tGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGcgg                                                                      >  1:831039/1‑150 (MQ=255)
                                                                                   cttcttTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGACTTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGcggac                                                                    <  2:4871437/150‑1 (MQ=255)
                                                                                      cttTTGCCCTTGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGTAATAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGcggacgga                                                                 <  1:7733618/150‑1 (MQ=255)
                                                                                          tGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGTAATAATTCAAAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATc                                                             <  1:2332375/150‑1 (MQ=255)
                                                                                           gCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATcc                                                            >  2:769437/1‑150 (MQ=255)
                                                                                            cccATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATccg                                                           <  2:8552257/150‑1 (MQ=255)
                                                                                             ccATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATccgc                                                          >  1:4271129/1‑150 (MQ=255)
                                                                                             ccATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATccgc                                                          >  2:605448/1‑150 (MQ=255)
                                                                                               aTGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATccgccg                                                        >  1:3688229/1‑150 (MQ=255)
                                                                                                  cATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTa                                                     <  1:5614427/150‑1 (MQ=255)
                                                                                                    tCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGc                                                   >  2:2888213/1‑150 (MQ=255)
                                                                                                      cGCACGGACTTCGCCGTCGAGATCGGCACGGAAGACTGCATCACCACCCCGACCGCCGCCGCCACGGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGt                                                 <  2:4386510/150‑1 (MQ=255)
                                                                                                       gcccGGCCTGCGCCGTCGAGATCGTCAAGTAAGAATACAAACCCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTg                                                <  1:6653201/147‑1 (MQ=255)
                                                                                                        cACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGc                                               >  1:3837909/1‑150 (MQ=255)
                                                                                                         aCGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGcc                                              >  1:547269/1‑150 (MQ=255)
                                                                                                          cGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCg                                             <  2:7101306/150‑1 (MQ=255)
                                                                                                               ttCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGt                                        >  1:257313/1‑150 (MQ=255)
                                                                                                               tgcGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGt                                        <  1:4152231/148‑1 (MQ=255)
                                                                                                                tCGCCGTCGAGATCGGCAAGGAATAATTAATACCCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTg                                       <  1:7404311/150‑1 (MQ=255)
                                                                                                                 cGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTgg                                      <  2:1190698/150‑1 (MQ=255)
                                                                                                                   ccGTCGAGATCGGCAAGGAAGAATTCATAACCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGcc                                    <  1:8657388/150‑1 (MQ=255)
                                                                                                                   ccGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGcc                                    <  2:6556286/150‑1 (MQ=255)
                                                                                                                    cGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCg                                   >  2:6225169/1‑150 (MQ=255)
                                                                                                                    cGTCGAGATCGGAAAGTAAGAATTCATAACCAAAACGACCGCTGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCg                                   <  1:5306241/150‑1 (MQ=255)
                                                                                                                       cGAGATCGGCACGGAAGAATTCATACCCACCACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCaa                                <  1:1013336/150‑1 (MQ=255)
                                                                                                                       cGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCaa                                <  1:7614904/150‑1 (MQ=255)
                                                                                                                       cGAGATCGGCAAGGAAGAATACATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCaa                                >  1:6383578/1‑150 (MQ=255)
                                                                                                                       cGAGATCGGCAAGGAAGAATACATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCaa                                <  2:3429490/150‑1 (MQ=255)
                                                                                                                          gaTCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCa                             >  1:298360/1‑150 (MQ=255)
                                                                                                                          gaTCGGAAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCa                             <  2:6852360/150‑1 (MQ=255)
                                                                                                                             cgtcatggAAGAATTAATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCaaaa                          <  1:4532084/144‑1 (MQ=255)
                                                                                                                              ggCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaat                         >  2:5645886/1‑150 (MQ=255)
                                                                                                                              ggCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaat                         >  2:3835715/1‑150 (MQ=255)
                                                                                                                              ggCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaat                         >  2:8142728/1‑150 (MQ=255)
                                                                                                                               gCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaatt                        >  1:1241749/1‑150 (MQ=255)
                                                                                                                                 aaGGAAGATTTCATACCCAACACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaattaa                      <  2:3983542/150‑1 (MQ=255)
                                                                                                                                 aaGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaattaa                      <  1:8472902/150‑1 (MQ=255)
                                                                                                                                  aGGAAGAATTCATAACCAAAACGACCGCCGGCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACGAGCGGACGGAGATCCGGCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaattaat                     <  2:6666659/150‑1 (MQ=255)
                                                                                                                                  aGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAaattaat                     <  1:8221043/150‑1 (MQ=255)
                                                                                                                                    gaagaaTTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTc                   >  2:3015842/1‑150 (MQ=255)
                                                                                                                                        aaTTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAg               >  1:4887967/1‑150 (MQ=255)
                                                                                                                                         aTTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGa              >  1:346208/1‑150 (MQ=255)
                                                                                                                                           tCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGTCGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGaca            <  2:5753351/150‑1 (MQ=255)
                                                                                                                                           tCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGaca            <  1:677972/150‑1 (MQ=255)
                                                                                                                                            cATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGacac           >  2:139753/1‑150 (MQ=255)
                                                                                                                                            aaTAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGacac           >  2:8365200/2‑150 (MQ=255)
                                                                                                                                              tAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGAcacca         >  2:150204/1‑150 (MQ=255)
                                                                                                                                                  cAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGACACCACTGa     <  2:5213871/150‑1 (MQ=255)
                                                                                                                                                    aaaCGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCACCCAAAATTAATTCATATACACCCCTGAgg   >  1:1976342/1‑150 (MQ=255)
                                                                                                                                                    aaaCGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGACACCACTGAgg   >  2:2162492/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACGCTCACGGAAAGACGCACCAGCGGACGGACATCCGCCGCTAGCGTGCCGTCGGTGGCCGCAACCAAAATTAATTCATAGACACCACTGAGGC  >  NC_000913/3391258‑3391556

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