Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 326,966 T→C 4.9% V598A (GTC→GCC)  betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273326,9660TC 4.9% 243.9 / 2.0 89V598A (GTC→GCC) betTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base C (2/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATC  >  CP009273/326821‑327114
                                                                                                                                                 |                                                                                                                                                    
gTCTGGTACCCGGCAATGGGAGAAGTGGCGCCGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACCGCACTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTa                                                                                                                                                   <  1:44026/149‑1 (MQ=255)
       aCCCGGCAATGGAAGAAGTGGCGCAGGGGGTGCGGTTGCGCGTCGCGTACGTGGAGCTAAAAAGACCGCCACCGGCAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTCGTCTATCAGAtt                                                                                                                                            >  1:572148/1‑149 (MQ=255)
        cccGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAGAGCTTTGTCTATc                                                                                                                                                 >  2:1080254/1‑143 (MQ=255)
        cccGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGAAAGAGGGACAACAGTTGAGTCATCTGGATTTGTTGGGGCATATAGGCGAAGAGCAAAACTTTGTCTATc                                                                                                                                                 >  2:304649/1‑143 (MQ=255)
         ccGGCAATGGAAGAAGTGGCGCAGGTGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTg                                                                                                                                          >  1:652667/1‑149 (MQ=255)
           ggCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGCTTCGTTGGTGCATATGGGCGAAGAGCAAAACTTGGTCTATCCGa                                                                                                                                              >  2:431555/1‑143 (MQ=255)
                  gaagaaGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGCAGCGCAAAACTTTGTCTATCAGATTTGGCCgcagca                                                                                                                                 >  1:267102/1‑149 (MQ=255)
                  gaagaaGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCgcagca                                                                                                                                 >  1:299167/1‑149 (MQ=255)
                   aagaagTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGCCTATCAGATTTGGCCGCAGCaa                                                                                                                                >  1:68906/1‑149 (MQ=255)
                   aagaagTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGCTAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGCCTATCAGATTTGGCCGCAGCaa                                                                                                                                >  1:872558/1‑149 (MQ=255)
                     gaagGGGCGCAGGGGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCgca                                                                                                                                    <  2:720123/143‑1 (MQ=255)
                      aagTGGCGTAGGCGTTGCGGTTGCGCGGCGCGTACGTGGAGCCAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCgcag                                                                                                                                   <  2:890960/143‑1 (MQ=255)
                        ggggggcAGGGGGTGCGGGTGCGCGGGGCGTACGTGGGGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCgcagca                                                                                                                                 <  2:218934/138‑1 (MQ=255)
                        ggggCGCAGGGGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCgcagca                                                                                                                                 <  2:720711/141‑1 (MQ=255)
                         tGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCGGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATAGTCg                                                                                                                          >  1:1278510/1‑149 (MQ=255)
                          ggCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTCGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAat                                                                                                                               >  2:621618/1‑143 (MQ=255)
                            cgcATGATTTGCGGTTGCGCGGCGCGTACGCGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTCGGGTCATCTGGATTTGTAGGTGCATATGGGCGAAGAGCAACACTTTGTCTATCAGATTTGGCCGCAGTAATATTCGGTg                                                                                                                       <  1:95931/149‑1 (MQ=255)
                             gcAGGAGTTGCGGTTGCGCGGCGCGTCCCCAGCGGTTATAGGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGc                                                                                                                      >  1:1286831/1‑149 (MQ=255)
                             gcAGGAGTTGCGGTAGCGCGGCGCCTACGTGGAGCTAAAACGCCTCCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATAtt                                                                                                                            <  2:967327/143‑1 (MQ=255)
                                   gttgcggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGc                                                                                                                      >  2:843940/1‑143 (MQ=255)
                                   gttgcggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTCGAGTCATCTGGATTTGTTGGTGCAGATGGGCGAAGAGCACACCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGtcc                                                                                                                      >  2:583608/1‑141 (MQ=255)
                                    ttgcggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCtt                                                                                                               <  1:648238/149‑1 (MQ=255)
                                     tgcggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCg                                                                                                                    <  2:687655/143‑1 (MQ=255)
                                        ggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACGGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGc                                                                                                                 <  2:1039475/143‑1 (MQ=255)
                                        ggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTtacc                                                                                                           >  1:719400/1‑149 (MQ=255)
                                            gcgcggcgggTACGTGGGCGTAAAAGGCCTGCCCCCGGAAGGGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTta                                                                                                             <  2:1013139/143‑1 (MQ=255)
                                               cggcgcgTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTAGCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACcgcg                                                                                                    >  1:1162553/1‑149 (MQ=255)
                                                  cgcgTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTtacctacc                                                                                                       >  2:339384/1‑143 (MQ=255)
                                                    cgTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGcgcacgc                                                                                               >  1:112459/1‑149 (MQ=255)
                                                      tACCTGGAGCTAAACAGCCTGCCCCCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAg                                                                                             <  1:169596/149‑1 (MQ=255)
                                                       aCGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGcgca                                                                                                  <  2:847218/143‑1 (MQ=255)
                                                         gTGGAGCTAAAAAGCCTGCCACCGGAAGCGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCgg                                                                                          <  1:132494/149‑1 (MQ=255)
                                                         gTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCGTCTGGGTTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTAGCTGACGCGCGCGCAGCgg                                                                                          >  1:1121903/1‑149 (MQ=255)
                                                             aGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTaaa                                                                                      <  1:1193810/149‑1 (MQ=255)
                                                             aGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTaaa                                                                                      <  1:147889/149‑1 (MQ=255)
                                                               cTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATc                                                                                    >  1:216541/1‑149 (MQ=255)
                                                                   agaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCCATCAGCTTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGAcc                                                                                <  1:368828/147‑1 (MQ=255)
                                                                   acaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTaaa                                                                                      <  2:867868/141‑1 (MQ=255)
                                                                    aaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCGGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAAt                                                                                     >  2:7536/1‑143 (MQ=255)
                                                                      gCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCg                                                                                   <  2:585078/143‑1 (MQ=255)
                                                                        atGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACctacta                                                                           <  1:547819/148‑1 (MQ=255)
                                                                         tGCCGCCGGAAGAGGGACAGCGGTTGGGTTATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAACCTTTGGCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCtactac                                                                          <  1:89103/149‑1 (MQ=255)
                                                                         tGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCtactac                                                                          <  1:193179/149‑1 (MQ=255)
                                                                           ccACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGCAGGGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCg                                                                        >  1:690777/1‑149 (MQ=255)
                                                                               cGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCtactac                                                                          <  2:1111450/143‑1 (MQ=255)
                                                                                  aaGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATGTGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCgg                                                                       <  2:999724/143‑1 (MQ=255)
                                                                                   agagGGACAGCGGTTGGGTCATCTGGATTTGCTGGTGGATAGGGCCGACGAGCAAAACTTTGTCTATAGAACGTTGCCACAGGCCTATTCGGTGCCGTGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGaaa                                                                <  1:629947/149‑1 (MQ=255)
                                                                                     agGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCATGCAGCGGTAAATCGACCTACTACCGGCTGGAAAcc                                                              <  1:1234873/149‑1 (MQ=255)
                                                                                      gggACAGCAGTTGGGTCTTCTGGATTTGTTGGTGCATATGGGCGTAGAGCAAAACTTTGTCTATCAGATCTCGCCGCAGCAATTTTCGGTGCCGGGCTTTACCTTCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTgg                                                                   <  2:827714/143‑1 (MQ=255)
                                                                                         aCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTcc                                                          >  1:1137898/1‑149 (MQ=255)
                                                                                          catcagTTGTGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTCTGTCTATCAGATTTGGCCGCAGCAATATTCGTTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCt                                                         <  1:371252/146‑1 (MQ=255)
                                                                                          cagcagTTGGGTCATCCGAATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCt                                                         <  1:258923/149‑1 (MQ=255)
                                                                                           agcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTg                                                        <  1:305560/149‑1 (MQ=255)
                                                                                           agcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGGCAATCTTCCTACTACCGGCTGGAAACCTTCCTg                                                        <  1:959357/149‑1 (MQ=255)
                                                                                            gcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCt                                                             >  2:420636/1‑143 (MQ=255)
                                                                                              agTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAGCTTTGTCGATCAGATTTGTCCGCAGCAATGTTCGGTGCCGGCCTTTACATCCGGCGCAGCCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTc                                                           >  2:556063/1‑143 (MQ=255)
                                                                                                ttGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGGGCAAACCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGGGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCGACTACCGGCTGGAAACCTTCCt                                                         <  2:528772/143‑1 (MQ=255)
                                                                                                  gggTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAg                                                 <  1:722404/149‑1 (MQ=255)
                                                                                                  gggTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAg                                                 <  1:464524/149‑1 (MQ=255)
                                                                                                  gggTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTATCCCTGTCCGATGTGGCCGCAGCAATAATCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAg                                                 <  1:38006/149‑1 (MQ=255)
                                                                                                     tCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAg                                                    >  2:981989/1‑143 (MQ=255)
                                                                                                        tCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGcc                                           >  1:325570/1‑149 (MQ=255)
                                                                                                         cTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAAACTTCCTGTTCGAAgg                                                >  2:395853/1‑143 (MQ=255)
                                                                                                              tgtgttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCa                                     <  1:387079/147‑1 (MQ=255)
                                                                                                              tgtggtgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCa                                     <  1:1093383/144‑1 (MQ=255)
                                                                                                                tgttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAAc                                   <  1:1219744/149‑1 (MQ=255)
                                                                                                                 gttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGAAATACTAGAAATCGTTCTTTCCAGTACGTTCACGCAGCGGTAAATCGACCTACTACCGGCTGTATACCTTCCTGTTAGAAGGTAGCCAGGGCAACg                                  >  1:630195/1‑149 (MQ=255)
                                                                                                                  ttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGa                                 >  1:799027/1‑149 (MQ=255)
                                                                                                                     gcgcagctGGGCGAAGAGCAAACCTCTTTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCCTTCCCTACCGCGCACGCAGCCGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCaa                                    <  2:490903/136‑1 (MQ=255)
                                                                                                                     gTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCt                              >  1:987890/1‑149 (MQ=255)
                                                                                                                     gTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGGGCCGGGCATTACCTGCCGCGAACGCAACGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAGGGCGGCCAGGACAACGACCt                              >  1:130995/1‑149 (MQ=255)
                                                                                                                        cATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGa                                 >  2:524282/1‑143 (MQ=255)
                                                                                                                          tatTGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCCACGAcc                               >  2:473274/1‑143 (MQ=255)
                                                                                                                          tatGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGAAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGagtg                         >  1:1119141/1‑146 (MQ=255)
                                                                                                                              ggCGAAGAGCTAAACTTTGCCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGAt                           <  2:758214/143‑1 (MQ=255)
                                                                                                                              ggCGAAGAGCAAAACTTTGCCTATCCGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGAt                           <  2:58539/143‑1 (MQ=255)
                                                                                                                                cGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACACAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTCGAAGGCAGCCAGGGCAACGACCTGATGGCCTACa                   >  1:925610/1‑149 (MQ=255)
                                                                                                                                 gAAGAGCAAAACTTTGTCGATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCACACGAAGCGGTAACTCGACCTACAACCAGCTGGAAACCTTCCTGTTAGAAGGCAGCCGGGGCAACCACCTGGTGGa                        >  2:20136/1‑143 (MQ=255)
                                                                                                                                  aaGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCCGGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGCTGGACGACAGc                 >  1:6450/1‑149 (MQ=255)
                                                                                                                                   agagCAAAACTTTGTCTATCAGAATTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACt                      >  2:771898/1‑143 (MQ=255)
                                                                                                                                     agCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCaaa              <  1:772545/149‑1 (MQ=255)
                                                                                                                                     agCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCaaa              <  1:710770/149‑1 (MQ=255)
                                                                                                                                        aaaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAggag           >  1:316307/1‑149 (MQ=255)
                                                                                                                                         aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCa                >  2:602293/1‑143 (MQ=255)
                                                                                                                                          aaCTTTGTCTATCAGATTTGGCCGCAGCACTATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCCCCGCCCTGATGGACTACAGCACAGGGCa         >  1:345576/1‑149 (MQ=255)
                                                                                                                                           cccttGTCTATGAGACTTGTCCGCAGCAATATTCGGTGCCGGGCTTTACCTTCCGCGCACGCAGCGGTCACCCGACCTACTTCCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCaaa              <  2:125643/140‑1 (MQ=255)
                                                                                                                                           aCTTTGTCTATCAGAGTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCaaa              <  2:611022/143‑1 (MQ=255)
                                                                                                                                            cTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAgg       <  1:711816/149‑1 (MQ=255)
                                                                                                                                             tttGTCTAGCAGATTTGGCCGCAGCAATACTCGGTGCCGGTCTTTACCTACCGCGCACGCAGCGGGAAATCGCCCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGTCCTGATGGACTACAGCAAaga            <  2:341496/143‑1 (MQ=255)
                                                                                                                                                 tCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCAACGGCCACGACCGGATGGACTACAGCAAAGAGCAGGTGATc  <  1:1286310/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATC  >  CP009273/326821‑327114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: