Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,885 G→T 11.6% Q36K (CAG→AAG)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,8850GT11.6% 797.5 / 52.4 241Q36K (CAG→AAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (90/123);  new base T (12/16);  total (102/139)
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.56e-01

Only 100 of 243 total aligned reads displayed.

GCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAG  >  CP009273/3958738‑3959029
                                                                                                                                                   |                                                                                                                                                 
  tGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGaa                                                                                                                                                 <  2:3227053/148‑1 (MQ=255)
   gATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAAt                                                                                                                                                <  2:1130643/148‑1 (MQ=255)
   gATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAAt                                                                                                                                                >  1:926657/1‑148 (MQ=255)
   gATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAAt                                                                                                                                                <  2:46214/148‑1 (MQ=255)
     ttCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATgg                                                                                                                                              <  2:2157190/148‑1 (MQ=255)
      tCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGAGCCTGAATggg                                                                                                                                             <  2:916901/148‑1 (MQ=255)
       cACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGc                                                                                                                                            >  1:1568974/1‑148 (MQ=255)
         ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGt                                                                                                                                          <  1:2135108/148‑1 (MQ=255)
         ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGt                                                                                                                                          <  1:2134432/148‑1 (MQ=255)
           ttGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTAc                                                                                                                                        >  1:1048444/1‑148 (MQ=255)
            tGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACa                                                                                                                                       <  1:3264648/148‑1 (MQ=255)
              cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGt                                                                                                                                     <  2:2678089/148‑1 (MQ=255)
                   ggCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGGCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGa                                                                                                                                <  2:604881/148‑1 (MQ=255)
                       aTCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAcc                                                                                                                            <  1:375429/148‑1 (MQ=255)
                       aTCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAcc                                                                                                                            <  1:279924/148‑1 (MQ=255)
                       aTCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAcc                                                                                                                            <  1:303289/148‑1 (MQ=255)
                        tCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAccc                                                                                                                           <  1:243795/148‑1 (MQ=255)
                        tCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAccc                                                                                                                           <  2:1156454/148‑1 (MQ=255)
                         cGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTAACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCt                                                                                                                         >  2:3175657/1‑148 (MQ=255)
                           cAGGATGAGAGAGACGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                                                                       >  2:601715/1‑148 (MQ=255)
                           cAGGATGAGAGAGAAGATAATTAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                                                                       <  2:1825012/148‑1 (MQ=255)
                           cAGGATGAGAGAGAAGATAATTAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                                                                       <  2:2507008/148‑1 (MQ=255)
                           cAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCttt                                                                                                                       >  1:1190808/1‑148 (MQ=255)
                             ggATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCttttt                                                                                                                     <  2:2754974/148‑1 (MQ=255)
                             ggATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttt                                                                                                                     >  1:1973363/1‑148 (MQ=255)
                               aTGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttttt                                                                                                                   <  1:751622/148‑1 (MQ=255)
                               aTGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGAGCCTGAATGGGCGTACAGTTATGAAACCCttttttt                                                                                                                   >  1:2748478/1‑148 (MQ=255)
                                   gagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCaa                                                                                                               <  2:151264/148‑1 (MQ=255)
                                       gaagaTAATGAAACGTGGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGc                                                                                                           <  2:1925144/148‑1 (MQ=255)
                                        aagaTAATGAAACGTTGACGTAAGAAACGACATCGTTTTCCCGGTAGCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCt                                                                                                          <  1:3235729/148‑1 (MQ=255)
                                           aTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTc                                                                                                        >  2:49843/1‑148 (MQ=255)
                                             aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                      >  2:1608304/1‑148 (MQ=255)
                                             aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                      >  2:1066804/1‑148 (MQ=255)
                                             aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                                                                      >  1:826396/1‑148 (MQ=255)
                                             aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTAc                                                                                                     >  2:1466344/1‑148 (MQ=255)
                                             aaTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTAc                                                                                                     >  1:368743/1‑148 (MQ=255)
                                              aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTAc                                                                                                     >  2:2919742/1‑148 (MQ=255)
                                              aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACa                                                                                                    >  2:1367316/1‑148 (MQ=255)
                                               tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                                                                    <  2:796028/148‑1 (MQ=255)
                                               tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                                                                    <  1:1470121/148‑1 (MQ=255)
                                               tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACaa                                                                                                   >  1:1111961/1‑148 (MQ=255)
                                                gAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAAc                                                                                                  <  1:655869/148‑1 (MQ=255)
                                                 aaaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                                                                                  <  2:1777700/148‑1 (MQ=255)
                                                 aaaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTGTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                                                                                  <  2:286485/148‑1 (MQ=255)
                                                    cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCtt                                                                                               <  1:2221247/148‑1 (MQ=255)
                                                    cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTc                                                                                              <  2:331/148‑1 (MQ=255)
                                                        gACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGAt                                                                                          <  1:1221879/148‑1 (MQ=255)
                                                         aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                         >  2:2588898/1‑148 (MQ=255)
                                                         aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                         <  2:71872/148‑1 (MQ=255)
                                                         aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTGTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                         <  2:1746747/148‑1 (MQ=255)
                                                          cGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCNGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                         >  1:1678146/1‑148 (MQ=255)
                                                           gTAAGAAACGCCCTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                                                                       <  2:2177993/148‑1 (MQ=255)
                                                            tAAGAAACGCAATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAg                                                                                      >  2:1979281/1‑148 (MQ=255)
                                                                      cATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGAGGCAGGGCGAAGTCg                                                                            >  2:3268443/1‑148 (MQ=255)
                                                                       aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                                                                           >  2:1552370/1‑148 (MQ=255)
                                                                        tCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGa                                                                          <  1:872295/148‑1 (MQ=255)
                                                                          gTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaa                                                                        <  1:2951148/148‑1 (MQ=255)
                                                                           ttttCCCGGGACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                       <  2:1679769/148‑1 (MQ=255)
                                                                               cccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                                                                   <  1:1703396/148‑1 (MQ=255)
                                                                               cccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                                                                   <  1:597125/148‑1 (MQ=255)
                                                                                ccGGTACCGGTTTGCGCCTGCCCGGCTACGCCACGACCCGCCAGCGTCAGCGGAAGGGCCAGGGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct                                                                  <  1:2873550/148‑1 (MQ=255)
                                                                                ccGGTACCGGTTGGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct                                                                  <  1:2873970/148‑1 (MQ=255)
                                                                                 cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg                                                                 >  1:2656235/1‑148 (MQ=255)
                                                                                    tACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc                                                              <  1:1116814/148‑1 (MQ=255)
                                                                                    tACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGAGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc                                                              >  1:1026701/1‑148 (MQ=255)
                                                                                      ccGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct                                                             <  1:3178264/148‑1 (MQ=255)
                                                                                       cGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                            <  1:2343702/148‑1 (MQ=255)
                                                                                         gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt                                                          >  2:2041872/1‑148 (MQ=255)
                                                                                          tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa                                                        >  1:3126368/1‑148 (MQ=255)
                                                                                          tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa                                                        >  1:77872/1‑148 (MQ=255)
                                                                                           ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa                                                       <  2:224908/148‑1 (MQ=255)
                                                                                              cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtg                                                     <  2:2313694/148‑1 (MQ=255)
                                                                                              cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgt                                                    <  2:1869298/148‑1 (MQ=255)
                                                                                                ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                  <  1:2004957/148‑1 (MQ=255)
                                                                                                  ggCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                 <  1:227874/147‑1 (MQ=255)
                                                                                                    cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTg                                              <  1:2482460/148‑1 (MQ=255)
                                                                                                           aCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                       <  2:1048376/148‑1 (MQ=255)
                                                                                                            cGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgt                                      >  1:814726/1‑148 (MQ=255)
                                                                                                             gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                     <  2:964650/148‑1 (MQ=255)
                                                                                                             gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                     <  2:3052777/148‑1 (MQ=255)
                                                                                                             gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                     <  2:2618625/148‑1 (MQ=255)
                                                                                                               cACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGtgtg                                   >  1:2034016/1‑148 (MQ=255)
                                                                                                                aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg                                  >  2:3261124/1‑148 (MQ=255)
                                                                                                                    cccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGa                              >  2:2958948/1‑148 (MQ=255)
                                                                                                                    cccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGa                              >  2:1438644/1‑148 (MQ=255)
                                                                                                                      cGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAat                            <  1:1426691/148‑1 (MQ=255)
                                                                                                                      cGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAat                            >  2:3235654/1‑148 (MQ=255)
                                                                                                                        ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAata                           <  2:826396/148‑1 (MQ=255)
                                                                                                                            cGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAg                       <  1:3087971/148‑1 (MQ=255)
                                                                                                                            cGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGc                      >  2:239894/1‑148 (MQ=255)
                                                                                                                              tCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCt                     <  1:1066804/148‑1 (MQ=255)
                                                                                                                              tCAGCGGAAGGGCCAGAGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCtt                    <  2:1096720/148‑1 (MQ=255)
                                                                                                                              tCAGCGGAAGGGCCAGAGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCtt                    <  2:2748478/148‑1 (MQ=255)
                                                                                                                               cAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTa                   >  2:647141/1‑148 (MQ=255)
                                                                                                                                    gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTAt               <  1:1440264/148‑1 (MQ=255)
                                                                                                                                     aaGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTg             <  2:126041/148‑1 (MQ=255)
                                                                                                                                              tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaa    >  1:561099/1‑148 (MQ=255)
                                                                                                                                                ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAg  >  1:136720/1‑148 (MQ=255)
                                                                                                                                                ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAg  >  1:143538/1‑148 (MQ=255)
                                                                                                                                                ccTGAATGGGCGTACAGTTATGAAACACTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAg  <  1:1610162/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAG  >  CP009273/3958738‑3959029

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