Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,750 T→A 42.9% K81* (AAG→TAG)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,7500TA42.9% 99.1 / 248.3 182K81* (AAG→TAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (71/33);  new base A (41/37);  total (112/70)
Fisher's exact test for biased strand distribution p-value = 4.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 184 total aligned reads displayed.

TCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACA  >  CP009273/3958606‑3958897
                                                                                                                                                |                                                                                                                                                   
      gTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGc                                                                                                                                            >  2:1221787/1‑148 (MQ=255)
             gTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCa                                                                                                                                     <  2:2046351/148‑1 (MQ=255)
             gTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCa                                                                                                                                     <  2:1119142/148‑1 (MQ=255)
             gTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCa                                                                                                                                     <  2:2892377/148‑1 (MQ=255)
              tAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                    >  1:3414046/1‑148 (MQ=255)
                gCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGa                                                                                                                                  >  1:969959/1‑148 (MQ=255)
                   aTCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATga                                                                                                                               <  1:2873044/148‑1 (MQ=255)
                    tCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgag                                                                                                                              >  2:1973859/1‑148 (MQ=255)
                    tCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgag                                                                                                                              >  2:3112200/1‑148 (MQ=255)
                    tCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATgag                                                                                                                              >  1:698729/1‑148 (MQ=255)
                    tCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATgag                                                                                                                              >  1:697842/1‑148 (MQ=255)
                     caccacCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATgaga                                                                                                                             <  1:1058658/148‑1 (MQ=255)
                      accaccGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATgagag                                                                                                                            <  1:3416187/148‑1 (MQ=255)
                      accaccGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATgagag                                                                                                                            >  2:264343/1‑148 (MQ=255)
                        caccGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagagag                                                                                                                          <  2:3557879/148‑1 (MQ=255)
                        caccGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagagag                                                                                                                          >  2:3240423/1‑148 (MQ=255)
                             tAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGAt                                                                                                                     <  2:458752/148‑1 (MQ=255)
                               aGCCAGACCCAGCTTCAGGCCAGTCGCTTCCGCCAGCGGTTCTTCGTCGGCCTGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATaa                                                                                                                   <  2:3073360/148‑1 (MQ=255)
                                gCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGNGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAAt                                                                                                                  <  1:1802783/148‑1 (MQ=255)
                                gCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAAt                                                                                                                  >  2:47180/1‑148 (MQ=255)
                                  cAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGa                                                                                                                <  1:574877/148‑1 (MQ=255)
                                        cAGCTTCAGGCCAGTAGCTTCCGCCAGCTGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGtt                                                                                                          <  2:654753/148‑1 (MQ=255)
                                        cAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGtt                                                                                                          >  2:445571/1‑148 (MQ=255)
                                         aGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGCGAAGATAATGAAACGTTg                                                                                                         >  1:434469/1‑148 (MQ=255)
                                         aGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTg                                                                                                         >  2:3172528/1‑148 (MQ=255)
                                          gCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGa                                                                                                        <  1:165923/148‑1 (MQ=255)
                                            ttCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACg                                                                                                      <  2:474927/148‑1 (MQ=255)
                                                  ccAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGaa                                                                                                >  2:2441206/1‑148 (MQ=255)
                                                     gTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACg                                                                                             >  2:1793265/1‑148 (MQ=255)
                                                      tAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGc                                                                                            >  2:991738/1‑148 (MQ=255)
                                                       aGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCACTCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGcc                                                                                           >  2:3162841/1‑148 (MQ=255)
                                                       aGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGcc                                                                                           >  2:1017284/1‑148 (MQ=255)
                                                       aGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGcc                                                                                           <  2:2209945/148‑1 (MQ=255)
                                                        gCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCa                                                                                          <  2:2695447/148‑1 (MQ=255)
                                                        gCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCa                                                                                          <  2:2693636/148‑1 (MQ=255)
                                                         cTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCAt                                                                                         >  2:3311483/1‑148 (MQ=255)
                                                           tCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCg                                                                                       >  1:2918551/1‑148 (MQ=255)
                                                             cGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCCGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGtt                                                                                     <  2:992622/148‑1 (MQ=255)
                                                              gCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGttt                                                                                    <  1:3251724/148‑1 (MQ=255)
                                                              gCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGttt                                                                                    <  1:2957481/148‑1 (MQ=255)
                                                              gCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGttt                                                                                    >  2:1303489/1‑148 (MQ=255)
                                                              gCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGttt                                                                                    >  2:1303369/1‑148 (MQ=255)
                                                               ccAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGtttt                                                                                   <  2:2840493/148‑1 (MQ=255)
                                                                  gCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTccc                                                                                >  2:905533/1‑148 (MQ=255)
                                                                   cGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCg                                                                               >  2:515470/1‑148 (MQ=255)
                                                                   cGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCg                                                                               >  1:2289832/1‑148 (MQ=255)
                                                                   cGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCg                                                                               >  1:2289878/1‑148 (MQ=255)
                                                                      ttCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTa                                                                            >  2:2951729/1‑148 (MQ=255)
                                                                      ttCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTa                                                                            >  2:1935563/1‑148 (MQ=255)
                                                                      ttCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTa                                                                            >  2:1935138/1‑148 (MQ=255)
                                                                      ttCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTa                                                                            >  1:2947668/1‑148 (MQ=255)
                                                                      ttCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTa                                                                            >  1:3244260/1‑148 (MQ=255)
                                                                       tCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTAc                                                                           <  1:1119796/148‑1 (MQ=255)
                                                                        cTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTAcc                                                                          >  1:1758079/1‑148 (MQ=255)
                                                                            gTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGtt                                                                      <  2:881697/148‑1 (MQ=255)
                                                                            gTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGtt                                                                      <  1:2006594/148‑1 (MQ=255)
                                                                            gTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGtt                                                                      <  2:1557885/148‑1 (MQ=255)
                                                                              cGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTg                                                                    <  1:3098854/148‑1 (MQ=255)
                                                                                 cATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCTCGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTgcgc                                                                 >  2:2361515/1‑148 (MQ=255)
                                                                                  aTGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGcc                                                                <  1:964650/148‑1 (MQ=255)
                                                                                     gATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                             >  2:2012813/1‑148 (MQ=255)
                                                                                      aTCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGcc                                                            <  2:838036/148‑1 (MQ=255)
                                                                                          gCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGc                                                        >  2:632990/1‑148 (MQ=255)
                                                                                             cGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTAc                                                     >  1:784005/1‑148 (MQ=255)
                                                                                             cGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTAc                                                     >  1:2452998/1‑148 (MQ=255)
                                                                                                cAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTc                                                  >  2:3189589/1‑148 (MQ=255)
                                                                                                cAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACgcc                                                  >  1:660069/1‑146 (MQ=255)
                                                                                                cAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTc                                                  >  1:1570577/1‑148 (MQ=255)
                                                                                                   gTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACg                                               >  1:137478/1‑148 (MQ=255)
                                                                                                   gTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACg                                               <  1:615085/148‑1 (MQ=255)
                                                                                                   gTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACg                                               <  1:2398015/148‑1 (MQ=255)
                                                                                                        cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGTGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                          >  1:781794/1‑148 (MQ=255)
                                                                                                        cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAATATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                          >  1:2564475/1‑148 (MQ=255)
                                                                                                        cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                          >  2:3247139/1‑148 (MQ=255)
                                                                                                            tCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAg                                      >  2:2623813/1‑148 (MQ=255)
                                                                                                             cGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGc                                     >  2:2869558/1‑148 (MQ=255)
                                                                                                              ggTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCACCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCg                                    <  2:2494696/148‑1 (MQ=255)
                                                                                                               gTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGt                                   >  1:2181171/1‑148 (MQ=255)
                                                                                                               gTGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGt                                   <  1:89674/148‑1 (MQ=255)
                                                                                                                tGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTc                                  >  1:2976556/1‑148 (MQ=255)
                                                                                                                tGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTc                                  <  2:1688197/148‑1 (MQ=255)
                                                                                                                tGCCATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTc                                  <  2:786210/148‑1 (MQ=255)
                                                                                                                   cATAATTAAGGCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGc                               >  2:1323457/1‑148 (MQ=255)
                                                                                                                    aTAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCg                              >  2:3110044/1‑148 (MQ=255)
                                                                                                                    aTAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCg                              >  1:2862932/1‑148 (MQ=255)
                                                                                                                     tAATTAAGGCACGCGGCTGATTCACCGAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCgg                             <  2:2160236/148‑1 (MQ=255)
                                                                                                                           aGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTGTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGc                       <  1:403214/148‑1 (MQ=255)
                                                                                                                            ggCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGcc                      >  2:656853/1‑148 (MQ=255)
                                                                                                                            ggCACGCGGCTGATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGcc                      >  2:2788251/1‑148 (MQ=255)
                                                                                                                                   ggCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCt               <  2:1939593/148‑1 (MQ=255)
                                                                                                                                     cTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGa             >  1:2307176/1‑148 (MQ=255)
                                                                                                                                     cTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGa             <  1:1702503/148‑1 (MQ=255)
                                                                                                                                       gATTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAAt           >  2:2143780/1‑148 (MQ=255)
                                                                                                                                        aTTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATg          >  2:2606635/1‑148 (MQ=255)
                                                                                                                                        aTTCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATg          <  2:1809026/148‑1 (MQ=255)
                                                                                                                                         ttCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATgg         >  2:2375885/1‑148 (MQ=255)
                                                                                                                                          tCACCTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATggg        >  2:226659/1‑148 (MQ=255)
                                                                                                                                             ccTAGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGt     >  2:1385748/1‑148 (MQ=255)
                                                                                                                                                tGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACa  <  2:2452998/148‑1 (MQ=255)
                                                                                                                                                tGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACa  <  2:784005/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACA  >  CP009273/3958606‑3958897

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: