Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,433,555 G→T I278I (ATC→ATA)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,433,5550GT98.5% 192.5 / NA 67I278I (ATC→ATA) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (31/34);  total (32/35)
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.98e-01

ATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAA  >  CP009273/3433412‑3433699
                                                                                                                                               |                                                                                                                                                
aTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCCCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCtatata                                                                                                                                             <  1:169611/149‑1 (MQ=255)
      gCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCAAGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAAGCTCGGTACGCTGTACCAGATCGCCTATATAGTGGAt                                                                                                                                       >  1:674504/1‑149 (MQ=255)
        cGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCCGATCACCTATATAGTGGATAg                                                                                                                                     >  1:662907/1‑149 (MQ=255)
          gTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCCCAGTAAGAGACTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATGTAGGGGATAGCt                                                                                                                                   >  2:935180/1‑149 (MQ=255)
              ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAGTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTg                                                                                                                               >  2:497096/1‑149 (MQ=255)
               cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACTATGTTTTGGCGTCTTAAGGAGCTCAACCTCGGTACGCTGTACCAGAGGACCTATATAGTGGATAGCTTCTGc                                                                                                                              >  2:932389/1‑149 (MQ=255)
                 gtcttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCtt                                                                                                                            <  1:321848/146‑1 (MQ=255)
                  ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCCGATAACCAATATAGTGGATAACTTCTggctt                                                                                                                           >  1:359699/1‑145 (MQ=255)
                   tttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTa                                                                                                                          <  1:955981/149‑1 (MQ=255)
                     tCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTCAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGGGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAg                                                                                                                        >  2:33283/1‑149 (MQ=255)
                      ccAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAgg                                                                                                                       <  1:925546/149‑1 (MQ=255)
                        aGGCGCAGGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCCCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                                     <  2:464679/149‑1 (MQ=255)
                                     gagaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg                                                                                                        <  2:817085/149‑1 (MQ=255)
                                      agaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCggg                                                                                                       >  2:570114/1‑148 (MQ=255)
                                        aCAGTTCACGGGGAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                                     <  2:179696/149‑1 (MQ=255)
                                        aCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                                     <  2:785382/149‑1 (MQ=255)
                                             ccACGGGAAGCCAGCACGGCTTTAATCCCAGGAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                                                <  2:229418/149‑1 (MQ=255)
                                                 gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                                            <  2:482932/149‑1 (MQ=255)
                                                    acggcagtctggcTTTAATCTTAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGGACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGa                                                                                         <  1:460057/137‑1 (MQ=255)
                                                    aaGCCCGCACTTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCata                                                                                         <  2:878056/149‑3 (MQ=255)
                                                      gCCAGCACGTCTTTAATCTCAGTCAAAGATTTTTTACCAAGGGTAGGCGTGTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAAAGCAGACAGTCAATTCCAGATc                                                                                       >  2:483977/1‑149 (MQ=255)
                                                      gCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                                       <  2:318622/149‑1 (MQ=255)
                                                         aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGTATCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTc                                                                                    >  1:747339/1‑149 (MQ=255)
                                                            aCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCacc                                                                                 >  2:293411/1‑147 (MQ=255)
                                                              gTCTTTGATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGTAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCCGATCGTCAACAg                                                                               >  2:252547/1‑149 (MQ=255)
                                                              gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTCCCAGATCACCTATATAGTGGATAGCTTCTGCTATAAGGCAGTTAACAGAGCAGACAGTCAATTCCAGATCGTCAACAg                                                                               >  2:895458/1‑149 (MQ=255)
                                                              gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAg                                                                               <  2:616833/149‑1 (MQ=255)
                                                              gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTGGAGGCAGTTAGCAGAGCGGACAGTCAATTCCCAATCGTCAACAg                                                                               >  1:414899/1‑149 (MQ=255)
                                                                 tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                            <  1:160794/149‑1 (MQ=255)
                                                                 tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                            >  2:731254/1‑149 (MQ=255)
                                                                    aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGcaca                                                                         <  1:919387/149‑1 (MQ=255)
                                                                    aaTCCCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTTGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgca                                                                         <  2:896017/149‑1 (MQ=255)
                                                                        tcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGcagcag                                                                     >  1:752229/1‑149 (MQ=255)
                                                                          aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTGAGGCAGTAAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAgga                                                                   >  1:707116/1‑149 (MQ=255)
                                                                          aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAgga                                                                   >  1:874112/1‑149 (MQ=255)
                                                                          aGTAAGAGAGTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAgga                                                                   <  1:416805/149‑1 (MQ=255)
                                                                                  aTTTTTTACGAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGgatcggatcg                                                           >  2:886502/1‑149 (MQ=255)
                                                                                          ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTgg                                                   <  2:312906/149‑1 (MQ=255)
                                                                                              ggTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtc                                               <  1:272616/149‑1 (MQ=255)
                                                                                               gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAAGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtct                                              <  1:39390/149‑1 (MQ=255)
                                                                                                ttcggctgtttAAGGAGCTCAACCTCGGTACGCTGCACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCTAGATCGTCAACAGGGCGCAGCTGGATCGGATCGAACTCTGGTTtctc                                             <  1:533341/141‑1 (MQ=255)
                                                                                                   ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGCTCGGATCGAACTCTGGTTTCtcttc                                          <  1:479761/149‑1 (MQ=255)
                                                                                                    gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCCCCTATCTAGCGGCTAGCTTCTGCTCTAAGGCAGTTAGCAGAGCGGACAGTCAATGCCAGATCGTCAACAGGGCGCAGCAGGCTCGCACCGAACACTGGGTCCtctcct                                         >  2:780876/1‑149 (MQ=255)
                                                                                                    gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                         <  2:774549/149‑1 (MQ=255)
                                                                                                       ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                                      <  2:28894/149‑1 (MQ=255)
                                                                                                        tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                                     <  2:430663/149‑1 (MQ=255)
                                                                                                                gCTCAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGc                             <  1:823289/149‑1 (MQ=255)
                                                                                                                 cgcAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATGCCAGATCGTCAACAGGGCGCGGCAGAGTCGGGTCGGCCGCTGAGTTCGCTTCGTTCACGGCAGGCt                            >  2:566087/3‑149 (MQ=255)
                                                                                                                   cAACCTCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGa                          <  1:832778/149‑1 (MQ=255)
                                                                                                                    aaCCTCGGTACGCTGTACCAGATCACCTAGATAGTGGATAGCTTCTGCTTTAAGGCAGTTATCAGAGAGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGATTCTCTTCTTTCACTTCAGGCTGAc                         >  1:584571/1‑149 (MQ=255)
                                                                                                                       cTCGGTACGCTGTACCAGAGCACCTAGATAGTGGATAGCTTGTGCTTTAAGGGAGTTAGCAGAGCGGACAGTCAATTCCAGATAGTCGACAGGGCGCAGCAGGATCGGATCGGACGTCGGAGTGTCTTCTTTCACTTCAAGATGACgga                      >  2:181476/1‑147 (MQ=255)
                                                                                                                        tCGGTACGCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc                     >  1:739997/1‑149 (MQ=255)
                                                                                                                             aggCTGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCAc                >  2:583868/3‑149 (MQ=255)
                                                                                                                               gcagcaacAGATCACCTATATGGTGGATAGCTTCTGCTTTCAGGCAGTTAGCAGAGCGGGCAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTGCGCTTCTTTCACTTCAGGCTGACGTACATCACGt              >  2:877998/8‑149 (MQ=255)
                                                                                                                                 tGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTaa            <  1:944948/149‑1 (MQ=255)
                                                                                                                                 tGTACCAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTGTCAGTTCAGGCTGACGTACATCACGTaa            <  1:715165/149‑1 (MQ=255)
                                                                                                                                 tGTACCAGATCACCTATATAGTGGATAGATTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTaa            >  1:526790/1‑149 (MQ=255)
                                                                                                                                  gTACCAGATCACCTATATAGTCGCTAGCTTCTGCTTTAAGGCAGTTAGCAGCGCGGACAGTCAATTTCAGATCGTCAACAGCGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAg           <  2:495611/149‑1 (MQ=255)
                                                                                                                                     ccAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACGTCACGTAAGTCa        >  1:86324/1‑149 (MQ=255)
                                                                                                                                      cAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCaa       <  2:888228/149‑1 (MQ=255)
                                                                                                                                      cAGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTCAGTATCTTTTACTGTCACTTCACGACCGGCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCaa       >  1:943125/1‑149 (MQ=255)
                                                                                                                                       aGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCCCGTAAGCCAAc      >  2:169526/1‑149 (MQ=255)
                                                                                                                                       aGATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAAc      <  1:932389/149‑1 (MQ=255)
                                                                                                                                        gATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACg     >  2:772749/1‑149 (MQ=255)
                                                                                                                                        gATCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACg     >  2:806752/1‑149 (MQ=255)
                                                                                                                                         aTCACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGa    <  2:780698/149‑1 (MQ=255)
                                                                                                                                           cACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaaa  <  1:731254/149‑1 (MQ=255)
                                                                                                                                           cACCTATATAGTGGATAGCTTCCGCTTGAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaaa  <  2:377190/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAA  >  CP009273/3433412‑3433699

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: ‑