Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,390,993 T→G 17.6% E566D (GAA→GAC)  bvgS ← Virulence sensor protein BvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,390,9930TG17.6% 155.6 / 2.7 110E566D (GAA→GAC) bvgSVirulence sensor protein BvgS
Reads supporting (aligned to +/- strand):  ref base T (52/34);  new base G (17/6);  total (70/40)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GGCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCA  >  ACB143/2390864‑2391125
                                                                                                                                 |                                                                                                                                     
ggCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCTGGCAACCCTCGCGGGCGCGCACATACATGGGGGGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGgcg                                                                                                                                 <  4:13283/136‑1 (MQ=255)
ggCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATTAATGCCAACTGGTCGTTCAGCGCGCGTTCGgcg                                                                                                                                 <  7:16209/136‑1 (MQ=255)
 gCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACTGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGAAGCGCGCGGTGGgccc                                                                                                                                >  5:117004/1‑134 (MQ=255)
    tGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAAGGCGCGGTCGGCgcgct                                                                                                                             >  6:364826/1‑136 (MQ=255)
       aCGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCgcgctggc                                                                                                                          >  2:101033/1‑136 (MQ=255)
             tCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGcgctggcgctgg                                                                                                                     <  2:204296/135‑1 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGcgctggcgcggggg                                                                                                                   >  2:15738/1‑134 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCgctggcgctggc                                                                                                                    >  8:187726/1‑135 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGctggcgctggcg                                                                                                                   >  6:373865/1‑136 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGctggcgcgggcg                                                                                                                   >  1:45826/1‑136 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCCCGTTCGGCGCGctggcgctggcg                                                                                                                   >  8:193594/1‑136 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACAGGTCGTTAAGCCCGCGTTCGGCGCgcggtcgcgggcc                                                                                                                   >  3:129870/1‑135 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGAAGGGCGCGAATGAATGCCAACTGGTCGTTCATCGCGCGTTCGGCGCGctggcgctggcg                                                                                                                   >  6:165002/1‑136 (MQ=255)
              cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGGGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGcgctggcgctgg                                                                                                                     <  1:24405/134‑1 (MQ=255)
                  gTAGCTGTAGTTGCAACTTTTCCGGCAACCCTCGCGGTCGCGCACATACAGGGGGGGGGGCGTACCGTTGAGCAGGGCGCGCAGGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCGGGCGCTGGCGAATc                                                                                                               <  4:244604/136‑1 (MQ=255)
                        gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGCGCCCGTGGGGCCCGCGGGCGCGGGGGGATC‑Tggcgg                                                                                                        >  1:391975/1‑132 (MQ=255)
                        gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGCCCGCGGTAGGGGCGCGGGCGGGGGCGAATA‑Tgccgg                                                                                                        >  3:358541/1‑136 (MQ=255)
                        gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTACGCGCGCGGTCGGCGCGCTGGGGCTGGCGAAAA‑Tgcagg                                                                                                        >  7:287594/1‑136 (MQ=255)
                           gTTGCAAGTTTTCAGGCAACCCTCGCGGTGGCGCACATTCTTGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCTGCGCGCGTTGGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcg                                                                                                      <  2:1633/135‑1 (MQ=255)
                               cacattttCAGGACACCCTCCAGGTCGCGAAAATACCTGGGGTGGGTGTTAAGGTTCAGCAGGGCGCTCATTAATCCCAACTGGTCGTTCCGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggagcaggc                                                                                                 <  8:335992/132‑1 (MQ=255)
                               cAACTTTTCAGGCACCCCTCGCGGTCCCGCACATACATGGGGTGCGGCGGACCGTTGAGCAGGCCGCGCATGAATCCAAACTGGTCGTTCAGCGCGCGGTCGGCGCGGTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  4:157488/136‑1 (MQ=255)
                               cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  7:284107/136‑1 (MQ=255)
                               cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  5:110909/136‑1 (MQ=255)
                               cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGGGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  3:357317/136‑1 (MQ=255)
                                aaCTTTTCAGGAAACCCTCGGGGTCGCCCAAATACTGGGGGGGGGGCGTACGGTTGAGCTGGGCGGGCTTGAATCCCAAGTGGGGGTTCTGGGCGGGTTCGGCGCGCGGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  4:5148/135‑1 (MQ=255)
                                     ttCAGGCAACCCTCGCGGTCGCGCACATACATCGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCa                                                                                            <  5:721/135‑1 (MQ=255)
                                           cAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCaga                                                                                     <  6:391030/136‑1 (MQ=255)
                                                   gcgGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGGCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGcc                                                                             <  1:17507/136‑1 (MQ=255)
                                                   gcgGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCAGGAAGGCCAACTGGTGGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGcc                                                                             <  4:128223/136‑1 (MQ=255)
                                                             acatacttGGGGGGGGGCGTACCGTTGCGCCGGGCGCGCATGAATGCCAACTGGTGGTTCTGCGGGCGTTGGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACaaaa                                                                   <  3:51346/136‑1 (MQ=255)
                                                                  ctggggggggggCGTACGTTTGCGCTGGGCGCGCATTAATCCCACGTGGTCGTTCCGCGGGGGGGCGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGAAACAAACGCACg                                                              <  6:358418/133‑1 (MQ=255)
                                                                        gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCCCGCGATAGGCGCGCCGGCGCTGGCGAATC‑TGCACGCGCACGCGGGCCTTCCAGAGCCGGGCCAGCAACAAAAGCACGGCCGTg                                                        >  6:83088/1‑136 (MQ=255)
                                                                        gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCCGCGCGCTGGCGCTGGCGAATA‑TGAAGGCGCAGGCGGGCATTCCAGAGCAGGGCCCGCCACAAAAGCACGGCCCTg                                                        >  1:288252/1‑136 (MQ=255)
                                                                        gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCCGTCGGCGCGCTGGCGCGGGCGAATA‑TGCAGGCGCAGGCGGGCATTCCAGAGCACGGCCCCCAAAAAAAGCAccgccgag                                                        >  5:66057/1‑134 (MQ=255)
                                                                        gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCAAACTGGTCGTTAAGCGCGCGTTCGGCGCGCTTGCCCAGCCGAAAA‑TGCAGGACCACGACGGCATTACAGAGAAGTGCCACCAAAAAAAACAAGGACgat                                                        >  6:362033/1‑134 (MQ=255)
                                                                          ggggCGTACCGTTGAGCAGGGCGCGCATGAATGCAAACTGGTCGTTCAGCCCGCGTTCGGCGCGCTGGCGATGGCGAAAC‑TGAAGGCGCCGGCGGGCATTACAGAGCAGGGCAAACAAAAAAATCCCGGACGTgga                                                      >  2:121507/1‑135 (MQ=255)
                                                                           ggccGTCCCGTGGCGCAGGCCGCGCATCAAGCCAAACGGGGCGTGCCCCGCGGGGGCGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGACCAGGGCCAGCAACAAAAGCACGGCCGTGGcc                                                     <  1:322177/133‑1 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCCTTACAGAGCAGGGGCAGCCACAAAAGCCCGGCCGTGGccc                                                    >  1:52889/1‑135 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCCGCAACAAAAGCACGGCCGGGGccc                                                    >  6:347976/1‑135 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGccc                                                    >  3:110790/1‑135 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCa                                                    >  8:29284/1‑136 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTACAGAGCAGGGCCAGCCACAAAAGCACCGCCGTGGccc                                                    >  3:151513/1‑135 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTAGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCCGCAACAAAAGCACGGCCGTGGccc                                                    >  5:167214/1‑135 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCCGGCGGGCATTGCAGAGCAGGGCCAGCAACAAAAGCCCGGCCGTGGCCa                                                    >  6:346315/1‑136 (MQ=255)
                                                                            ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATA‑TGCAGGCGCAGGGGGGCCGTACAGAGCCGGGCCAGCAACAAAAGCACGGCCGTGGccc                                                    >  6:5549/1‑135 (MQ=255)
                                                                              cGTACCGTTGCGCAGGGCGCGCATGAATGCCAACTGGTCGTGCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGc                                                  <  7:89711/136‑1 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGAGCGCGTTCGGCGCGCGGGCGCTGGGGATTA‑TGCAGGCGCATGGGGGCAGTACAAACCATGGGCAGCCAGAAAAGAACGGCCGTTGCCAGc                                                  >  4:27018/1‑136 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGCATC‑TGCAGGCGCAGGCGGGCCTTCCAGAGCAGGGCCAGCAACAAAAGCACCGCCGTGGCCCGc                                                  >  3:128162/1‑136 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGCGCACGGCCCGCAACAAAATCACCGCCGTGGCCacc                                                  >  3:329213/1‑134 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCAGGCGCTGGCGAATC‑TTCAGGCGCAGGCGGGCATTCCACAGCAGGGCCAGAAACAAAAGCACCGCCGGGGcccag                                                  >  3:239894/1‑133 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCCCGCGGTCGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCCGGCCGGCATTCCAGAGCCGGGCCAGCAACAAAAGCCCCGCCCTGGcccg                                                   >  2:55159/1‑133 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCACCGCGCGTTCGGCGCGCTGGCGCTGGCGAATA‑AGCACGCGCAGGCGGGCATTACACAGCAGGGCCAGCAACAAAAGCACCGCCGTGGCCacc                                                  >  1:251032/1‑134 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCACCGCGCGTTCGGCGCGCTGGCGCGGGCGAATC‑TGCAGGCGCCGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCCCGGCCGTGGCCAGc                                                  >  3:125513/1‑136 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGGAGCGCGCGTTGGGCGCGCGGGCGCTGGGGGCTC‑TGCCGGCGCCGGCGGGGCGTGCACAGCAGGGGCCTCAACAAAACCCAGGCCGTgggccgc                                                  >  6:44678/1‑131 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGGAGAGCCCGTTGGGCGCGCTGGGGCGGGGGAATCTTGCA‑GCGCGCGCGGGGCGTCAACCGCACGGCCAGCCACAAAACCACGGCCGTGGcgcgc                                                  >  4:204583/1‑132 (MQ=255)
                                                                              cGTACCGTTGAGCAGGGCGCGCATGAATGACAACTGGTCGTTGAGCGCCCGGGCGGCCCGCGGGCGGGCGGCCCGC‑TGCGCGCGCAGGCGGGCTTGCCACAGCAGGGCCCGCCACAAAAGCACGGCCCTGGccccc                                                  >  4:63591/1‑133 (MQ=255)
                                                                               gTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTAGGCGCGCGGGCGCTGGCGAATC‑AGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCACCAACAAAAGCACGGCCGTGGCCCGca                                                 >  1:245367/1‑136 (MQ=255)
                                                                                tACCGTTGCGCAGGGCGCGCATGAATGCCAACTGGTCGTTCCGCGCGCGTTGGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGcac                                                <  7:70455/136‑1 (MQ=255)
                                                                                   ggtggAGCAGGGCGCGCATGAATGCCAACTGGTCGTGCAGCGCGCGTGCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGcaccac                                             <  5:212802/132‑1 (MQ=255)
                                                                                    gTTGCGCAGGGCGCGCTTTAATCCAAACTGGGGGTTCAGCGCGCGTGCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAACGCACGGCCGTGGCCAGCaccacc                                            <  5:315832/136‑1 (MQ=255)
                                                                                    gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGCGCAGGGCCAGCAACAAAAGCCCCGCCGTGGCCAGcacccc                                             >  2:44531/1‑133 (MQ=255)
                                                                                    gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACaaaa                                                                   <  5:258657/113‑1 (MQ=255)
                                                                                    gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATG‑TGCAGGCGCAGGCGGGGATTCCAGAGCACGGCCCGCAACAAAAGCCCGGCCGTGGCCCGCaccccc                                            >  3:31043/1‑136 (MQ=255)
                                                                                    gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACaaaa                                                                   >  6:258657/1‑113 (MQ=255)
                                                                                        aGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTGGGCGCCCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCCCCGCCGTGGCCAGCACCACCTg                                          >  4:385957/1‑134 (MQ=255)
                                                                                        aGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGGTAGGCGCGCGGGCGCGGGCGGAAA‑TGCAGGCGCAGGCGGGCCTTCCAGACCAGGGCCAGCAACAAAAGCACGGCCGTGGCCCCCCCCACCAGCa                                        >  6:147926/1‑136 (MQ=255)
                                                                                           aGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGCATC‑TGCCGGCGCCGGCGGGGATTGCAGACGAGGGCCCGCAACAAAAGCACCGCCCTGGCCAGCCCCACCAGCAggg                                     >  4:53904/1‑136 (MQ=255)
                                                                                           aGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAggg                                     >  7:289625/1‑136 (MQ=255)
                                                                                           aGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCCCTGGCGCTGGCGAATC‑TGCAGGCCCCGGCGGGCATTCCAGAGCACGGCCCGCAACAAAAGCACCGCCGTGGCCAGCACACCCAGCAggg                                     >  2:15338/1‑136 (MQ=255)
                                                                                           aGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTACAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAggg                                     >  2:32956/1‑136 (MQ=255)
                                                                                             ggCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCCCCACCTGCAGGGc                                    >  7:89778/1‑135 (MQ=255)
                                                                                             ggCGCGCATGAATGCCAACTGGTCGTGCAGCGCGCGTGCGGGGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGcc                                   <  6:73627/136‑1 (MQ=255)
                                                                                             ggCGCGCATGAATGCCAACTGGGCGTTCAGCGCGCGGTCGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAACGCACGGCCGTGGCCAGCACCACCTGCAGGGcc                                   <  1:317573/136‑1 (MQ=255)
                                                                                             ggCGCGCATGAATGCCAACGGGGCGTCCAGCGCGCGTGCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGcc                                   <  3:56321/136‑1 (MQ=255)
                                                                                             ggCGCGCATGAAGCCCACCGGGCCGTCCCGCGCGCGGCGGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGcc                                   <  3:242831/136‑1 (MQ=255)
                                                                                             ggCGCGCATCACGCCCACGTGGTGGGTCCGCGCGCGGGCGCGGCGCGGGCGCTGGCGAAGC‑TGCAGGCGCAGGCGGGCATTCCTGACCCGGGCCAGCAACAAACGCACGGCCGTGGCCAGCACCACCTGCAGGGcc                                   <  3:98852/136‑1 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCCGGCGCCGGCGGGCATTCCACAGCAGGGCCAGCAACAAAACCACGGCCGTGGCCCGCCCCCCCGGCCGGGCCacg                                >  1:109759/1‑134 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGCGCAGGGCCAGCAACAAAAGCACCGCCGTGGCCCGCACCACCTGCAGGGCCCgg                                >  8:307900/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCGGGCGCATC‑TGCAGGCGCCGGCGGGCATTCCAGAGCAGGGCCCGCAACAAAAGCACGGCCCGGGCCCGCCCCCCCAGCCGGGCCCgg                                >  2:296660/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCGGGCGAATC‑TGCAGGCGCAGGCGGGCATTACAGAGCAGGGCCAGCAACAAAAGCCCGGCCGTGGCCAGCCCCCCCAGGCGGGCCAgg                                >  5:97891/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCCGGGCCAGCAACAAAAGCACGGCCGTGGCCCGCCCCCCCTGCAGGGCCAgg                                >  2:217973/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGGTGGGCGCGCTGGCGCGGGCGAATA‑TGCAGGCGCCGGCGGGCAGTCCAGAGCCCGGCCAGCAACAAAAGCACCGCCCTGGCCCGCCCCCCCCGcccggcccgg                                >  4:17189/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCGGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCCGCCCCACCAGCCGGGCCCgg                                >  6:62135/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTCAACCCGCGGTCGGCGCGCAGGCGCGGGCGAATC‑TTCACGCGCAGGCGGGGATTCCAGAGCAGGGCCAGCAACAAAAGCACGCCCGTGGCCACCACCACCAGCCGGGCCCgg                                >  4:284327/1‑136 (MQ=255)
                                                                                                gcgcATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGCAGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGGATTCCAGCGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCCCCCCCGGCAGGGCCCgg                                >  2:147380/1‑136 (MQ=255)
                                                                                                  gcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTACAGAGCAGGGCCCGCAACAAAAGCACGGCCGTGGCCAGCCCCACCTGCAGGGCCAggcg                              >  5:326236/1‑136 (MQ=255)
                                                                                                  gcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCGGGCGAATC‑TGCAGGCGCAGGCGGGCAGTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGGGGCCCGCCCCACCAGCAGGGCCAggcg                              >  4:66169/1‑136 (MQ=255)
                                                                                                    aTGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAggcggc                            <  7:35684/136‑1 (MQ=255)
                                                                                                    aTGAATCCCAACTGTGCGTCCTGCGCGGGTGCGGGGGGCGGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATACCAGAGCAGGGCCAGCAACAAACGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAggcggc                            <  6:117004/136‑1 (MQ=255)
                                                                                                      gAATGCCACCTGGTCGTTCAGCGCGCGTGCGGCGCGCTGGCGCGGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGAAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCt                           <  4:366030/135‑1 (MQ=255)
                                                                                                       aaTGCCAACGGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGt                         <  6:215274/136‑1 (MQ=255)
                                                                                                           ccAACTGGGCGTTCAGCGCGGGTGCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGAAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGc                      <  2:356571/135‑1 (MQ=255)
                                                                                                               cTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 >  6:354094/1‑56 (MQ=255)
                                                                                                               cTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑Tgcaggcgcaggc                                                                                                 <  5:354094/56‑1 (MQ=255)
                                                                                                                 ggTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTg                 >  7:88855/1‑134 (MQ=255)
                                                                                                                    cGTTAAGCGCGCGTTAGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCCGGGCCCGCAACAAAAGCACCGCCCTGGCCAGCACCACCAGCATGGCCAGGCGGGTGTAGCCGTTGGAGtt            >  4:357365/1‑136 (MQ=255)
                                                                                                                             cgggTTGGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCTGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTc   <  3:93497/133‑1 (MQ=255)
                                                                                                                             cgcgTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTc   <  1:279304/136‑1 (MQ=255)
                                                                                                                              gcgTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCa  <  5:347976/136‑1 (MQ=255)
                                                                                                                              gcgTTCGGCGCGCTGGCGCAGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAAAAAAAGCACGGCCGTGGCCAGCCCCACCTGCAGGGCCAGGCGGCGGTAGCCTTTGGAGTTGTCAACcaca  >  6:344955/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                                                     
GGCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATC‑TGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCA  >  ACB143/2390864‑2391125

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