Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,324,066 C→A 39.8% E350* (GAA→TAA)  obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,324,0660CA39.8% 277.8 / 259.1 226E350* (GAA→TAA) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (81/55);  new base A (60/30);  total (141/85)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 227 total aligned reads displayed.

TGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAG  >  CP009273/3323917‑3324212
                                                                                                                                                     |                                                                                                                                                  
tGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTc                                                                                                                                                    <  1:2095413/150‑1 (MQ=255)
    ttGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACt                                                                                                                                                <  1:2996295/150‑1 (MQ=255)
    ttGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACt                                                                                                                                                <  2:3385808/150‑1 (MQ=255)
     tGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACtt                                                                                                                                               <  1:1000734/150‑1 (MQ=255)
      ggCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACttt                                                                                                                                              >  1:742579/1‑150 (MQ=255)
       gCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTtc                                                                                                                                             <  1:75870/150‑1 (MQ=255)
       gCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTtc                                                                                                                                             <  2:1063443/150‑1 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTtctc                                                                                                                                           >  1:288562/1‑150 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTtctc                                                                                                                                           <  2:2471021/150‑1 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTtctc                                                                                                                                           <  1:727182/150‑1 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTtctc                                                                                                                                           <  2:389056/150‑1 (MQ=255)
          tttATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTctct                                                                                                                                          >  1:1024361/1‑150 (MQ=255)
            tATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTgg                                                                                                                                        >  1:1316487/1‑150 (MQ=255)
             atcatcAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGc                                                                                                                                       >  1:1117250/1‑150 (MQ=255)
              tcatcaGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCt                                                                                                                                      >  1:1566901/1‑150 (MQ=255)
               catcaGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTg                                                                                                                                     >  1:2652568/1‑150 (MQ=255)
                 tcaGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGtt                                                                                                                                   <  2:1862848/150‑1 (MQ=255)
                   aGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTc                                                                                                                                 <  2:593679/150‑1 (MQ=255)
                     agATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGACCCAGTCGTCATCCCAGTCATCATCATCCTCTTCAGCAATCTCTTCAATCTGCTGGCGATGATAATCATACCACATGAATTAGACTTTCACTGGCTGTTTAGc                                                                                                                               >  2:244472/2‑150 (MQ=255)
                      gATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGct                                                                                                                              >  1:2754295/1‑150 (MQ=255)
                             gCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAg                                                                                                                       >  1:3538133/1‑150 (MQ=255)
                               ttGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGcc                                                                                                                     >  1:878333/1‑150 (MQ=255)
                                tGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCt                                                                                                                    <  1:1163948/150‑1 (MQ=255)
                                 gTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTg                                                                                                                   <  2:237454/150‑1 (MQ=255)
                                 gTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTg                                                                                                                   <  1:384835/150‑1 (MQ=255)
                                  tAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGc                                                                                                                  >  1:67337/1‑150 (MQ=255)
                                   aaaTGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCa                                                                                                                 >  2:1214823/1‑150 (MQ=255)
                                   aaaTGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCa                                                                                                                 <  1:1673840/150‑1 (MQ=255)
                                     aTGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACg                                                                                                               >  1:2522568/1‑150 (MQ=255)
                                     aTGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACg                                                                                                               >  1:2504626/1‑150 (MQ=255)
                                       gAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGAc                                                                                                             >  1:1580262/1‑150 (MQ=255)
                                        aaCTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGAcc                                                                                                            >  2:697457/1‑150 (MQ=255)
                                              aCGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGttt                                                                                                      >  1:4027573/1‑150 (MQ=255)
                                               cGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGtttt                                                                                                     <  2:298179/150‑1 (MQ=255)
                                                    tcttcGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATg                                                                                                >  1:3147796/1‑150 (MQ=255)
                                                        cgtcgtCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATaa                                                                                            <  2:327211/150‑1 (MQ=255)
                                                         gtcgtcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATaaa                                                                                           >  1:1681766/1‑150 (MQ=255)
                                                          tcgtcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAg                                                                                          <  2:2111499/150‑1 (MQ=255)
                                                            gtcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGt                                                                                        >  1:1517747/1‑150 (MQ=255)
                                                            gtcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGt                                                                                        <  2:1901277/150‑1 (MQ=255)
                                                            gtcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGt                                                                                        >  1:4051274/1‑150 (MQ=255)
                                                             tcTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGtc                                                                                       <  2:1396474/150‑1 (MQ=255)
                                                               ttCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGtcat                                                                                     <  1:1721270/150‑1 (MQ=255)
                                                                  gTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCACGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcac                                                                                  <  1:9479/150‑1 (MQ=255)
                                                                  gTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcac                                                                                  >  1:1916589/1‑150 (MQ=255)
                                                                  gTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcac                                                                                  >  1:2406441/1‑150 (MQ=255)
                                                                  gTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAAACTCTTCAAGCTGCTGGCGATGATAATAATCACACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATaac                                                                                  >  1:2309788/1‑150 (MQ=255)
                                                                   tCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcaca                                                                                 >  1:1120459/1‑150 (MQ=255)
                                                                     ccAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATc                                                                               <  2:3793613/150‑1 (MQ=255)
                                                                      cAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATcc                                                                              >  2:318147/1‑150 (MQ=255)
                                                                              aTCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCagag                                                                      >  1:2779696/1‑150 (MQ=255)
                                                                                   aGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCtt                                                                 <  2:1387196/150‑1 (MQ=255)
                                                                                    gTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCttt                                                                >  1:1148249/1‑150 (MQ=255)
                                                                                     tCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACCTGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTc                                                               >  2:2705189/1‑150 (MQ=255)
                                                                                     tCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTc                                                               >  2:1240584/1‑150 (MQ=255)
                                                                                      cTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCa                                                              >  1:774708/1‑150 (MQ=255)
                                                                                      cTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCa                                                              >  2:2012769/1‑150 (MQ=255)
                                                                                         catcatCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGc                                                           >  1:3776589/1‑150 (MQ=255)
                                                                                          atcatcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGcc                                                          >  2:3753222/1‑150 (MQ=255)
                                                                                          atcatcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGcc                                                          >  1:3558778/1‑150 (MQ=255)
                                                                                            catcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCa                                                        <  2:2381918/150‑1 (MQ=255)
                                                                                             atcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAg                                                       >  1:2354361/1‑150 (MQ=255)
                                                                                                   ttCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACt                                                 <  1:2862203/150‑1 (MQ=255)
                                                                                                    tCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTc                                                <  1:59219/150‑1 (MQ=255)
                                                                                                    tCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTc                                                <  1:98780/150‑1 (MQ=255)
                                                                                                     cAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcg                                               <  1:866422/150‑1 (MQ=255)
                                                                                                      aGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgc                                              >  1:3629874/1‑150 (MQ=255)
                                                                                                       gCAATCTCTTCACGCTGCTGGCGATGATAATCATCCCACATTAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgcc                                             <  1:912711/150‑1 (MQ=255)
                                                                                                             tctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCagag                                       <  2:3631976/150‑1 (MQ=255)
                                                                                                              ctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAgaga                                      >  2:2319955/1‑150 (MQ=255)
                                                                                                              ctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAgaga                                      >  1:2567860/1‑150 (MQ=255)
                                                                                                                  aaGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAg                                  >  1:2173274/1‑150 (MQ=255)
                                                                                                                  aaGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAg                                  >  1:2957369/1‑150 (MQ=255)
                                                                                                                  aaGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAg                                  >  2:2521969/1‑150 (MQ=255)
                                                                                                                   aGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGa                                 >  1:924156/1‑150 (MQ=255)
                                                                                                                      tgctgGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGataa                              >  1:1163473/1‑150 (MQ=255)
                                                                                                                            cGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTa                        <  1:4003489/150‑1 (MQ=255)
                                                                                                                            cGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTa                        <  2:182896/150‑1 (MQ=255)
                                                                                                                                gatAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCtt                    <  2:2156352/150‑1 (MQ=255)
                                                                                                                                 atAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTc                   <  2:4119610/150‑1 (MQ=255)
                                                                                                                                   aaTCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccc                 >  1:1036176/1‑150 (MQ=255)
                                                                                                                                   aaTCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATACCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccc                 >  1:2440956/1‑150 (MQ=255)
                                                                                                                                   aaTCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccc                 >  2:659048/1‑150 (MQ=255)
                                                                                                                                    atcatcCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccca                >  2:302488/1‑150 (MQ=255)
                                                                                                                                    atcatcCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccca                >  1:2242660/1‑150 (MQ=255)
                                                                                                                                     tcatcCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccag               >  1:4115637/1‑150 (MQ=255)
                                                                                                                                     tcatcCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccag               >  1:4115653/1‑150 (MQ=255)
                                                                                                                                       atcCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcc             >  2:417055/1‑150 (MQ=255)
                                                                                                                                          ccACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcccag          <  2:756215/150‑1 (MQ=255)
                                                                                                                                          ccACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcccag          >  2:609010/1‑150 (MQ=255)
                                                                                                                                            acaTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAgcg        >  1:423346/1‑150 (MQ=255)
                                                                                                                                            acaTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAgcg        >  2:1933205/1‑150 (MQ=255)
                                                                                                                                            acaTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAgcg        >  1:3765757/1‑150 (MQ=255)
                                                                                                                                             caTTAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAgcgc       >  1:228189/1‑150 (MQ=255)
                                                                                                                                              aTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGcc      >  1:2263365/1‑150 (MQ=255)
                                                                                                                                               tGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCt     <  2:2047213/150‑1 (MQ=255)
                                                                                                                                                 aaTTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCa   >  2:4053571/1‑150 (MQ=255)
                                                                                                                                                  aTTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAg  <  1:205247/150‑1 (MQ=255)
                                                                                                                                                  aTTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAg  >  2:3502492/1‑150 (MQ=255)
                                                                                                                                                  aTTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAg  >  1:3290669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAG  >  CP009273/3323917‑3324212

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: