Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,909:1 +T 100% coding (82/1266 nt) rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,9091.T95.4% 602.5 / 29.4 194coding (82/1266 nt)rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base . (7/2);  new base T (83/102);  total (90/104)
Fisher's exact test for biased strand distribution p-value = 8.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 230 total aligned reads displayed.

TCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAA  >  CP009273/3958762‑3959052
                                                                                                                                                    |                                                                                                                                               
  tctggATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtt                                                                                                                                               <  2:2872260/145‑1 (MQ=255)
  gCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtt                                                                                                                                               <  2:767337/148‑1 (MQ=255)
  gCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtt                                                                                                                                               <  2:3032455/148‑1 (MQ=255)
    aGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttt                                                                                                                                             <  2:208774/148‑1 (MQ=255)
      gATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttt                                                                                                                                           >  1:89869/1‑148 (MQ=255)
        tGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGACACCCtttttttt                                                                                                                                         >  2:644458/1‑148 (MQ=255)
         gagagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttttttt                                                                                                                                        <  1:3095097/148‑2 (MQ=255)
            agagaAAATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttttcaa                                                                                                                                     >  1:1773446/1‑144 (MQ=255)
             gagtAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttttcaag                                                                                                                                    <  2:1386265/148‑6 (MQ=255)
             gagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttttcaag                                                                                                                                    >  1:722107/1‑143 (MQ=255)
                 agaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCt                                                                                                                                 <  2:3152382/148‑1 (MQ=255)
                  gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCtt                                                                                                                                >  2:3496346/1‑148 (MQ=255)
                  gaTAATGAAACGTTGACGTAAGAAACGCAATCGTTTTCCCGGTACCTGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGTACAGTGCCTGAATGGGAGTACAGTTATGAAACCCTTTTTTTTTCAAGGTCtt                                                                                                                                <  2:2903855/148‑1 (MQ=255)
                     aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                                             >  1:1872513/1‑148 (MQ=255)
                     aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                                             >  1:3175993/1‑148 (MQ=255)
                     aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                                             >  1:3339755/1‑148 (MQ=255)
                     aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                                             <  1:1538683/148‑1 (MQ=255)
                      aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTAc                                                                                                                            >  1:639264/1‑148 (MQ=255)
                       tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                                                                                           >  2:147559/1‑148 (MQ=255)
                       tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                                                                                           <  1:2673884/148‑1 (MQ=255)
                        gAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTGCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACaa                                                                                                                          >  2:3037489/1‑148 (MQ=255)
                        gAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACAGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACatc                                                                                                                         >  2:551794/1‑146 (MQ=255)
                         aaaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTGTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                                                                                                         <  2:517220/148‑1 (MQ=255)
                          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTCCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAcc                                                                                                                        <  2:3153533/148‑1 (MQ=255)
                          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAcc                                                                                                                        <  1:719661/148‑1 (MQ=255)
                          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAcc                                                                                                                        <  2:1406923/148‑1 (MQ=255)
                           aCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCt                                                                                                                       >  1:411472/1‑148 (MQ=255)
                            cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCtt                                                                                                                      >  1:2464712/1‑148 (MQ=255)
                               tGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCgg                                                                                                                   >  2:3324144/1‑148 (MQ=255)
                               tGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCgg                                                                                                                   >  2:2355724/1‑148 (MQ=255)
                                gACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGa                                                                                                                  <  1:601721/148‑1 (MQ=255)
                                gACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGa                                                                                                                  <  1:1812469/148‑1 (MQ=255)
                                   gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                                                                                               <  1:2059722/148‑1 (MQ=255)
                                   gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                                                                                               <  1:2804526/148‑1 (MQ=255)
                                      aGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                                                                                            >  2:676362/1‑148 (MQ=255)
                                      aGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                                                                                            >  1:1868569/1‑148 (MQ=255)
                                       gAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc                                                                                                          >  1:1974727/1‑148 (MQ=255)
                                          aCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa                                                                                                        <  2:173785/148‑1 (MQ=255)
                                           cGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                                                                                       <  1:1027331/148‑1 (MQ=255)
                                           cGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                                                                                       >  1:3303238/1‑148 (MQ=255)
                                             ccATCGTTTGCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                                                                                                     <  1:2159044/148‑1 (MQ=255)
                                                tCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                                                                                                  >  1:3060548/1‑148 (MQ=255)
                                                    tttCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                                              <  2:834829/148‑1 (MQ=255)
                                                    tttCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                                              <  1:2762833/148‑1 (MQ=255)
                                                    tttCCCGCGACCGGTTTGCGCCTGCCCGGCTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                                              <  2:2151100/148‑1 (MQ=255)
                                                      tCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACt                                                                                            <  2:1152782/148‑1 (MQ=255)
                                                        ccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                                                                                          >  1:343510/1‑148 (MQ=255)
                                                         cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct                                                                                         <  1:1363896/148‑1 (MQ=255)
                                                         cGGCACCGGTGTGCGCCGGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGGGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct                                                                                         <  1:2136552/148‑1 (MQ=255)
                                                          ggTACCGTTTTGCGCCTGCCCGGATACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTATTGGGCGTTCAGTTATGCTACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGTATGCAGGGCGAATTCGGAAAACttctg                                                                                        <  1:1715218/148‑1 (MQ=255)
                                                           gTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                                       <  1:3421291/148‑1 (MQ=255)
                                                            tACCGGTTTGCGGCTGCCCATCTACGTCACGACCTGCCTTCGTCAGCGGAAGTGCCATTGCCTGAATGTGCGTACTGTTATGAATCCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                                      <  2:1810603/148‑1 (MQ=255)
                                                            tACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                                      <  2:1872513/148‑1 (MQ=255)
                                                             aCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc                                                                                     >  1:768928/1‑148 (MQ=255)
                                                             aCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc                                                                                     >  1:1511011/1‑148 (MQ=255)
                                                               cGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                                                   <  2:507325/148‑1 (MQ=255)
                                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa                                                                               >  1:3470756/1‑148 (MQ=255)
                                                                   ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa                                                                               <  2:2597407/148‑1 (MQ=255)
                                                                    tGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa                                                                              <  1:362512/148‑1 (MQ=255)
                                                                        ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgtg                                                                          >  2:2606113/1‑148 (MQ=255)
                                                                          tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTCAATGTGtt                                                                        >  2:2493730/1‑148 (MQ=255)
                                                                          tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                                        >  1:849947/1‑148 (MQ=255)
                                                                          tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                                        <  2:2649588/148‑1 (MQ=255)
                                                                          tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                                        <  1:1702355/148‑1 (MQ=255)
                                                                                  tACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAt                                                                <  1:1270692/148‑1 (MQ=255)
                                                                                  tACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAt                                                                <  1:1928454/148‑1 (MQ=255)
                                                                                   aCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTAAAGGGCTTAAAAAACCTTCCGTTTGCGGGCGAAGTAGGAAAAAATGAGTTGGTTTAAATTTGtttttgtgatt                                                               >  2:3462623/1‑141 (MQ=255)
                                                                                    cGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                                              <  1:174405/148‑1 (MQ=255)
                                                                                      tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                                            <  2:3298220/148‑1 (MQ=255)
                                                                                       cACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGtgtg                                                          >  1:2872865/1‑148 (MQ=255)
                                                                                        aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGtgtg                                                          <  2:611379/148‑1 (MQ=255)
                                                                                        aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGtgtg                                                          >  1:1060245/1‑148 (MQ=255)
                                                                                        aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGtgtg                                                          <  2:1956621/148‑1 (MQ=255)
                                                                                          gACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                                        >  1:2299521/1‑148 (MQ=255)
                                                                                                  atcGTCAGCGGAAGGGCCTGTGCATTAATGGGCTTACTGTTATTAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATc                                                <  2:1918067/146‑1 (MQ=255)
                                                                                                    cGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAg                                              >  1:435162/1‑148 (MQ=255)
                                                                                                         gCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTAc                                         >  1:3237868/1‑148 (MQ=255)
                                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                                        >  2:1656805/1‑148 (MQ=255)
                                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                                        >  2:3101760/1‑148 (MQ=255)
                                                                                                           ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTa                                       <  1:48017/148‑1 (MQ=255)
                                                                                                              aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTg                                    >  2:2248537/1‑148 (MQ=255)
                                                                                                                 gCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCtt                                 >  1:128843/1‑148 (MQ=255)
                                                                                                                  ccAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCttt                                <  1:2647696/148‑1 (MQ=255)
                                                                                                                  ccAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCttt                                <  1:1550680/148‑1 (MQ=255)
                                                                                                                   cAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTa                               <  2:1109518/148‑1 (MQ=255)
                                                                                                                    aGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTAc                              >  1:2399196/1‑148 (MQ=255)
                                                                                                                        ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaaa                          >  1:3008010/1‑148 (MQ=255)
                                                                                                                        ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaaa                          <  1:283569/148‑1 (MQ=255)
                                                                                                                        ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaaa                          <  1:1202416/148‑1 (MQ=255)
                                                                                                                        ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaaa                          <  1:936010/148‑1 (MQ=255)
                                                                                                                         cTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAg                         >  2:3058214/1‑148 (MQ=255)
                                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCg                       >  1:1683413/1‑148 (MQ=255)
                                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCg                       >  1:1651188/1‑148 (MQ=255)
                                                                                                                                  cGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCgg                <  1:1401904/148‑1 (MQ=255)
                                                                                                                                          tATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATaa        >  2:1295845/1‑148 (MQ=255)
                                                                                                                                            tGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAg      >  1:780437/1‑148 (MQ=255)
                                                                                                                                             gAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTCAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGACAGCGTAGCCGGTGAAATAAAGt     >  2:904381/1‑148 (MQ=255)
                                                                                                                                             gAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGt     >  2:1964296/1‑148 (MQ=255)
                                                                                                                                               taCCTTTTTTTTTAAAGTGCTTCTAGAAGCTTCGTATTCAGGTCGAAGTAGGAATTCTTGTGTTGTGTTAAATGTGTTTTGCTCATTTTGTGGTAGTATATCAGCTTACTATTGCTTTACGATAGCGTATCCGGTGAAATAAAGTCaa  <  2:3175993/147‑1 (MQ=255)
                                                                                                                                               aaCCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCa   <  1:1673250/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAA  >  CP009273/3958762‑3959052

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