Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,324,066 C→A 45.7% 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,0660CA45.7% 216.4 / inf 264E350* (GAA→TAA) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (67/76);  new base A (68/52);  total (135/129)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

TGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAGCGA  >  CP009273/3323917‑3324215
                                                                                                                                                     |                                                                                                                                                     
tGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTc                                                                                                                                                       <  2:1737248/150‑1 (MQ=255)
tGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTc                                                                                                                                                       >  1:3811774/1‑150 (MQ=255)
   tttGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGAc                                                                                                                                                    <  1:217922/150‑1 (MQ=255)
   tttGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGAc                                                                                                                                                    <  2:815687/150‑1 (MQ=255)
     tGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACtt                                                                                                                                                  <  2:3424382/150‑1 (MQ=255)
     tGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACtt                                                                                                                                                  >  1:805316/1‑150 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTtctc                                                                                                                                              >  2:1901541/1‑150 (MQ=255)
         gTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTtctc                                                                                                                                              >  1:4537143/1‑150 (MQ=255)
            tATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTgg                                                                                                                                           <  2:2475702/150‑1 (MQ=255)
            tATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTgg                                                                                                                                           >  1:2192941/1‑150 (MQ=255)
             atcatcAGTGATTAACGCTTGTAAATGAACTCACCGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGc                                                                                                                                          <  1:4254076/150‑1 (MQ=255)
             atcatcAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGc                                                                                                                                          <  1:530078/150‑1 (MQ=255)
                 tcaGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGtt                                                                                                                                      >  1:1284532/1‑150 (MQ=255)
                  caGTGATTAACTCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGttt                                                                                                                                     <  2:2397265/150‑1 (MQ=255)
                  caGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGttt                                                                                                                                     <  2:2364075/150‑1 (MQ=255)
                   aGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTc                                                                                                                                    <  1:3322054/150‑1 (MQ=255)
                    gTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCg                                                                                                                                   <  2:4373722/150‑1 (MQ=255)
                    gTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCg                                                                                                                                   <  1:1579419/150‑1 (MQ=255)
                       aTTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGctt                                                                                                                                <  1:4058316/150‑1 (MQ=255)
                       aTTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGctt                                                                                                                                >  1:211050/1‑150 (MQ=255)
                          aaCGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGcttctt                                                                                                                             <  1:3555785/150‑1 (MQ=255)
                          aaCGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGcttctt                                                                                                                             <  1:1608764/150‑1 (MQ=255)
                             gCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAg                                                                                                                          >  1:2235142/1‑150 (MQ=255)
                             gCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAg                                                                                                                          <  1:1305313/150‑1 (MQ=255)
                              cTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGc                                                                                                                         >  1:1076809/1‑150 (MQ=255)
                               ttGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACATTCTCTGGCTGTTTCGCTTCTTCAGcc                                                                                                                        >  2:4087791/1‑150 (MQ=255)
                               ttGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGcc                                                                                                                        <  2:2151074/150‑1 (MQ=255)
                                tGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCt                                                                                                                       <  1:984860/150‑1 (MQ=255)
                                   aaaTGAACTCAACGCCTTCTTCGTCGTCTACGTCCCAGTCGTCATCCCAGTCTTCATCCGCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCa                                                                                                                    <  1:767941/150‑1 (MQ=255)
                                      tGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGa                                                                                                                 <  2:3686727/150‑1 (MQ=255)
                                       gAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGAc                                                                                                                >  1:1874586/1‑150 (MQ=255)
                                          cTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCgg                                                                                                             <  2:3010688/150‑1 (MQ=255)
                                           tCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCggg                                                                                                            >  1:2120467/1‑150 (MQ=255)
                                            cAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGt                                                                                                           >  2:3308685/1‑150 (MQ=255)
                                                gCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTc                                                                                                       >  1:3323668/1‑150 (MQ=255)
                                                   ttcttcGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAAt                                                                                                    <  2:560638/150‑1 (MQ=255)
                                                    tcttcGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATg                                                                                                   >  1:4446103/1‑150 (MQ=255)
                                                      ttcGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGAt                                                                                                 >  1:483609/1‑150 (MQ=255)
                                                       tcgtcgTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCGGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATa                                                                                                <  2:1752011/150‑1 (MQ=255)
                                                                 cGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGtcatca                                                                                      <  1:4325224/150‑1 (MQ=255)
                                                                 cGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGtcatca                                                                                      >  2:2093231/1‑150 (MQ=255)
                                                                  gTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcac                                                                                     <  1:3889809/150‑1 (MQ=255)
                                                                   tCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATcaca                                                                                    <  1:1266790/150‑1 (MQ=255)
                                                                    cccAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCACTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACAt                                                                                   <  1:660466/150‑1 (MQ=255)
                                                                      cAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATcc                                                                                 <  2:3218229/150‑1 (MQ=255)
                                                                           gtcATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCcagca                                                                            <  1:2870825/150‑1 (MQ=255)
                                                                            tcATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCcagcag                                                                           <  1:3282703/150‑1 (MQ=255)
                                                                                    gTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACTTTAATTAGACTTTCTCTGGCTTTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCttt                                                                   <  2:3260440/150‑1 (MQ=255)
                                                                                    gTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCttt                                                                   <  2:2590453/150‑1 (MQ=255)
                                                                                     tCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTc                                                                  >  1:991400/1‑150 (MQ=255)
                                                                                         catcatCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGc                                                              >  1:2423438/1‑150 (MQ=255)
                                                                                          atcatcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGcc                                                             >  2:3140482/1‑150 (MQ=255)
                                                                                           tcatcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGccc                                                            <  2:4251771/150‑1 (MQ=255)
                                                                                            catcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCa                                                           <  2:2875725/150‑1 (MQ=255)
                                                                                             atcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATCATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGTTAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAg                                                          >  2:4541563/1‑150 (MQ=255)
                                                                                             atcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAg                                                          >  2:4214946/1‑150 (MQ=255)
                                                                                             atcCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAg                                                          >  1:4121515/1‑150 (MQ=255)
                                                                                                  ctTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCAc                                                     >  1:2028757/1‑150 (MQ=255)
                                                                                                   ttCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACt                                                    >  2:1971428/1‑150 (MQ=255)
                                                                                                    tCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTc                                                   <  2:2413974/150‑1 (MQ=255)
                                                                                                      aGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgc                                                 >  1:2835427/1‑150 (MQ=255)
                                                                                                      aGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGACATCACATCCCAGCAGAGATCTTTCACTCCCAGTCCACTcgc                                                 >  1:2835451/1‑150 (MQ=255)
                                                                                                       gaaaTCTCTTCATGCTTCTGGCGATGATAATCATCCCACATGAATTTTACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACAGGGTTTTCAATGATATAGGTCATCACATCCCAGCATAGATCTTTCACGCCCAGTCCACTcgcc                                                <  2:3958072/148‑1 (MQ=255)
                                                                                                       gCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgcc                                                <  1:3087292/150‑1 (MQ=255)
                                                                                                       gCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgcc                                                <  2:1634004/150‑1 (MQ=255)
                                                                                                       gCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAATCCACTcgcc                                                <  2:4330862/150‑1 (MQ=255)
                                                                                                        cAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTcgccg                                               >  2:1311001/1‑150 (MQ=255)
                                                                                                            ctctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCaga                                           <  2:4361777/150‑1 (MQ=255)
                                                                                                            ctctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCaga                                           <  2:2196949/150‑1 (MQ=255)
                                                                                                              ctTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAgaga                                         >  1:2733771/1‑150 (MQ=255)
                                                                                                                 cAAGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCa                                      >  2:3263754/1‑150 (MQ=255)
                                                                                                                  aaGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAg                                     <  2:1257849/150‑1 (MQ=255)
                                                                                                                  aaGCTGCTGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAg                                     >  2:4005402/1‑150 (MQ=255)
                                                                                                                    gctgctGGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGat                                   <  2:1486888/150‑1 (MQ=255)
                                                                                                                      tgctgGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGataa                                 >  1:3972831/1‑150 (MQ=255)
                                                                                                                      tgctgGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGataa                                 >  2:1317067/1‑150 (MQ=255)
                                                                                                                        ctgGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGataata                               >  2:3578351/1‑150 (MQ=255)
                                                                                                                        ctgGCGATGATAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGataata                               >  2:3551309/1‑150 (MQ=255)
                                                                                                                         tgGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAatat                              >  2:1423576/1‑150 (MQ=255)
                                                                                                                         tgGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAatat                              >  2:1423643/1‑150 (MQ=255)
                                                                                                                          gcctatgatAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATAtt                             <  2:2511292/146‑1 (MQ=255)
                                                                                                                            cGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTa                           >  1:3191881/1‑150 (MQ=255)
                                                                                                                                gatAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCtt                       >  2:2717521/1‑150 (MQ=255)
                                                                                                                                gatAATCATCCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCtt                       <  2:636297/150‑1 (MQ=255)
                                                                                                                                    atcatcCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTccca                   >  1:2611014/1‑150 (MQ=255)
                                                                                                                                     tcatcCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccag                  >  1:4233073/1‑150 (MQ=255)
                                                                                                                                       atcCCACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcc                <  2:1178922/150‑1 (MQ=255)
                                                                                                                                         cccACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagccca              <  2:3850087/150‑1 (MQ=255)
                                                                                                                                         cccACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagccca              >  2:1830545/1‑150 (MQ=255)
                                                                                                                                          ccACCTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcccag             <  1:3645277/150‑1 (MQ=255)
                                                                                                                                          ccACATGAATTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTcccagcccag             >  2:1548024/1‑150 (MQ=255)
                                                                                                                                           cacaTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCccagcccagc            >  2:2513947/1‑150 (MQ=255)
                                                                                                                                              aTGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGcc         >  2:1546049/1‑150 (MQ=255)
                                                                                                                                                gatttCTACTTTCTCTGGATGTTTCGCTTCTACAGCCTGCACGACCGGGTTTTCACTGATTACGGTCATCACATCCCCGCAGAGATCTTTCACGCCCAGTCCTCTCGCCGCAGCGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTc       <  2:1616411/147‑1 (MQ=255)
                                                                                                                                                gAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTc       <  2:868236/150‑1 (MQ=255)
                                                                                                                                                 aaTTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCa      >  1:1048693/1‑150 (MQ=255)
                                                                                                                                                 aaTTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCa      >  1:2543241/1‑150 (MQ=255)
                                                                                                                                                  aTTAGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAg     >  1:4395941/1‑150 (MQ=255)
                                                                                                                                                   ttCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAGc    <  2:2936224/150‑1 (MQ=255)
                                                                                                                                                     cGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAGCGa  >  1:1520700/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCGTCGTCTTCGTCCCAGTCGTCATCCCAGTCTTCATCATCCTCTTCAGCAATCTCTTCAAGCTGCTGGCGATGATAATCATCCCACATGAATTCGACTTTCTCTGGCTGTTTCGCTTCTTCAGCCTGCACGACCGGGTTTTCAATGATAAAGGTCATCACATCCCAGCAGAGATCTTTCACGCCCAGTCCACTCGCCGCAGAGATCAGATAATATTTATCTTCCCAGCCCAGCGCCTCAGCGA  >  CP009273/3323917‑3324215

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: