Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,055,246 T→A 8.4% L241Q (CTG→CAG)  yihM → putative sugar phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,055,2460TA 8.4% 750.2 / 45.6 234L241Q (CTG→CAG) yihMputative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (104/110);  new base A (12/8);  total (116/118)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAC  >  CP009273/4055097‑4055395
                                                                                                                                                     |                                                                                                                                                     
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       <  1:2781734/150‑1 (MQ=255)
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       >  2:272358/1‑150 (MQ=255)
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       >  2:2620140/1‑150 (MQ=255)
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       >  1:3187542/1‑150 (MQ=255)
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       >  1:1354908/1‑150 (MQ=255)
cGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCt                                                                                                                                                       >  2:1591390/1‑150 (MQ=255)
 gCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTg                                                                                                                                                      >  2:1220672/1‑150 (MQ=255)
    ttAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGt                                                                                                                                                   <  2:359903/150‑1 (MQ=255)
     tAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGtg                                                                                                                                                  <  1:2546970/150‑1 (MQ=255)
        aaaCGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGtgatg                                                                                                                                               <  1:1195999/150‑1 (MQ=255)
        aaaCGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGtgatg                                                                                                                                               >  2:4065059/1‑150 (MQ=255)
         aaCGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGtgatga                                                                                                                                              >  1:1656111/1‑150 (MQ=255)
            gATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGtgatgatga                                                                                                                                           >  1:1473551/1‑150 (MQ=255)
                       atatTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGa                                                                                                                                <  1:1322985/150‑1 (MQ=255)
                         atTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACg                                                                                                                              >  2:3555550/1‑150 (MQ=255)
                           tACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGc                                                                                                                            <  1:147169/150‑1 (MQ=255)
                           tACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGc                                                                                                                            >  1:813742/1‑150 (MQ=255)
                           tACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACGGc                                                                                                                            >  1:2488303/1‑150 (MQ=255)
                              ccAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCata                                                                                                                         >  2:1434614/1‑150 (MQ=255)
                              ccAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACGGCata                                                                                                                         <  2:809987/150‑1 (MQ=255)
                                aGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATg                                                                                                                       >  1:2665576/1‑150 (MQ=255)
                                 ggTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATgg                                                                                                                      <  2:2831555/150‑1 (MQ=255)
                                  gTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGc                                                                                                                     >  2:2598750/1‑150 (MQ=255)
                                  gTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGc                                                                                                                     >  2:515163/1‑150 (MQ=255)
                                     cATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTg                                                                                                                  >  2:2408860/1‑150 (MQ=255)
                                        atCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGgaa                                                                                                               >  2:2659429/1‑150 (MQ=255)
                                           aaaGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGgaagaa                                                                                                            >  1:66175/1‑150 (MQ=255)
                                            aaGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGaagaag                                                                                                           >  1:1025769/1‑150 (MQ=255)
                                               aTGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAgg                                                                                                        <  1:1665588/150‑1 (MQ=255)
                                                 gCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGtt                                                                                                      >  1:2367856/1‑150 (MQ=255)
                                                 gCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGtt                                                                                                      <  1:2865826/150‑1 (MQ=255)
                                                 gCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGtt                                                                                                      <  1:2865971/150‑1 (MQ=255)
                                                  ccTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTg                                                                                                     >  1:2046749/1‑150 (MQ=255)
                                                   cTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGa                                                                                                    <  1:1194120/150‑1 (MQ=255)
                                                   cTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGa                                                                                                    >  2:164739/1‑150 (MQ=255)
                                                     tGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGAtt                                                                                                  <  2:742412/150‑1 (MQ=255)
                                                     tGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGAtt                                                                                                  >  2:424113/1‑150 (MQ=255)
                                                      gATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTa                                                                                                 >  1:1075878/1‑150 (MQ=255)
                                                       aTCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTAc                                                                                                <  1:1885282/150‑1 (MQ=255)
                                                         cGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTa                                                                                              >  1:3674252/1‑150 (MQ=255)
                                                                    cGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGt                                                                                   >  1:431294/1‑150 (MQ=255)
                                                                    cGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGt                                                                                   <  1:1079636/150‑1 (MQ=255)
                                                                      ggTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTc                                                                                 >  1:180535/1‑150 (MQ=255)
                                                                       gTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTcc                                                                                <  2:2052870/150‑1 (MQ=255)
                                                                         ggCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCGGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGc                                                                              >  2:3399251/1‑150 (MQ=255)
                                                                          gCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCt                                                                             <  2:1375621/150‑1 (MQ=255)
                                                                            cTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCttt                                                                           <  2:4034726/150‑1 (MQ=255)
                                                                            cTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCttt                                                                           >  2:3088886/1‑150 (MQ=255)
                                                                              gggTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGa                                                                         <  2:1976040/150‑1 (MQ=255)
                                                                                gTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAg                                                                       <  2:980743/150‑1 (MQ=255)
                                                                                 tCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGa                                                                      <  1:2680327/150‑1 (MQ=255)
                                                                                    tAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGa                                                                   >  1:764033/1‑150 (MQ=255)
                                                                                     aaaGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGaa                                                                  <  1:2284723/150‑1 (MQ=255)
                                                                                     aaaGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGaa                                                                  <  1:1889721/150‑1 (MQ=255)
                                                                                         ccTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAgatg                                                              >  1:870099/1‑150 (MQ=255)
                                                                                         ccTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAgatg                                                              >  1:2708705/1‑150 (MQ=255)
                                                                                           tGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAgatgat                                                            <  2:3526950/150‑1 (MQ=255)
                                                                                           tGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAgatgat                                                            <  1:314521/150‑1 (MQ=255)
                                                                                           tGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAgatgat                                                            >  1:3645344/1‑150 (MQ=255)
                                                                                                 tcaggtcaggGTGATATGGCGTTCAAAGCGTTATTAAAGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGGCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCg                                                      <  2:2194262/150‑1 (MQ=255)
                                                                                                  taggtcaggGTGATATGACGTTCAAAGAGTTATTAACGCACCTTATCGGCCAGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCTTTCCGCTTTGAATACGATGATGATATTcctt                                                     <  2:1235301/149‑3 (MQ=255)
                                                                                                    ggtcaggGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTgg                                                   >  1:2231027/1‑150 (MQ=255)
                                                                                                    ggtcaggGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTgg                                                   <  2:3853239/150‑1 (MQ=255)
                                                                                                     gtcaggGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGa                                                  >  2:2362781/1‑150 (MQ=255)
                                                                                                      tcaggGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGAt                                                 <  2:3289384/150‑1 (MQ=255)
                                                                                                          ggTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCTCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCt                                             <  2:494033/150‑1 (MQ=255)
                                                                                                          ggTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCt                                             <  1:3934729/150‑1 (MQ=255)
                                                                                                          ggTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCt                                             <  1:2946147/150‑1 (MQ=255)
                                                                                                           gTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCtt                                            >  1:3462656/1‑150 (MQ=255)
                                                                                                           gTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCAGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCtt                                            <  1:2857969/150‑1 (MQ=255)
                                                                                                                 tGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc                                      <  2:4107309/150‑1 (MQ=255)
                                                                                                                    cGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGa                                   <  2:902234/150‑1 (MQ=255)
                                                                                                                     gTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGAt                                  >  2:2864314/1‑150 (MQ=255)
                                                                                                                      ttCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATg                                 >  1:408210/1‑150 (MQ=255)
                                                                                                                        cAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGGTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAg                               >  1:1799807/1‑150 (MQ=255)
                                                                                                                         aaaGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGt                              <  2:1993055/150‑1 (MQ=255)
                                                                                                                           aGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa                            >  2:3946269/1‑150 (MQ=255)
                                                                                                                           aGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa                            >  1:3021983/1‑150 (MQ=255)
                                                                                                                           aGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa                            >  1:3236719/1‑150 (MQ=255)
                                                                                                                             cGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGaaa                          <  2:1155643/150‑1 (MQ=255)
                                                                                                                                tattaACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAacac                       >  1:3723632/1‑150 (MQ=255)
                                                                                                                                   taACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCAt                    >  2:3679449/1‑150 (MQ=255)
                                                                                                                                    aaCGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCAtt                   <  1:91521/150‑1 (MQ=255)
                                                                                                                                      cGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGc                 <  2:3321758/150‑1 (MQ=255)
                                                                                                                                      cGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGc                 >  2:1785094/1‑150 (MQ=255)
                                                                                                                                      cGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGc                 >  2:1812644/1‑150 (MQ=255)
                                                                                                                                        cACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCa               >  2:3201482/1‑150 (MQ=255)
                                                                                                                                         aCCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAg              >  1:2958843/1‑150 (MQ=255)
                                                                                                                                           cTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaa            >  1:893051/1‑150 (MQ=255)
                                                                                                                                              aTCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAAt         >  2:3864973/1‑150 (MQ=255)
                                                                                                                                               tCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATc        <  2:607374/150‑1 (MQ=255)
                                                                                                                                                  gCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCAtt     >  1:2263688/1‑150 (MQ=255)
                                                                                                                                                  gCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCAtt     <  1:591673/150‑1 (MQ=255)
                                                                                                                                                   ccTGGGTGATGATTAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTAGATAACTATGCGCCGGCGTACCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATAGCCAGATGAGTGAAACACCATTTCCAGAAAATCAtat    >  2:69049/1‑148 (MQ=255)
                                                                                                                                                   ccTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCAttt    <  1:1738232/150‑1 (MQ=255)
                                                                                                                                                    cTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa   >  1:1433054/1‑150 (MQ=255)
                                                                                                                                                    cTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa   >  1:1433025/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc  >  2:1817164/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc  >  2:3454908/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc  >  2:1030016/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTTTAAAAACGATGGCACCGCATATTACCCAGGTCCATATCAAAGATGCCTTGATCGTTAAAGAACCGGGTGGCCTGGGTCATAAAGCCTGTATTTCAGGTCAGGGTGATATGCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAC  >  CP009273/4055097‑4055395

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: