Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 829,223 G→A 7.8% S233N (AGC→AAC)  dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273829,2230GA 7.8% 1228.1 / 64.5 345S233N (AGC→AAC) dinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (158/159);  new base A (15/12);  total (173/172)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 362 total aligned reads displayed.

CTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTG  >  CP009273/829074‑829371
                                                                                                                                                     |                                                                                                                                                    
 tGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCGAACCATGCGCTGGTGATGGCGGCGATGGAaagc                                                                                                                                                     >  1:618323/1‑150 (MQ=255)
      tcgtTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAaagcgaagc                                                                                                                                                <  2:2217852/150‑1 (MQ=255)
        gtTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCgg                                                                                                                                              <  1:1017225/150‑1 (MQ=255)
         tttAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGt                                                                                                                                             <  1:2769531/150‑1 (MQ=255)
           tAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTGTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTAt                                                                                                                                           >  1:268092/1‑150 (MQ=255)
             aGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTg                                                                                                                                         >  2:3566796/1‑150 (MQ=255)
              gTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGc                                                                                                                                        <  2:291572/150‑1 (MQ=255)
              gTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGc                                                                                                                                        <  2:2552785/150‑1 (MQ=255)
                    aCAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGAcc                                                                                                                                  <  2:2678567/150‑1 (MQ=255)
                    aCAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGAcc                                                                                                                                  >  1:1178895/1‑150 (MQ=255)
                      aaaGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCg                                                                                                                                <  2:4693142/150‑1 (MQ=255)
                      aaaGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCg                                                                                                                                <  2:3106492/150‑1 (MQ=255)
                        aGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGaa                                                                                                                              >  1:1444592/1‑150 (MQ=255)
                          ccAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGaaaa                                                                                                                            >  2:2896933/1‑150 (MQ=255)
                            aGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGaaaaaa                                                                                                                          >  1:4797982/1‑149 (MQ=255)
                            aGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAAt                                                                                                                          >  1:5079716/1‑150 (MQ=255)
                            aGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAAt                                                                                                                          >  1:4065674/1‑150 (MQ=255)
                              cTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAAttt                                                                                                                        >  1:1535137/1‑150 (MQ=255)
                               tGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTa                                                                                                                       >  1:2158886/1‑150 (MQ=255)
                                gCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTAc                                                                                                                      >  1:5108445/1‑150 (MQ=255)
                                 ccTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTAct                                                                                                                     <  2:4177189/150‑1 (MQ=255)
                                 ccTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTAct                                                                                                                     <  2:4102026/150‑1 (MQ=255)
                                  cTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTActg                                                                                                                    <  1:319242/150‑1 (MQ=255)
                                   tCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTActgc                                                                                                                   >  2:3423829/1‑150 (MQ=255)
                                     aaCCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTActgctg                                                                                                                 <  1:1864473/150‑1 (MQ=255)
                                       ccGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACtgctggt                                                                                                               <  1:1824392/150‑1 (MQ=255)
                                       ccGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACtgctggt                                                                                                               <  2:4757492/150‑1 (MQ=255)
                                       ccGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAACGAAGCGGTATTGCCTGACCCGAAAAATTTACtgctggt                                                                                                               >  1:1730517/1‑150 (MQ=255)
                                          cAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACtgctggtgct                                                                                                            >  1:2119684/1‑150 (MQ=255)
                                             ttGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGa                                                                                                         <  2:3747440/150‑1 (MQ=255)
                                               gTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACg                                                                                                       <  2:4985617/150‑1 (MQ=255)
                                                   ctactaTCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAgg                                                                                                   >  1:4989761/1‑150 (MQ=255)
                                                     actaTCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGcc                                                                                                 <  2:532168/150‑1 (MQ=255)
                                                            tGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCt                                                                                          >  2:2936843/1‑150 (MQ=255)
                                                              aaTGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGc                                                                                        >  2:3568833/1‑150 (MQ=255)
                                                              aaTGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGc                                                                                        >  1:323270/1‑150 (MQ=255)
                                                                 gCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCgg                                                                                     >  1:4623721/1‑150 (MQ=255)
                                                                  cccGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAACGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGa                                                                                    >  1:3695773/1‑150 (MQ=255)
                                                                   ccGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGAt                                                                                   >  1:546619/1‑150 (MQ=255)
                                                                     gTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGAtgt                                                                                 <  2:1669726/150‑1 (MQ=255)
                                                                         tttGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGgcg                                                                             >  2:1730168/1‑150 (MQ=255)
                                                                          ttGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGgcgc                                                                            >  2:4360010/1‑150 (MQ=255)
                                                                           tGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGcgcg                                                                           <  1:2849889/150‑1 (MQ=255)
                                                                            gTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCgg                                                                          <  1:2461527/150‑1 (MQ=255)
                                                                              cGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGa                                                                        <  2:2050033/150‑1 (MQ=255)
                                                                               gCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGAt                                                                       <  1:4345612/150‑1 (MQ=255)
                                                                                   gtcgGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATgcgc                                                                   >  1:4584255/1‑150 (MQ=255)
                                                                                   gtcgGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATgcgc                                                                   >  1:521605/1‑150 (MQ=255)
                                                                                    tcgGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCt                                                                  >  2:3632443/1‑150 (MQ=255)
                                                                                    tcgGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAACGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCt                                                                  <  1:853469/150‑1 (MQ=255)
                                                                                      gggAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTg                                                                <  2:4132123/150‑1 (MQ=255)
                                                                                      gggAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTg                                                                >  2:2479560/1‑150 (MQ=255)
                                                                                      gggAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAACGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTg                                                                >  2:425376/1‑150 (MQ=255)
                                                                                       ggAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTga                                                               >  1:3090140/1‑150 (MQ=255)
                                                                                       ggAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTga                                                               <  1:2369647/150‑1 (MQ=255)
                                                                                        gagaTTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTgag                                                              <  1:4930503/150‑1 (MQ=255)
                                                                                         agaTTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTgaga                                                             >  2:3968814/1‑150 (MQ=255)
                                                                                          gaTTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGAt                                                            >  2:726568/1‑150 (MQ=255)
                                                                                            ttCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGa                                                          <  2:5082589/150‑1 (MQ=255)
                                                                                             tCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAg                                                         <  1:3395361/150‑1 (MQ=255)
                                                                                               aGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAgcg                                                       <  1:4264955/150‑1 (MQ=255)
                                                                                                     cGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAACGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGaaa                                                 <  2:3296324/150‑1 (MQ=255)
                                                                                                       gAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATc                                               >  1:704204/1‑150 (MQ=255)
                                                                                                        aaGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCa                                              <  1:5085615/150‑1 (MQ=255)
                                                                                                        aaGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCa                                              >  1:1029303/1‑150 (MQ=255)
                                                                                                           tggtggtggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACcg                                           <  1:2102/150‑1 (MQ=255)
                                                                                                           tggtggtggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACcg                                           >  1:3330415/1‑150 (MQ=255)
                                                                                                           tggtggtggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACcg                                           >  1:763928/1‑150 (MQ=255)
                                                                                                               ggtggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGcc                                       <  1:3471851/150‑1 (MQ=255)
                                                                                                               ggtggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGcc                                       <  1:988373/150‑1 (MQ=255)
                                                                                                                 tggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGt                                     >  2:1918202/1‑150 (MQ=255)
                                                                                                                  ggCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTg                                    >  1:1876526/1‑150 (MQ=255)
                                                                                                                   gCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTgg                                   >  1:446426/1‑150 (MQ=255)
                                                                                                                    cAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGt                                  <  1:3665087/150‑1 (MQ=255)
                                                                                                                     aaaCCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTa                                 >  2:2138094/1‑150 (MQ=255)
                                                                                                                     aaaCCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGAGGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTa                                 >  1:4084941/1‑150 (MQ=255)
                                                                                                                       aCCATGCGCCGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATc                               >  2:4815359/1‑150 (MQ=255)
                                                                                                                            gcgcTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTa                          >  1:2002068/1‑150 (MQ=255)
                                                                                                                              gcTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACa                        <  2:1324553/150‑1 (MQ=255)
                                                                                                                               cTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAg                       <  2:4546090/150‑1 (MQ=255)
                                                                                                                                tGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGc                      >  2:3063489/1‑150 (MQ=255)
                                                                                                                                  gTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTg                    >  1:2716498/1‑150 (MQ=255)
                                                                                                                                   tGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTgg                   <  2:2282881/150‑1 (MQ=255)
                                                                                                                                    gATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGa                  >  1:4497017/1‑150 (MQ=255)
                                                                                                                                      tGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACt                >  1:1387148/1‑150 (MQ=255)
                                                                                                                                         cggcgATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGt             <  1:1549655/150‑1 (MQ=255)
                                                                                                                                          ggcgATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCTGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGtt            >  2:3704549/1‑150 (MQ=255)
                                                                                                                                          ggcgATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGtt            <  1:1337505/150‑1 (MQ=255)
                                                                                                                                            cgATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCa          >  1:4736582/1‑150 (MQ=255)
                                                                                                                                            cgATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCa          <  1:1416322/150‑1 (MQ=255)
                                                                                                                                             gATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCAc         >  1:3731860/1‑150 (MQ=255)
                                                                                                                                             gATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCAc         >  1:4210338/1‑150 (MQ=255)
                                                                                                                                             gATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCAc         <  1:440242/150‑1 (MQ=255)
                                                                                                                                              aTGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACg        <  2:1462265/150‑1 (MQ=255)
                                                                                                                                                ggAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGaa      >  2:11261/1‑150 (MQ=255)
                                                                                                                                                ggAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGaa      <  1:4002414/150‑1 (MQ=255)
                                                                                                                                                 gAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGaaa     >  1:3582890/1‑150 (MQ=255)
                                                                                                                                                 gAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGaaa     <  1:4380549/150‑1 (MQ=255)
                                                                                                                                                   aagcgaagcgGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACt   <  2:4488246/150‑1 (MQ=255)
                                                                                                                                                    agcgaagcgGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTg  >  2:1736076/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCTCGTCGGGAGATTCAGGAAGCGGAAGTGGTGGTGGCAAACCATGCGCTGGTGATGGCGGCGATGGAAAGCGAAGCGGTATTGCCTGACCCGAAAAATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTG  >  CP009273/829074‑829371

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: