Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,922,164 C→A 72.2% P311P (CCG→CCT)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,922,1640CA72.2% 392.5 / 299.3 328P311P (CCG→CCT) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base C (43/48);  new base A (124/113);  total (167/161)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 330 total aligned reads displayed.

ACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAAC  >  CP009273/1922015‑1922310
                                                                                                                                                     |                                                                                                                                                  
     ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAAt                                                                                                                                               <  2:1831105/150‑1 (MQ=255)
      gCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAAtt                                                                                                                                              <  2:3669740/150‑1 (MQ=255)
       ccGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTa                                                                                                                                             <  1:4619068/150‑1 (MQ=255)
        cGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTaa                                                                                                                                            >  1:1098735/1‑150 (MQ=255)
        cGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTaa                                                                                                                                            >  1:2928454/1‑150 (MQ=255)
         ggATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAAc                                                                                                                                           <  1:4811125/150‑1 (MQ=255)
         ggATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAAc                                                                                                                                           >  1:4511786/1‑150 (MQ=255)
            tcatcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAActct                                                                                                                                        >  2:5247512/1‑150 (MQ=255)
            tcatcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAActct                                                                                                                                        <  1:492523/150‑1 (MQ=255)
               tcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGAAGTTTTCGCGTGGAGGAATTAACTCATca                                                                                                                                     >  2:1274059/1‑149 (MQ=255)
                caAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTccc                                                                                                                                    <  2:2268437/150‑1 (MQ=255)
                  aaGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCAc                                                                                                                                  >  2:3728935/1‑150 (MQ=255)
                      cAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGtt                                                                                                                              >  1:5128488/1‑150 (MQ=255)
                       aaTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTg                                                                                                                             >  2:4257146/1‑150 (MQ=255)
                        aTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGc                                                                                                                            >  2:1951956/1‑150 (MQ=255)
                         tACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCt                                                                                                                           >  1:511432/1‑150 (MQ=255)
                          aCCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCtc                                                                                                                          >  2:2134167/1‑150 (MQ=255)
                          aCCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCtc                                                                                                                          >  2:4766174/1‑150 (MQ=255)
                             aaTGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCCGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAACTCTTCCCACTGTTGCtcatc                                                                                                                       >  2:258047/1‑150 (MQ=255)
                             aaTGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGAAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAACTCTTCCCACTGTTGCtcatc                                                                                                                       >  2:5049621/1‑150 (MQ=255)
                                gACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACg                                                                                                                    <  1:4241949/150‑1 (MQ=255)
                                   ttCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCAt                                                                                                                 <  1:2750498/150‑1 (MQ=255)
                                          gTCTTGCTGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTa                                                                                                          <  1:2591387/150‑1 (MQ=255)
                                          gTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTa                                                                                                          <  1:2992435/150‑1 (MQ=255)
                                           tCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTAc                                                                                                         <  1:4399255/150‑1 (MQ=255)
                                            cTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGATGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACg                                                                                                        <  2:2631871/150‑1 (MQ=255)
                                            cTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACg                                                                                                        <  1:1303965/150‑1 (MQ=255)
                                            cTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACg                                                                                                        >  1:1793122/1‑150 (MQ=255)
                                                cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtat                                                                                                    <  2:1606809/150‑1 (MQ=255)
                                                  gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATaa                                                                                                  >  2:2511528/1‑150 (MQ=255)
                                                  gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATaa                                                                                                  >  1:726915/1‑150 (MQ=255)
                                                  gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATaa                                                                                                  >  1:4964427/1‑150 (MQ=255)
                                                   ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAg                                                                                                 <  1:5225611/150‑1 (MQ=255)
                                                   ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAg                                                                                                 >  1:4644171/1‑150 (MQ=255)
                                                   ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAg                                                                                                 <  1:4696350/150‑1 (MQ=255)
                                                           gcgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGtt                                                                                         <  2:1375712/150‑1 (MQ=255)
                                                           gcgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGtt                                                                                         <  2:1718234/150‑1 (MQ=255)
                                                              cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGttttt                                                                                      <  2:3066285/150‑1 (MQ=255)
                                                              cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGttttt                                                                                      <  1:3593038/150‑1 (MQ=255)
                                                               aaaCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTg                                                                                     >  2:4144930/1‑150 (MQ=255)
                                                                    ggTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTg                                                                                <  1:3159643/150‑1 (MQ=255)
                                                                     gTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTgg                                                                               <  2:3292668/150‑1 (MQ=255)
                                                                     gTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTgg                                                                               >  1:1990267/1‑150 (MQ=255)
                                                                       tAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCg                                                                             <  2:4396160/150‑1 (MQ=255)
                                                                       tAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCg                                                                             <  1:2654932/150‑1 (MQ=255)
                                                                       tAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCg                                                                             >  1:1646549/1‑150 (MQ=255)
                                                                           ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                         >  1:2172716/1‑150 (MQ=255)
                                                                           ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                         <  2:1513372/150‑1 (MQ=255)
                                                                              cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                      <  1:1401175/150‑1 (MQ=255)
                                                                                cTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACg                                                                    <  1:1628595/150‑1 (MQ=255)
                                                                                cTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACg                                                                    >  2:4114454/1‑150 (MQ=255)
                                                                                cTGACGCTCTTCAACCACCAGCCAGCCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACg                                                                    <  2:2695324/150‑1 (MQ=255)
                                                                                   aCGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAg                                                                 <  2:1037607/150‑1 (MQ=255)
                                                                                     gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                               <  2:4024751/150‑1 (MQ=255)
                                                                                     gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                               >  1:631487/1‑150 (MQ=255)
                                                                                     gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                               <  1:226896/150‑1 (MQ=255)
                                                                                       tctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaa                                                             <  2:1442047/150‑1 (MQ=255)
                                                                                       tctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaa                                                             <  2:391982/150‑1 (MQ=255)
                                                                                        ctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGAGGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                            >  2:2737318/1‑150 (MQ=255)
                                                                                         ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGGGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                           >  2:4943032/1‑150 (MQ=255)
                                                                                           cAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                         <  2:823163/150‑1 (MQ=255)
                                                                                            aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACg                                                        >  1:326272/1‑150 (MQ=255)
                                                                                            aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACg                                                        <  1:876561/150‑1 (MQ=255)
                                                                                              ccaccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                      <  1:4100661/150‑1 (MQ=255)
                                                                                              ccaccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                      <  2:3191768/150‑1 (MQ=255)
                                                                                                 ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACTATGCt                                                   <  1:3727219/150‑1 (MQ=255)
                                                                                                    gccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGt                                                <  1:440120/150‑1 (MQ=255)
                                                                                                          cGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGAt                                          >  2:5202249/1‑150 (MQ=255)
                                                                                                           ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATc                                         >  1:1266215/1‑150 (MQ=255)
                                                                                                           ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATc                                         <  2:1669636/150‑1 (MQ=255)
                                                                                                            gTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCa                                        >  1:3118067/1‑150 (MQ=255)
                                                                                                             tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                       >  2:2826208/1‑150 (MQ=255)
                                                                                                             tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                       <  2:113129/150‑1 (MQ=255)
                                                                                                             tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                       >  2:2889208/1‑150 (MQ=255)
                                                                                                              aaaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                      <  2:1176856/150‑1 (MQ=255)
                                                                                                              aaaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                      >  2:4445749/1‑150 (MQ=255)
                                                                                                              aaaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                      >  2:4196124/1‑150 (MQ=255)
                                                                                                               aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                     <  1:2378601/150‑1 (MQ=255)
                                                                                                                  aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                  >  2:3539978/1‑150 (MQ=255)
                                                                                                                  aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                  <  1:69960/150‑1 (MQ=255)
                                                                                                                    cGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTa                                >  2:4895220/1‑150 (MQ=255)
                                                                                                                     gTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCTGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                               <  2:236118/150‑1 (MQ=255)
                                                                                                                     gTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                               >  1:1386614/1‑150 (MQ=255)
                                                                                                                      tAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAtt                              <  1:3851565/150‑1 (MQ=255)
                                                                                                                       aaaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                             >  1:1570287/1‑150 (MQ=255)
                                                                                                                       aaaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                             <  2:1793122/150‑1 (MQ=255)
                                                                                                                        aaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCg                            >  1:3414984/1‑150 (MQ=255)
                                                                                                                           ccTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                         >  2:3397646/1‑150 (MQ=255)
                                                                                                                            cTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGTTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGAt                        >  2:1978349/1‑150 (MQ=255)
                                                                                                                            cTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGAt                        <  2:3881407/150‑1 (MQ=255)
                                                                                                                             ttCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATc                       <  1:1174317/150‑1 (MQ=255)
                                                                                                                                cAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                    <  2:466018/150‑1 (MQ=255)
                                                                                                                                cAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                    <  1:4257146/150‑1 (MQ=255)
                                                                                                                                cAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                    <  2:510702/150‑1 (MQ=255)
                                                                                                                                 aGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTCGTGATCAAGGCTGTATTCGTGATCTTTg                   >  2:4615224/1‑150 (MQ=255)
                                                                                                                                   cATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcg                 >  1:1019068/1‑150 (MQ=255)
                                                                                                                                        tgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg            >  2:788725/1‑150 (MQ=255)
                                                                                                                                                cgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGAGATCTTTGCGGCGCGGCAGaaaa    >  1:452180/1‑150 (MQ=255)
                                                                                                                                                cgTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAAATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaa    >  2:680230/1‑150 (MQ=255)
                                                                                                                                                  tGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTGGCCGTGGCGGTTGGAACGCAGATACAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac  <  2:884236/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAAC  >  CP009273/1922015‑1922310

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