Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,889,739 C→A 8.9% M77I (ATG→ATT purP ← adenine permease, high affinity; adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,889,7390CA 8.9% 582.6 / 24.4 168M77I (ATG→ATTpurPadenine permease, high affinity; adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (81/72);  new base A (7/8);  total (88/80)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 175 total aligned reads displayed.

TAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCGGTCCG  >  CP009273/3889592‑3889885
                                                                                                                                                   |                                                                                                                                                  
tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatc                                                                                                                                                    >  1:272031/1‑148 (MQ=255)
     gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataat                                                                                                                                               <  2:1247276/148‑1 (MQ=255)
     gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataat                                                                                                                                               <  2:1922977/148‑1 (MQ=255)
      cagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataata                                                                                                                                              <  2:1269569/148‑1 (MQ=255)
       agACCTATCGCGCCCCAGAAGAGTGCGGCCATCGCGACCTGCCATTGCAAGCCCATCTGCTGTACATCGACAAAAGCGAAGAACGCATTCATGCCCATAGCGGGTTCCAGTGCTACTGGCTGGTTAGCAAACTGTCCCATCATTATAc                                                                                                                                             <  2:2895884/148‑1 (MQ=255)
         aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTg                                                                                                                                           >  2:593477/1‑148 (MQ=255)
         aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTg                                                                                                                                           >  2:2255950/1‑148 (MQ=255)
         aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTg                                                                                                                                           <  2:1625884/148‑1 (MQ=255)
         aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATAATAATACTg                                                                                                                                           >  1:2530412/1‑148 (MQ=255)
            tATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCg                                                                                                                                        >  1:891478/1‑148 (MQ=255)
              tCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGaa                                                                                                                                      <  1:3027603/148‑1 (MQ=255)
                 cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGc                                                                                                                                   >  1:396538/1‑148 (MQ=255)
                 cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGc                                                                                                                                   <  2:1043068/148‑1 (MQ=255)
                        gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                                                                                            >  2:2030996/1‑148 (MQ=255)
                         aagaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                                                                                           >  2:730693/1‑148 (MQ=255)
                           gaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGa                                                                                                                         >  1:261863/1‑148 (MQ=255)
                            aTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGAc                                                                                                                        <  2:232039/148‑1 (MQ=255)
                             ttGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACa                                                                                                                       <  1:1736748/148‑1 (MQ=255)
                             ttGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACa                                                                                                                       >  2:599368/1‑148 (MQ=255)
                              tGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                                                                                                      <  2:3183676/148‑1 (MQ=255)
                                   ccATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGt                                                                                                                 >  1:2661660/1‑148 (MQ=255)
                                       cccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACg                                                                                                             >  1:1552106/1‑148 (MQ=255)
                                       cccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACg                                                                                                             >  1:3109919/1‑148 (MQ=255)
                                        ccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGa                                                                                                            >  2:1581322/1‑148 (MQ=255)
                                        ccGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGa                                                                                                            <  2:1750733/148‑1 (MQ=255)
                                         cGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGaa                                                                                                           >  1:422807/1‑148 (MQ=255)
                                         cGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGaa                                                                                                           <  2:2821866/148‑1 (MQ=255)
                                         cGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGaa                                                                                                           >  1:3075201/1‑148 (MQ=255)
                                           aCCAGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGAAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGACGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                                                                                         >  2:985995/1‑148 (MQ=255)
                                             gtGCCACGGCAAGCCCATCGCCTGTACAACGACAACAGCGAAGAACGCATTCAGGGCCATAGCGGGTGCCTGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATCCTGACGAATGCAGCGATCAGACAGGTAGTGACGAAGACg                                                                                                       <  1:1087561/147‑1 (MQ=255)
                                              tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgg                                                                                                      <  2:1232646/148‑1 (MQ=255)
                                                     gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGt                                                                                               >  1:2759804/1‑148 (MQ=255)
                                                      cAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTa                                                                                              >  2:3306783/1‑148 (MQ=255)
                                                        aGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATc                                                                                            >  1:1802478/1‑148 (MQ=255)
                                                        aGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATc                                                                                            <  1:1919711/148‑1 (MQ=255)
                                                         gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATcc                                                                                           <  1:1482135/148‑1 (MQ=255)
                                                         gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATcc                                                                                           >  1:1816411/1‑148 (MQ=255)
                                                         gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATcc                                                                                           >  1:1640810/1‑148 (MQ=255)
                                                          cccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCa                                                                                          <  2:2947230/148‑1 (MQ=255)
                                                          cccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCa                                                                                          >  2:774664/1‑148 (MQ=255)
                                                             aTCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGc                                                                                       >  1:1924526/1‑148 (MQ=255)
                                                             aTCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGc                                                                                       <  1:240104/148‑1 (MQ=255)
                                                               cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCa                                                                                     >  2:815637/1‑148 (MQ=255)
                                                               cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCa                                                                                     >  1:2481104/1‑148 (MQ=255)
                                                                 ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGc                                                                                   >  1:1126409/1‑148 (MQ=255)
                                                                  cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATCGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCa                                                                                  <  2:482144/148‑1 (MQ=255)
                                                                  cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCa                                                                                  <  1:1340073/148‑1 (MQ=255)
                                                                   tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCaa                                                                                 >  2:519988/1‑148 (MQ=255)
                                                                      aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGc                                                                              <  1:1833444/148‑1 (MQ=255)
                                                                       cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGcc                                                                             >  2:2220702/1‑148 (MQ=255)
                                                                       cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGcc                                                                             >  1:2899955/1‑148 (MQ=255)
                                                                              aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAAtt                                                                      <  2:3031295/148‑1 (MQ=255)
                                                                               aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAAttt                                                                     >  2:2577546/1‑148 (MQ=255)
                                                                                  gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                                                  >  1:1616742/1‑148 (MQ=255)
                                                                                     aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCTCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGt                                                               <  2:931772/148‑1 (MQ=255)
                                                                                        aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGTCAAGAATTTGCGGGTTaa                                                            >  2:2568860/1‑148 (MQ=255)
                                                                                        aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTaa                                                            <  2:2450833/148‑1 (MQ=255)
                                                                                             aTTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaa                                                       <  1:2944027/148‑1 (MQ=255)
                                                                                               tCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAAc                                                     >  2:2917824/1‑148 (MQ=255)
                                                                                               tCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAAc                                                     >  2:1723972/1‑148 (MQ=255)
                                                                                               tCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAAc                                                     >  1:3149806/1‑148 (MQ=255)
                                                                                                    cccATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                                                <  1:264566/148‑1 (MQ=255)
                                                                                                     ccATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTa                                               <  2:3267562/148‑1 (MQ=255)
                                                                                                      cATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaa                                              >  2:1891820/1‑148 (MQ=255)
                                                                                                      cATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaa                                              <  1:2761229/148‑1 (MQ=255)
                                                                                                       aTAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaaa                                             <  1:2868514/148‑1 (MQ=255)
                                                                                                       aTAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaaa                                             <  1:1074389/148‑1 (MQ=255)
                                                                                                       aTAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaaa                                             <  1:1074391/148‑1 (MQ=255)
                                                                                                        tAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                            >  2:374990/1‑148 (MQ=255)
                                                                                                        tAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                            >  2:549095/1‑148 (MQ=255)
                                                                                                         aGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAcc                                           <  1:2969684/148‑1 (MQ=255)
                                                                                                          gCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCa                                          >  2:522526/1‑148 (MQ=255)
                                                                                                           cGGGTGCCAGTGCAACTGGCAGTTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAATACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                                         >  2:792131/1‑148 (MQ=255)
                                                                                                           cGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                                         <  2:3083000/148‑1 (MQ=255)
                                                                                                           cGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                                         >  1:2545115/1‑148 (MQ=255)
                                                                                                           cGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                                         <  2:1906005/148‑1 (MQ=255)
                                                                                                              gTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt                                      <  2:369193/148‑1 (MQ=255)
                                                                                                                      gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAgg                              >  1:2224320/1‑148 (MQ=255)
                                                                                                                       cAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggt                             <  1:3205610/148‑1 (MQ=255)
                                                                                                                       cAACTGGCAGGTTAGCAAACAGTCCCATAATAATNCTGCCGAATGCAGCGATCAGACAGGNAGTGACGAAGACGGCGCTGGNATNCATGCCAGCAACGCCAAGAATTTGNGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggt                             <  1:738289/148‑1 (MQ=255)
                                                                                                                             gCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaa                       <  1:2040391/148‑1 (MQ=255)
                                                                                                                              cAGGTTAGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaaa                      <  2:2530412/148‑1 (MQ=255)
                                                                                                                                  ttAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCgg                  >  2:632959/1‑148 (MQ=255)
                                                                                                                                  ttAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCgg                  <  2:1649138/148‑1 (MQ=255)
                                                                                                                                   tAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGc                 >  1:3256030/1‑148 (MQ=255)
                                                                                                                                   tAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGc                 >  1:1656205/1‑148 (MQ=255)
                                                                                                                                     gCAAACAGTCCCATCATAATACTGTCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGa               <  2:2121688/148‑1 (MQ=255)
                                                                                                                                     gCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGa               >  1:2482744/1‑148 (MQ=255)
                                                                                                                                       aaaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATc             <  1:577071/148‑1 (MQ=255)
                                                                                                                                        aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCa            >  1:1320105/1‑148 (MQ=255)
                                                                                                                                        aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCa            >  2:526825/1‑148 (MQ=255)
                                                                                                                                        aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCa            <  1:878889/148‑1 (MQ=255)
                                                                                                                                         aCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCAc           >  2:2679939/1‑148 (MQ=255)
                                                                                                                                          cAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACt          >  1:2007236/1‑148 (MQ=255)
                                                                                                                                           aGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACtt         >  2:2496360/1‑148 (MQ=255)
                                                                                                                                             tCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCg       >  1:2306272/1‑148 (MQ=255)
                                                                                                                                              cccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCgg      >  1:1948024/1‑148 (MQ=255)
                                                                                                                                                catcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCGGTc    >  1:1920495/1‑148 (MQ=255)
                                                                                                                                                 atcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCGGTcc   >  2:1454221/1‑148 (MQ=255)
                                                                                                                                                  tcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCGGTCCg  >  2:2641974/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATCACTTCGGTCCG  >  CP009273/3889592‑3889885

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: