Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,020,762 C→A 17.3% A279E (GCG→GAG)  helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,020,7620CA17.3% 731.7 / 107.2 243A279E (GCG→GAG) helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base C (88/113);  new base A (25/17);  total (113/130)
Fisher's exact test for biased strand distribution p-value = 8.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 264 total aligned reads displayed.

TGCAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGC  >  CP009273/1020613‑1020911
                                                                                                                                                     |                                                                                                                                                     
tGCAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggc                                                                                                                                                       >  1:4322672/1‑150 (MQ=255)
tGCAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggc                                                                                                                                                       >  1:269014/1‑150 (MQ=255)
 gCAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggcg                                                                                                                                                      >  2:3895989/1‑150 (MQ=255)
  cAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggagc                                                                                                                                                     >  1:3960660/1‑150 (MQ=255)
    ggCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAGGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggcgctg                                                                                                                                                   >  2:970685/1‑150 (MQ=255)
    ggCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATgcgctggcgctg                                                                                                                                                   <  1:2923166/150‑1 (MQ=255)
      cTGGTTGATGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCa                                                                                                                                                 <  1:1708232/150‑1 (MQ=255)
        ggTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCata                                                                                                                                               >  2:1612766/1‑150 (MQ=255)
        ggTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCata                                                                                                                                               >  1:3374352/1‑150 (MQ=255)
         gTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCatat                                                                                                                                              >  1:2379001/1‑150 (MQ=255)
                gcgcGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTcag                                                                                                                                       <  1:1838283/150‑1 (MQ=255)
                  gcgTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTcagca                                                                                                                                     <  2:4018185/150‑1 (MQ=255)
                    gtggtgAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAgg                                                                                                                                   <  2:1843638/150‑1 (MQ=255)
                       gtgAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGgca                                                                                                                                >  1:1038688/1‑150 (MQ=255)
                       gtgAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGgca                                                                                                                                <  1:3726754/150‑1 (MQ=255)
                        tgAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGgcag                                                                                                                               >  2:1711304/1‑150 (MQ=255)
                           aGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCaa                                                                                                                            >  2:3784168/1‑150 (MQ=255)
                               tCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGGTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCaaaaaa                                                                                                                        >  1:3428955/1‑150 (MQ=255)
                                 cccTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGt                                                                                                                      <  1:872735/150‑1 (MQ=255)
                                  ccTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGtt                                                                                                                     <  1:928766/150‑1 (MQ=255)
                                  ccTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGtt                                                                                                                     <  2:3013454/150‑1 (MQ=255)
                                   cTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTc                                                                                                                    <  1:3547543/150‑1 (MQ=255)
                                   cTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTc                                                                                                                    >  1:3532637/1‑150 (MQ=255)
                                     gAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCg                                                                                                                  <  1:4090856/150‑1 (MQ=255)
                                      aGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGa                                                                                                                 >  2:3211885/1‑150 (MQ=255)
                                         aaaTTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCTCGCAAGTTTCATGCGCTGGCGCTGCATATTTTTCAGCAGGGCAGCAAAAATTTTCCTATAg                                                                                                              <  2:562707/150‑1 (MQ=255)
                                         aaaTTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAg                                                                                                              >  1:2995143/1‑150 (MQ=255)
                                         aaaTTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAg                                                                                                              <  2:4310035/150‑1 (MQ=255)
                                          aaTTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGt                                                                                                             >  1:3992586/1‑150 (MQ=255)
                                           aTTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTc                                                                                                            <  1:1177666/150‑1 (MQ=255)
                                            ttttATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCa                                                                                                           >  1:4418939/1‑150 (MQ=255)
                                             tttATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAg                                                                                                          <  2:724199/150‑1 (MQ=255)
                                             tttATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAg                                                                                                          <  2:2335534/150‑1 (MQ=255)
                                               tATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCa                                                                                                        <  1:2014447/150‑1 (MQ=255)
                                                 ttGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCaaa                                                                                                      >  2:2094821/1‑150 (MQ=255)
                                                  tGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAAc                                                                                                     <  2:11740/150‑1 (MQ=255)
                                                  tGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAAc                                                                                                     <  1:4109442/150‑1 (MQ=255)
                                                     tGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTgg                                                                                                  <  1:250983/150‑1 (MQ=255)
                                                     tGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTgg                                                                                                  <  1:490596/150‑1 (MQ=255)
                                                      ggCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGa                                                                                                 <  1:4166507/150‑1 (MQ=255)
                                                      ggCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGa                                                                                                 <  2:3176985/150‑1 (MQ=255)
                                                         gTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGaaaa                                                                                              >  1:4237164/1‑150 (MQ=255)
                                                         gTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGaaaa                                                                                              >  2:2621032/1‑150 (MQ=255)
                                                           ttGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATg                                                                                            <  1:3764020/150‑1 (MQ=255)
                                                            tGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGa                                                                                           >  1:2633642/1‑150 (MQ=255)
                                                            tGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGa                                                                                           >  2:4489115/1‑150 (MQ=255)
                                                            tGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGa                                                                                           >  1:371156/1‑150 (MQ=255)
                                                                cGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATAcc                                                                                       <  2:56105/150‑1 (MQ=255)
                                                                cGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATAcc                                                                                       <  1:214905/150‑1 (MQ=255)
                                                                   aaaGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCt                                                                                    >  2:1094233/1‑150 (MQ=255)
                                                                   aaaGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCt                                                                                    >  1:2855031/1‑150 (MQ=255)
                                                                    aaGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTg                                                                                   >  1:2281114/1‑150 (MQ=255)
                                                                    aaGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTg                                                                                   >  1:2833346/1‑150 (MQ=255)
                                                                      gCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGcc                                                                                 <  1:1409825/150‑1 (MQ=255)
                                                                        cGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCg                                                                               <  2:3241330/150‑1 (MQ=255)
                                                                         gCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGt                                                                              >  1:1315172/1‑150 (MQ=255)
                                                                          cTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTc                                                                             >  1:2213466/1‑150 (MQ=255)
                                                                            gAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCAt                                                                           >  1:1235017/1‑150 (MQ=255)
                                                                             aaGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATg                                                                          >  2:3515060/1‑150 (MQ=255)
                                                                              agagATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGa                                                                         >  2:2538581/1‑150 (MQ=255)
                                                                               gagaTGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGaa                                                                        >  2:3751793/1‑150 (MQ=255)
                                                                                  aTGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAActc                                                                     >  2:1223567/1‑150 (MQ=255)
                                                                                  aTGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAActc                                                                     >  2:4157722/1‑150 (MQ=255)
                                                                                   tGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAActct                                                                    <  1:83030/150‑1 (MQ=255)
                                                                                    ggACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCACACTGGAAAATGATACCGCTGCCCGTCATGAACTCtt                                                                   >  2:3762648/1‑150 (MQ=255)
                                                                                     gACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCttt                                                                  <  1:4026281/150‑1 (MQ=255)
                                                                                        gAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTAtt                                                               <  1:678745/150‑1 (MQ=255)
                                                                                            ggATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTg                                                           >  2:1945939/1‑150 (MQ=255)
                                                                                             gATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGa                                                          >  1:3590687/1‑150 (MQ=255)
                                                                                              aTTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAg                                                         <  2:1409301/150‑1 (MQ=255)
                                                                                               ttCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGt                                                        <  1:3112157/150‑1 (MQ=255)
                                                                                               ttCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGt                                                        <  2:3336490/150‑1 (MQ=255)
                                                                                               ttCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGt                                                        <  2:2715962/150‑1 (MQ=255)
                                                                                                    gAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGgcgc                                                   <  1:3784259/150‑1 (MQ=255)
                                                                                                        ggCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCgcaagc                                               <  2:3629138/150‑1 (MQ=255)
                                                                                                         gCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCgcaagca                                              <  1:988272/150‑1 (MQ=255)
                                                                                                          cTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCgcaagcaa                                             <  1:247552/150‑1 (MQ=255)
                                                                                                          cTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCgcaagcaa                                             >  1:186579/1‑150 (MQ=255)
                                                                                                            aCATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATCGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATg                                           <  2:874932/150‑1 (MQ=255)
                                                                                                            aCATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATg                                           <  1:979069/150‑1 (MQ=255)
                                                                                                                aCCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGc                                       <  1:2461745/150‑1 (MQ=255)
                                                                                                                aCCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGc                                       <  1:195230/150‑1 (MQ=255)
                                                                                                                 ccGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCg                                      >  1:3700391/1‑150 (MQ=255)
                                                                                                                 ccGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCg                                      <  1:582823/150‑1 (MQ=255)
                                                                                                                  cGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGa                                     >  2:813503/1‑150 (MQ=255)
                                                                                                                       aCATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAaaaga                                <  2:434653/150‑1 (MQ=255)
                                                                                                                        cATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAaaagaa                               >  2:247251/1‑150 (MQ=255)
                                                                                                                        cATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAaaagaa                               >  1:1453869/1‑150 (MQ=255)
                                                                                                                           tACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAgc                            <  1:4379535/150‑1 (MQ=255)
                                                                                                                           tACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAgc                            >  2:1329528/1‑150 (MQ=255)
                                                                                                                            aCCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCACTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAgcg                           <  2:1768614/150‑1 (MQ=255)
                                                                                                                                cacgcacgTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAgg                       >  2:1451934/1‑150 (MQ=255)
                                                                                                                                     acgTTTCATGCGCTGGAGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGacg                  >  2:169088/1‑148 (MQ=255)
                                                                                                                                      cgTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAgg                 >  2:4025105/1‑150 (MQ=255)
                                                                                                                                      cgTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAgg                 >  2:1745875/1‑150 (MQ=255)
                                                                                                                                       gTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAggg                >  1:4104111/1‑150 (MQ=255)
                                                                                                                                               cgctggagctgCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTggcggc        >  1:3902132/1‑150 (MQ=255)
                                                                                                                                                  tggcgctgCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAAt     >  1:3255336/1‑150 (MQ=255)
                                                                                                                                                   ggcgctgCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATg    <  1:100492/150‑1 (MQ=255)
                                                                                                                                                     cgctgCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGc  <  2:2960656/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAGGCTGGTTGCTGGCGCGTGGTGAAGCGTCCCCTGAGCAAATTTTATTGCTGGCGTTTGGTCGCAAAGCCGCTGAAGAGATGGACGAGCGGATTCGCGAACGGCTACATACCGAAGACATTACCGCACGCACGTTTCATGCGCTGGCGCTGCATATTATTCAGCAGGGCAGCAAAAAAGTTCCGATAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGC  >  CP009273/1020613‑1020911

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: ‑