Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,051,413 T→C K369E (AAG→GAG)  PA14_RS13990 ← SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,051,4130TC99.8% 1968.4 / NA 500K369E (AAG→GAG) PA14_RS13990SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (284/215);  minor base A (0/1);  total (284/216)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 500 total aligned reads displayed.

GCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051266‑3051559
                                                                                                                                                   |                                                                                                                                                  
 cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                                 <  2:5930679/150‑1 (MQ=255)
  tCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCg                                                                                                                                                <  1:8080661/150‑1 (MQ=255)
   cATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGa                                                                                                                                               >  1:7052105/1‑150 (MQ=255)
     tCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGc                                                                                                                                             <  1:10208545/150‑1 (MQ=255)
     tCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGc                                                                                                                                             >  2:9775381/1‑150 (MQ=255)
        ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGCACTCGGCGAGCGTg                                                                                                                                          >  1:1587838/1‑150 (MQ=255)
           gcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcg                                                                                                                                       >  1:3399682/1‑150 (MQ=255)
            ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      >  2:3079040/1‑150 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTCGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  1:1225879/150‑1 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  2:10521130/150‑1 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  2:6815123/150‑1 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  2:4020098/150‑1 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  1:5364818/150‑1 (MQ=255)
                 cAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTGGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTg                                                                                                                                 >  1:3685545/1‑150 (MQ=255)
                   gATCCATGAGTGACGCGCCGGACGCGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGcc                                                                                                                               <  1:5286403/150‑1 (MQ=255)
                     tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             >  2:3363303/1‑150 (MQ=255)
                     tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             >  1:5888240/1‑150 (MQ=255)
                     tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             >  1:6884375/1‑150 (MQ=255)
                            gTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtc                                                                                                                      <  1:10153831/150‑2 (MQ=255)
                             tGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtcg                                                                                                                     >  1:1932821/1‑148 (MQ=255)
                                cgcgCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAag                                                                                                                  <  2:11109346/150‑1 (MQ=255)
                                cgcgCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAag                                                                                                                  >  2:7664895/1‑150 (MQ=255)
                                  cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                                >  2:2567422/1‑150 (MQ=255)
                                  cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                                >  1:7872367/1‑150 (MQ=255)
                                   gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                               >  1:2515867/1‑150 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              <  1:10428386/150‑1 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              >  2:1588724/1‑150 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              >  1:4207059/1‑150 (MQ=255)
                                     cGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCtt                                                                                                             >  1:6773375/1‑150 (MQ=255)
                                       gACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtct                                                                                                           >  2:7914992/1‑150 (MQ=255)
                                          gggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGc                                                                                                        <  1:8859525/150‑1 (MQ=255)
                                           ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                                       >  1:3254786/1‑150 (MQ=255)
                                           ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                                       >  1:7098733/1‑150 (MQ=255)
                                             cgcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTgg                                                                                                     >  2:270328/1‑150 (MQ=255)
                                                cGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcgg                                                                                                  >  1:8422892/1‑150 (MQ=255)
                                                  ttCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                >  1:10588026/1‑150 (MQ=255)
                                                  ttCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                >  1:2965628/1‑150 (MQ=255)
                                                  ttCAGGTAGTAGTCGCCGACCGCGAGGCCCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                <  2:4144422/150‑1 (MQ=25)
                                                       gtagtaGGCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCt                                                                                           <  2:9224688/150‑1 (MQ=255)
                                                         agtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGa                                                                                         <  2:5082227/150‑1 (MQ=255)
                                                           tagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtc                                                                                       >  2:5806757/1‑150 (MQ=255)
                                                            agtCGCCGACCGCCTGGACCTTCCAGCGCACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtct                                                                                      <  2:6969827/150‑1 (MQ=25)
                                                             gtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGCGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCGGAtctc                                                                                     <  2:10353176/150‑1 (MQ=38)
                                                              tCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCg                                                                                    >  1:9287228/1‑150 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           >  2:3840508/1‑150 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           >  2:421049/1‑150 (MQ=255)
                                                                       gcgTGGACCTGCCAGCGTACCGGGTCGAGCAGGGTGTGGGAGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGACCTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCacc                                                                           <  1:6052678/150‑3 (MQ=255)
                                                                           ggACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTACCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCt                                                                       <  1:991486/150‑1 (MQ=255)
                                                                             aCCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTg                                                                     >  2:6329085/1‑150 (MQ=255)
                                                                             aCCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTg                                                                     <  2:7194899/150‑1 (MQ=255)
                                                                                ttCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGcc                                                                  <  2:223430/150‑1 (MQ=255)
                                                                                ttCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGcc                                                                  <  2:3644099/150‑1 (MQ=255)
                                                                                ttCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGcc                                                                  <  1:4422061/150‑1 (MQ=255)
                                                                                  ccAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTc                                                                <  1:9756623/150‑1 (MQ=255)
                                                                                   cAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCg                                                               <  1:6553298/150‑1 (MQ=255)
                                                                                    aGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCgg                                                              >  1:9524074/1‑150 (MQ=255)
                                                                                       gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCa                                                           >  1:7314739/1‑150 (MQ=255)
                                                                                       gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCa                                                           <  1:2369975/150‑1 (MQ=255)
                                                                                       gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCa                                                           >  2:7578606/1‑150 (MQ=255)
                                                                                        tACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGGGCCAGTGCCGGTCGAGGTTGAAGTCGGCGTGCGTGGCGCGGCTGCCGGCCTGCTCGATGTGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGGACCTTGGGCTCGGcccg                                                          >  2:10994620/1‑148 (MQ=37)
                                                                                         aCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACg                                                         <  1:11000233/150‑1 (MQ=255)
                                                                                          ccGGGTCGAGCAGGGTGTGTGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTAGGCGAGCGGGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTGGGCCTCGGCCACgg                                                        <  2:10870126/150‑1 (MQ=255)
                                                                                             ggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGa                                                     <  1:4868262/150‑1 (MQ=255)
                                                                                             ggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGa                                                     >  1:3962842/1‑150 (MQ=255)
                                                                                              gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGCACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  1:5029063/150‑1 (MQ=255)
                                                                                              gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCGGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  2:4396602/150‑1 (MQ=255)
                                                                                               tCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATc                                                   >  1:5506749/1‑150 (MQ=255)
                                                                                                 gTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGa                                                 >  1:7693305/1‑150 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGTGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                >  2:138220/1‑150 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                >  2:3149683/1‑150 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                >  2:8741410/1‑150 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                >  1:4144912/1‑150 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGGGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                <  1:6308036/150‑1 (MQ=255)
                                                                                                    cAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgc                                              >  1:2777097/1‑150 (MQ=255)
                                                                                                     aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcg                                             <  2:8817313/150‑1 (MQ=255)
                                                                                                         tgtgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTg                                         >  2:5294547/1‑150 (MQ=255)
                                                                                                          gtgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgc                                        >  2:4711516/1‑150 (MQ=255)
                                                                                                           tgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       >  2:1965321/1‑150 (MQ=255)
                                                                                                           tgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       <  2:7960234/150‑1 (MQ=255)
                                                                                                           ggggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       <  1:6018273/149‑1 (MQ=255)
                                                                                                             ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                                     <  2:4966327/150‑1 (MQ=255)
                                                                                                                gCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGa                                  <  1:569584/150‑1 (MQ=255)
                                                                                                                  gggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACt                                >  2:5379874/1‑150 (MQ=255)
                                                                                                                  gggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACt                                <  1:4605493/150‑1 (MQ=255)
                                                                                                                   ggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCACAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTc                               <  1:5148802/150‑1 (MQ=255)
                                                                                                                    gCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCg                              >  1:1399592/1‑150 (MQ=255)
                                                                                                                       ttGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCg                           >  2:6652463/1‑150 (MQ=255)
                                                                                                                         gcgcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTc                         <  1:2282238/150‑1 (MQ=255)
                                                                                                                            ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc                      >  2:9917539/1‑150 (MQ=255)
                                                                                                                             cAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAcc                     <  1:9156707/150‑1 (MQ=255)
                                                                                                                             cAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAcc                     >  2:10623842/1‑150 (MQ=255)
                                                                                                                                     gtcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGTTCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgc             >  1:8737927/4‑150 (MQ=255)
                                                                                                                                        gaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGc          >  1:5804333/1‑150 (MQ=255)
                                                                                                                                        gaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAGCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGc          >  2:8386055/1‑150 (MQ=255)
                                                                                                                                         aggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGcc         >  2:7788832/1‑150 (MQ=255)
                                                                                                                                         aggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGcc         >  1:10014859/1‑150 (MQ=255)
                                                                                                                                          ggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCa        >  2:5530667/1‑150 (MQ=255)
                                                                                                                                           gttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCAc       <  1:8150319/150‑1 (MQ=255)
                                                                                                                                           gttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCAc       <  2:8657797/150‑1 (MQ=255)
                                                                                                                                             tgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACtt     <  2:10576218/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051266‑3051559

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