Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,424,363 T→A 100% T138S (ACC→TCC)  tsaC ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,424,3630TA99.1% 721.2 / ‑2.5 217T138S (ACC→TCC) tsaCtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (108/107);  total (109/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 224 total aligned reads displayed.

TTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTG  >  CP009273/3424220‑3424506
                                                                                                                                               |                                                                                                                                               
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                             <  1:1062328/148‑1 (MQ=255)
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                             >  2:3713458/1‑148 (MQ=255)
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                             >  2:1220887/1‑148 (MQ=255)
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                             >  2:650466/1‑148 (MQ=255)
 tcccccGTCCGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACACCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                            <  1:1369450/145‑1 (MQ=255)
 tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                            <  1:1311296/148‑1 (MQ=255)
  cACCCGTCAGGGCATCGCGTATTTCTTAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAAc                                                                                                                                           <  1:1471716/148‑1 (MQ=255)
    cccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGTACTGGAAGAAACCa                                                                                                                                         >  2:3576670/1‑148 (MQ=255)
     ccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                        <  1:3563235/148‑1 (MQ=255)
        tCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCgg                                                                                                                                     <  1:2415815/148‑1 (MQ=255)
         cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCACTCCACTCAAGTTGGCACTGGAAGAAACCAGCGGt                                                                                                                                    >  2:3696484/1‑148 (MQ=255)
         cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGt                                                                                                                                    <  2:3509650/148‑1 (MQ=255)
            ggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTa                                                                                                                                 >  1:559385/1‑148 (MQ=255)
              cATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTAcc                                                                                                                               <  2:1036444/148‑1 (MQ=255)
               aTCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCa                                                                                                                              <  2:1367205/148‑1 (MQ=255)
                tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             >  1:1927317/1‑148 (MQ=255)
                tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             >  1:676955/1‑148 (MQ=255)
                 cgcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATa                                                                                                                            <  2:1047724/148‑1 (MQ=255)
                   cgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAg                                                                                                                          >  1:421729/1‑148 (MQ=255)
                    tgATTTCTGATGGATTTAAACGCCGCGCCGTTTCACCATGCACATCGGGGAACGCCGCGCCTATGTGTGCGCGAACTTCGTCTACTGTTCGACAAGTTGGCTATCCTCTCAAGTTGGCACTGTAAGATACGAGGGGTTTACCATaatc                                                                                                                         <  2:3402829/147‑3 (MQ=255)
                     gATTTCTGAAGGATTTAAACGCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCCGCGGTTTACCATAAGCCt                                                                                                                       >  2:1638840/1‑148 (MQ=255)
                      aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCt                                                                                                                       >  2:3770166/1‑148 (MQ=255)
                          cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGca                                                                                                                   <  1:848122/148‑1 (MQ=255)
                            gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                                                                                 >  2:2008921/1‑148 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                >  1:2053949/1‑148 (MQ=255)
                              aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                                                                               >  1:2254514/1‑148 (MQ=255)
                              aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                                                                               >  1:2145445/1‑148 (MQ=255)
                              aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                                                                               <  1:2765388/148‑1 (MQ=255)
                              aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                                                                               <  1:847691/148‑1 (MQ=255)
                                 aTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCa                                                                                                            <  1:2213577/148‑1 (MQ=255)
                                  tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCaa                                                                                                           >  1:2351068/1‑148 (MQ=255)
                                   ttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAac                                                                                                          >  1:97386/1‑148 (MQ=255)
                                    tAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGGGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCATGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAacc                                                                                                         <  1:3337543/148‑1 (MQ=255)
                                     aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAacca                                                                                                        >  2:2926738/1‑148 (MQ=255)
                                       aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                      >  1:1908111/1‑148 (MQ=255)
                                       aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                      >  2:1341889/1‑148 (MQ=255)
                                            cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACgg                                                                                                 <  2:1321682/148‑1 (MQ=255)
                                            cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACgg                                                                                                 >  1:1768053/1‑148 (MQ=255)
                                                cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                                                                             >  1:2631812/1‑148 (MQ=255)
                                                cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                                                                             <  2:1201374/148‑1 (MQ=255)
                                                  tttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTc                                                                                           >  2:3010511/1‑148 (MQ=255)
                                                   ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                          >  2:494015/1‑148 (MQ=255)
                                                   ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                          >  1:2139734/1‑148 (MQ=255)
                                                    tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCgg                                                                                         <  2:1124698/148‑1 (MQ=255)
                                                      aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                       >  1:3680234/1‑148 (MQ=255)
                                                      aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                       <  1:3690162/148‑1 (MQ=255)
                                                        cAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAc                                                                                     <  2:780660/148‑1 (MQ=255)
                                                        cAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAc                                                                                     <  1:2361225/148‑1 (MQ=255)
                                                         aGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACt                                                                                    >  2:2145181/1‑148 (MQ=255)
                                                           gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACCAGGTGGCAATCCACTAAAGTTGGCAATGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                  >  2:243153/1‑148 (MQ=255)
                                                           gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAANTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                  >  1:1938929/1‑148 (MQ=255)
                                                           gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                  >  2:3120202/1‑148 (MQ=255)
                                                            cacaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                                                                                 <  2:2005468/148‑1 (MQ=255)
                                                                aCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAg                                                                             >  1:3401077/1‑148 (MQ=255)
                                                                    ggAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAg                                                                         <  1:3814933/148‑1 (MQ=255)
                                                                    acgaCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAg                                                                         >  2:778114/4‑148 (MQ=255)
                                                                     gAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACAAACGGATGGTCGGTGACTCGTACAGCAAGc                                                                        >  1:2677035/1‑148 (MQ=255)
                                                                      aaCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCg                                                                       <  2:1540207/148‑1 (MQ=255)
                                                                      aaCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCg                                                                       >  1:3728852/1‑148 (MQ=255)
                                                                        cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGaa                                                                     <  2:2230840/148‑1 (MQ=255)
                                                                        cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGaa                                                                     <  1:3296719/148‑1 (MQ=255)
                                                                              gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAg                                                               <  1:314459/148‑1 (MQ=255)
                                                                                cAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCg                                                             >  2:633073/1‑148 (MQ=255)
                                                                                  aaTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGc                                                           <  2:2515877/148‑1 (MQ=255)
                                                                                           cgAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAAcc                                                  >  1:2372026/1‑148 (MQ=255)
                                                                                               cTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAgcg                                              >  1:165909/1‑148 (MQ=255)
                                                                                                 tCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGcgcg                                            >  1:888523/1‑148 (MQ=255)
                                                                                                 tCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGcgcg                                            <  2:1977561/148‑1 (MQ=255)
                                                                                                 tCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGcgcg                                            >  2:3689881/1‑148 (MQ=255)
                                                                                                 tCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGcgcg                                            <  2:45970/148‑1 (MQ=255)
                                                                                                    tCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtg                                         <  2:1839829/148‑1 (MQ=255)
                                                                                                    tCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtg                                         <  2:421729/148‑1 (MQ=255)
                                                                                                       aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      >  2:2731683/1‑148 (MQ=255)
                                                                                                       aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      >  2:3103505/1‑148 (MQ=255)
                                                                                                       aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      <  1:269669/148‑1 (MQ=255)
                                                                                                       aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      >  1:3183614/1‑148 (MQ=255)
                                                                                                        cTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                     >  1:3209563/1‑148 (MQ=255)
                                                                                                          gTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCg                                   >  2:2267723/1‑148 (MQ=255)
                                                                                                           ttCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGc                                  >  2:714110/1‑148 (MQ=255)
                                                                                                            tCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCa                                 <  2:2999045/148‑1 (MQ=255)
                                                                                                            tCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCa                                 >  2:3380944/1‑148 (MQ=255)
                                                                                                            acGACAAGGTGCCAATCCACTCAAGTGCGCACTGGAAGAAACCCGCGGTTTACCATAAGCCTGGCACAAAGCACCCACCAACGGATGGACGGTGACTCGTACAGCAAGCGAATCAACGCGGCCCGTCAACCAGCGCGGTGTTGTCGCa                                 <  2:3210836/147‑1 (MQ=255)
                                                                                                                   ggTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCgg                          >  1:3344299/1‑148 (MQ=255)
                                                                                                                   ggTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCgg                          >  2:1719588/1‑148 (MQ=255)
                                                                                                                       gCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGaaa                      >  2:1299670/1‑148 (MQ=255)
                                                                                                                        cAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAg                     <  2:3658145/148‑1 (MQ=255)
                                                                                                                         aaTCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAACGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGa                    <  1:1220887/148‑1 (MQ=255)
                                                                                                                            ccACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACaa                 <  2:614603/148‑1 (MQ=255)
                                                                                                                               cTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAgg              >  1:3106938/1‑148 (MQ=255)
                                                                                                                                 cAAGTTGGCACTGGAAGAAACCAGCTGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTg            <  2:3311592/148‑1 (MQ=255)
                                                                                                                                 cAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTg            <  2:3680234/148‑1 (MQ=255)
                                                                                                                                  aaGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGa           <  2:1593371/148‑1 (MQ=255)
                                                                                                                                    gTTTGCACTGGAAGAAACCAGCGTTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa         <  2:2199386/148‑1 (MQ=255)
                                                                                                                                    gTTGGCACTGTAAGAAACCAGGGTTTTACAATACCCCTGGCAAAAAGAAACCACCAAGGGTTGGTCGGTGACTCGTACACAAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa         <  1:1634990/148‑1 (MQ=255)
                                                                                                                                       ggCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGTATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGa      <  2:794425/148‑1 (MQ=255)
                                                                                                                                        gCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAc     <  2:2585033/148‑1 (MQ=255)
                                                                                                                                         cACTGGAAGAAACCAGCGTTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc    <  2:572663/148‑1 (MQ=255)
                                                                                                                                         cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc    <  1:2926745/148‑1 (MQ=255)
                                                                                                                                         cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc    <  2:766101/148‑1 (MQ=255)
                                                                                                                                           cTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTg  >  1:1300057/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTG  >  CP009273/3424220‑3424506

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