Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 327,051 A→T 18.0% E626D (GAA→GAT betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273327,0510AT18.0% 501.6 / 104.2 211E626D (GAA→GATbetTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base A (87/86);  new base T (27/11);  total (114/97)
Fisher's exact test for biased strand distribution p-value = 3.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 219 total aligned reads displayed.

AGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCG  >  CP009273/326904‑327198
                                                                                                                                                   |                                                                                                                                                   
agagGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGaa                                                                                                                                                     >  1:127149/1‑148 (MQ=255)
agagGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGaa                                                                                                                                                     >  1:2058137/1‑148 (MQ=255)
  agGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAAc                                                                                                                                                   <  2:2717294/148‑1 (MQ=255)
  agGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAAc                                                                                                                                                   <  2:190261/148‑1 (MQ=255)
     gACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCtt                                                                                                                                                >  1:117200/1‑148 (MQ=255)
      aCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTc                                                                                                                                               >  1:509292/1‑148 (MQ=255)
      aCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTc                                                                                                                                               <  1:872878/148‑1 (MQ=255)
       cagcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTcc                                                                                                                                              <  2:1237575/148‑1 (MQ=255)
        agcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCt                                                                                                                                             >  2:2315492/1‑148 (MQ=255)
          cagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGt                                                                                                                                           <  1:3552929/148‑1 (MQ=255)
           agTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGtt                                                                                                                                          <  2:891987/148‑1 (MQ=255)
               gggTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGaa                                                                                                                                      <  1:1578061/148‑1 (MQ=255)
                 gTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTATAAgg                                                                                                                                    >  1:3050301/1‑148 (MQ=255)
                       tGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCa                                                                                                                              <  1:2313178/148‑1 (MQ=255)
                         gATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAgg                                                                                                                            >  2:1095326/1‑148 (MQ=255)
                          aTTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAggg                                                                                                                           >  1:2767630/1‑148 (MQ=255)
                          aTTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAggg                                                                                                                           >  1:1905230/1‑148 (MQ=255)
                                tgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGa                                                                                                                     <  2:171184/148‑1 (MQ=255)
                                 ggTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTNCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGAc                                                                                                                    <  1:1937082/148‑1 (MQ=255)
                                       tatGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATg                                                                                                              >  1:2298396/1‑148 (MQ=255)
                                        atGGGCGAAGAGCAAAACTTTTTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCATCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCATGGCAACGACCTGATgg                                                                                                             >  2:2826735/1‑148 (MQ=255)
                                         tGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGa                                                                                                            <  2:3760007/148‑1 (MQ=255)
                                          gggCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGAc                                                                                                           >  1:2096441/1‑148 (MQ=255)
                                          gggCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGAc                                                                                                           <  2:203648/148‑1 (MQ=255)
                                          gggCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATAGTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGAc                                                                                                           <  2:203595/148‑1 (MQ=255)
                                           ggCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACt                                                                                                          <  1:2920672/148‑1 (MQ=255)
                                           ggCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACt                                                                                                          >  1:2617639/1‑148 (MQ=255)
                                            gCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTCCCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTa                                                                                                         <  2:2490881/148‑1 (MQ=255)
                                            gCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACGACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTa                                                                                                         <  2:570599/148‑1 (MQ=255)
                                             cGAAGAGCAAAACTTTGTCTATCAGATTTGGACGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTAc                                                                                                        <  1:686405/148‑1 (MQ=255)
                                               aaGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAg                                                                                                      >  1:3438750/1‑148 (MQ=255)
                                                agagCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGc                                                                                                     >  2:401724/1‑148 (MQ=255)
                                                  agCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCaa                                                                                                   >  1:10001/1‑148 (MQ=255)
                                                    cAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAag                                                                                                 <  2:577763/148‑1 (MQ=255)
                                                      aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAagag                                                                                               >  1:1413019/1‑148 (MQ=255)
                                                      aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAagag                                                                                               >  1:1020713/1‑148 (MQ=255)
                                                       aaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGc                                                                                              >  1:611091/1‑148 (MQ=255)
                                                          tttGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg                                                                                           <  2:3796894/148‑1 (MQ=255)
                                                          tttGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg                                                                                           <  2:3594773/148‑1 (MQ=255)
                                                          tttGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg                                                                                           <  2:3594106/148‑1 (MQ=255)
                                                          tttGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg                                                                                           <  2:3486701/148‑1 (MQ=255)
                                                          tttGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg                                                                                           >  2:2069591/1‑148 (MQ=255)
                                                              tCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGAt                                                                                       <  1:396260/148‑1 (MQ=255)
                                                              tCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGAt                                                                                       <  1:2909963/148‑1 (MQ=255)
                                                              tCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGAt                                                                                       >  2:291056/1‑148 (MQ=255)
                                                                 aTCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCAc                                                                                    <  1:1526547/148‑1 (MQ=255)
                                                                 aTCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCAc                                                                                    <  1:1017644/148‑1 (MQ=255)
                                                                   cAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCg                                                                                  >  1:1055520/1‑148 (MQ=255)
                                                                     gATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGat                                                                                <  1:3120882/148‑1 (MQ=255)
                                                                      aTTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGata                                                                               <  1:1923729/148‑1 (MQ=255)
                                                                      aTTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGata                                                                               >  1:3724873/1‑148 (MQ=255)
                                                                        ttGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATAtt                                                                             <  2:1463876/148‑1 (MQ=255)
                                                                         tGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTc                                                                            >  1:44625/1‑148 (MQ=255)
                                                                         tGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTc                                                                            <  2:627126/148‑1 (MQ=255)
                                                                          ggCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCt                                                                           <  1:1733809/148‑1 (MQ=255)
                                                                          ggCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCt                                                                           <  1:969776/148‑1 (MQ=255)
                                                                             cGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGa                                                                        <  1:547558/148‑1 (MQ=255)
                                                                               cagcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAcc                                                                      >  1:2092215/1‑148 (MQ=255)
                                                                               cagcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAcc                                                                      >  2:651379/1‑148 (MQ=255)
                                                                                agcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCa                                                                     <  2:3795276/148‑1 (MQ=255)
                                                                                 gcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAg                                                                    >  1:1141568/1‑148 (MQ=255)
                                                                                 gcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAg                                                                    >  1:120222/1‑148 (MQ=255)
                                                                                  caATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                   >  1:3318190/1‑148 (MQ=255)
                                                                                  caATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                   >  2:1730500/1‑148 (MQ=255)
                                                                                  caATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                   >  1:639840/1‑148 (MQ=255)
                                                                                  caATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                   >  1:1417835/1‑148 (MQ=255)
                                                                                   aaTATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTa                                                                  <  2:3679656/148‑1 (MQ=255)
                                                                                    atatTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTAc                                                                 >  1:225869/1‑148 (MQ=255)
                                                                                     tatTCGGTGCCGGGCTTTACCGACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACg                                                                <  1:3147408/148‑1 (MQ=255)
                                                                                      atTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGa                                                               <  1:2661228/148‑1 (MQ=255)
                                                                                       ttCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACTACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAg                                                              >  1:2175254/1‑148 (MQ=255)
                                                                                       ttCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAg                                                              >  1:3517230/1‑148 (MQ=255)
                                                                                        tCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGc                                                             >  2:1959079/1‑148 (MQ=255)
                                                                                            tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCa                                                         <  1:2123865/148‑1 (MQ=255)
                                                                                             gCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCAc                                                        >  1:3239924/1‑148 (MQ=255)
                                                                                                   cTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAAc                                                  >  1:547821/1‑148 (MQ=255)
                                                                                                      tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACttt                                               <  1:1976563/148‑1 (MQ=255)
                                                                                                         ctaccGCGCACGCAGCGGTAAATCGACCTACTACCGACTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTAtt                                            >  1:1772112/1‑148 (MQ=255)
                                                                                                          taccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCACCGACCTGATGGACTACAGCACAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTc                                           >  2:2822578/1‑148 (MQ=255)
                                                                                                          taccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTc                                           >  2:3456810/1‑148 (MQ=255)
                                                                                                          taccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTc                                           >  1:1458270/1‑148 (MQ=255)
                                                                                                            ccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCAt                                         >  1:318347/1‑148 (MQ=255)
                                                                                                                  acgcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCAt                                   <  2:385765/148‑1 (MQ=255)
                                                                                                                  acgcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCAt                                   <  2:299415/148‑1 (MQ=255)
                                                                                                                   cgcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATc                                  >  2:2326080/1‑148 (MQ=255)
                                                                                                                    gcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCg                                 >  1:1856695/1‑148 (MQ=255)
                                                                                                                    gcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCg                                 <  2:1367962/148‑1 (MQ=255)
                                                                                                                     caGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGt                                >  1:1607040/1‑148 (MQ=255)
                                                                                                                     caGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGt                                >  1:2432925/1‑148 (MQ=255)
                                                                                                                              aTCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCg                       >  1:3299687/1‑148 (MQ=255)
                                                                                                                              aTCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCg                       >  2:2187633/1‑148 (MQ=255)
                                                                                                                                     tactacCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATa                >  1:3584260/1‑148 (MQ=255)
                                                                                                                                     tactacCGGCTGGATACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATa                >  1:1688338/1‑148 (MQ=255)
                                                                                                                                      actacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAg               <  2:3318190/148‑1 (MQ=255)
                                                                                                                                         acCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACGGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGt            >  1:3231441/1‑148 (MQ=255)
                                                                                                                                         acCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGt            >  1:2944250/1‑148 (MQ=255)
                                                                                                                                              cTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGt       >  1:3497662/1‑148 (MQ=255)
                                                                                                                                                ggAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGTCCAGTACGAGCGGCACCTTAACTTTAGTCATCTCCATCGTGTAGCGCCGGGCCATAGCGTGTTGTTc     >  2:1972975/1‑148 (MQ=255)
                                                                                                                                                ggAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTc     <  1:2315492/148‑1 (MQ=255)
                                                                                                                                                   aaCCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCg  >  1:1750054/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCG  >  CP009273/326904‑327198

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: