Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,306,969 T→A 10.8% N23Y (AAC→TAC)  phnM ← ribophosphonate triphosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,306,9690TA10.8% 797.2 / 75.7 249N23Y (AAC→TAC) phnMribophosphonate triphosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (88/134);  new base A (10/17);  total (98/151)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 250 total aligned reads displayed.

AGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGC  >  CP009273/4306820‑4307112
                                                                                                                                                     |                                                                                                                                               
 ggCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTc                                                                                                                                                <  1:858727/150‑1 (MQ=255)
  gCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCt                                                                                                                                               <  1:4878709/150‑1 (MQ=255)
  gCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACt                                                                                                                                               <  1:1492209/150‑1 (MQ=255)
   ccAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTg                                                                                                                                              <  1:2380142/150‑1 (MQ=255)
      gTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAc                                                                                                                                           >  2:1350612/1‑150 (MQ=255)
      gTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAc                                                                                                                                           <  1:225608/150‑1 (MQ=255)
       tCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCAcc                                                                                                                                          <  1:1267410/150‑1 (MQ=255)
         aaCTTTCGGGCGCGGGGTGAAGAATTTCTCCAGATGATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTc                                                                                                                                        <  2:459086/150‑1 (MQ=255)
               cGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGa                                                                                                                                  >  2:2438913/1‑150 (MQ=255)
                  gcgcGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAcc                                                                                                                               <  1:1334674/150‑1 (MQ=255)
                   cgcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                              >  2:4845931/1‑150 (MQ=255)
                   cgcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                              <  1:4336269/150‑1 (MQ=255)
                   cgcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                              <  2:3945929/150‑1 (MQ=255)
                   cgcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCATTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                              <  1:2224507/150‑1 (MQ=255)
                        ggTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTac                                                                                                                         <  2:3009275/150‑1 (MQ=255)
                           gaagaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTaccac                                                                                                                      >  1:3212286/1‑150 (MQ=255)
                             agaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCt                                                                                                                    >  1:1123769/1‑150 (MQ=255)
                               aaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCg                                                                                                                  <  2:1063393/150‑1 (MQ=255)
                                aTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGt                                                                                                                 >  2:4978183/1‑150 (MQ=255)
                                  ttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                                                               <  1:2502834/150‑1 (MQ=255)
                                    aTCCAGATTATCGGTATGCAGCTCAATCATCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTc                                                                                                             >  2:4943344/1‑150 (MQ=255)
                                    aTCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTc                                                                                                             <  1:2682438/150‑1 (MQ=255)
                                      ccAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCa                                                                                                           >  1:2316211/1‑150 (MQ=255)
                                      ccAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCa                                                                                                           >  1:1086164/1‑150 (MQ=255)
                                        aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGc                                                                                                         >  2:1931640/1‑150 (MQ=255)
                                        aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGc                                                                                                         >  2:1931785/1‑150 (MQ=255)
                                         gATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa                                                                                                        <  2:1743627/150‑1 (MQ=255)
                                          aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAc                                                                                                       <  1:3373257/150‑1 (MQ=255)
                                           ttATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAcc                                                                                                      <  1:4692839/150‑1 (MQ=255)
                                            tATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCa                                                                                                     >  2:1897917/1‑150 (MQ=255)
                                            tATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCa                                                                                                     <  2:4185131/150‑1 (MQ=255)
                                                  tATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTaa                                                                                               <  2:3546520/150‑1 (MQ=255)
                                                  tATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTaa                                                                                               >  2:3958343/1‑150 (MQ=255)
                                                    tGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg                                                                                             <  1:680050/150‑1 (MQ=255)
                                                    tGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg                                                                                             <  2:3664165/150‑1 (MQ=255)
                                                     gCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGt                                                                                            >  1:1481187/1‑150 (MQ=255)
                                                     gCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGt                                                                                            <  2:601266/150‑1 (MQ=255)
                                                      cAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGtt                                                                                           <  2:2399861/150‑1 (MQ=255)
                                                       aGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                                                                                          >  2:2984967/1‑150 (MQ=255)
                                                       aGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                                                                                          >  1:2353304/1‑150 (MQ=255)
                                                        gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTAt                                                                                         <  2:4715382/150‑1 (MQ=255)
                                                              tCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCCGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATaa                                                                                   >  1:2856445/1‑150 (MQ=255)
                                                               cAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAAt                                                                                  <  1:5036878/150‑1 (MQ=255)
                                                               cAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAAt                                                                                  <  2:201140/150‑1 (MQ=255)
                                                               cAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGAGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAAt                                                                                  >  2:2863618/1‑150 (MQ=255)
                                                                aGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATc                                                                                 <  2:1188239/150‑1 (MQ=255)
                                                                 gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCa                                                                                <  1:1999933/150‑1 (MQ=255)
                                                                 gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCa                                                                                <  2:968857/150‑1 (MQ=255)
                                                                 gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCa                                                                                <  1:4135610/150‑1 (MQ=255)
                                                                 gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTGAACGGTCGTGATAATCa                                                                                >  2:3449021/1‑150 (MQ=255)
                                                                     cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGaa                                                                            >  2:5119005/1‑150 (MQ=255)
                                                                      ggCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAag                                                                           <  2:2314188/150‑1 (MQ=255)
                                                                        cagcagCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAagag                                                                         <  2:3025430/150‑1 (MQ=255)
                                                                        cagcagCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAagag                                                                         <  1:3058614/150‑1 (MQ=255)
                                                                          gcagcCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGc                                                                       >  2:1658939/1‑150 (MQ=255)
                                                                              ccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTccc                                                                   <  2:4665348/150‑1 (MQ=255)
                                                                                agccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCAt                                                                 >  1:1851971/1‑150 (MQ=255)
                                                                                agccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTAGTGATAATCATGAAGAGGCTCCCAt                                                                 >  2:3103945/1‑150 (MQ=255)
                                                                                   cgccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTgg                                                              <  1:3732230/150‑1 (MQ=255)
                                                                                    gccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTggg                                                             >  2:3729859/1‑150 (MQ=255)
                                                                                    gccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTggg                                                             <  2:611945/150‑1 (MQ=255)
                                                                                     ccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGt                                                            <  1:3407633/150‑1 (MQ=255)
                                                                                       ttCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGc                                                          <  1:5132198/150‑1 (MQ=255)
                                                                                       ttCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGc                                                          >  1:3975942/1‑150 (MQ=255)
                                                                                        tCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCa                                                         <  1:833999/150‑1 (MQ=255)
                                                                                          gCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAgg                                                       <  1:1578672/150‑1 (MQ=255)
                                                                                           ccGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTGAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGc                                                      >  2:1643960/1‑150 (MQ=255)
                                                                                             gTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCgg                                                    <  2:1015804/150‑1 (MQ=255)
                                                                                              tCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGt                                                   <  1:1653754/150‑1 (MQ=255)
                                                                                               ccATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTGAACGTTATTGATAATCATGAAGAGGCTCCCAGTGGGTGCAGGCGGTc                                                  >  2:4376241/1‑150 (MQ=255)
                                                                                               ccATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACAGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGATGAGGCTCCCATTGGGTGCAGGCGTTc                                                  <  1:4187152/150‑1 (MQ=255)
                                                                                                cATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCg                                                 >  1:2890881/1‑150 (MQ=255)
                                                                                                     ccTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGAc                                            <  2:3498428/150‑1 (MQ=255)
                                                                                                      cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACg                                           <  1:1033296/150‑1 (MQ=255)
                                                                                                      cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACg                                           >  2:4163890/1‑150 (MQ=255)
                                                                                                         cGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATTATCATGAAGATGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCa                                        <  1:4914846/150‑1 (MQ=255)
                                                                                                          ggCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCAt                                       >  1:3637279/1‑150 (MQ=255)
                                                                                                             aGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTAc                                    <  1:3237786/150‑1 (MQ=255)
                                                                                                                 ggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTAc                                <  1:3512678/150‑1 (MQ=255)
                                                                                                                 ggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTAc                                <  1:1012454/150‑1 (MQ=255)
                                                                                                                  gcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACa                               >  2:275800/1‑150 (MQ=255)
                                                                                                                    tctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGc                             <  2:225091/150‑1 (MQ=255)
                                                                                                                     ctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCt                            >  2:1124350/1‑150 (MQ=255)
                                                                                                                              gtcAAAGGCGCGGATTTCGCCGTTCTTCACCTCCAGCGAACCGCTTACCTCCTCGTTTTCCAGCACCAGCTTAACGTTTTTGTTAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGGGACGTCATTACGTACAGCTTCGTCATgg                   <  1:4385711/148‑1 (MQ=255)
                                                                                                                              ggCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATgg                   >  1:5173882/1‑150 (MQ=255)
                                                                                                                              ggCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATgg                   >  2:2475063/1‑150 (MQ=255)
                                                                                                                              ggCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATgg                   >  1:5173892/1‑150 (MQ=255)
                                                                                                                                 aaaGGCGCGGATATCGCCGTTCTGCACCTCCAGCCAACCGCTTACCACCTCGTTTACCAGCACCAGCATAACGGTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGGCAATACGAACAGCTTCGTCATGGAAg                >  2:4509193/1‑150 (MQ=255)
                                                                                                                                   aGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGAt              >  1:391652/1‑150 (MQ=255)
                                                                                                                                   aGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGAt              >  2:785501/1‑150 (MQ=255)
                                                                                                                                   aGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGAt              <  1:582470/150‑1 (MQ=255)
                                                                                                                                   aGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGAt              <  1:2754879/150‑1 (MQ=255)
                                                                                                                                     gcgcGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGc            <  1:4678849/150‑1 (MQ=255)
                                                                                                                                          gATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACg       <  1:4642110/150‑1 (MQ=255)
                                                                                                                                          gATTTCGCCGTTCTGCACCTCCAGCGAACCGCTGACCACCTCGTTTTCCAGCCCCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCAGTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACg       <  2:4866526/150‑1 (MQ=255)
                                                                                                                                          gATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTaag       <  1:311606/150‑3 (MQ=255)
                                                                                                                                             ttCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATg    <  1:4037018/150‑1 (MQ=255)
                                                                                                                                              tCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATgg   <  2:635184/150‑1 (MQ=255)
                                                                                                                                               cGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGc  <  1:4534233/150‑1 (MQ=255)
                                                                                                                                               cGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGc  <  2:2349797/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGC  >  CP009273/4306820‑4307112

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: