Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 158,967 C→A 19.0% H126N (CAT→AAT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273158,9670CA19.0% 790.3 / 146.2 293H126N (CAT→AAT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (112/125);  new base A (20/36);  total (132/161)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 302 total aligned reads displayed.

CAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATC  >  CP009273/158818‑159116
                                                                                                                                                     |                                                                                                                                                     
cAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTc                                                                                                                                                       <  1:4964697/150‑1 (MQ=255)
cAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGAATGAACCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGACTTTc                                                                                                                                                       >  1:1890862/1‑150 (MQ=255)
 aGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCa                                                                                                                                                      <  1:48923/150‑1 (MQ=255)
 aGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCa                                                                                                                                                      <  1:2221862/150‑1 (MQ=255)
    cGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGa                                                                                                                                                   <  1:301588/150‑1 (MQ=255)
    cGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGa                                                                                                                                                   <  1:3250349/150‑1 (MQ=255)
    cGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCGTGATGAATTTAATGa                                                                                                                                                   >  1:1444916/1‑150 (MQ=255)
          cGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAg                                                                                                                                             >  1:3727091/1‑150 (MQ=255)
              ggCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTg                                                                                                                                         >  2:3885279/1‑150 (MQ=255)
               gCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGc                                                                                                                                        >  2:3753687/1‑150 (MQ=255)
                cTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCa                                                                                                                                       >  1:1820194/1‑150 (MQ=255)
                 tACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAg                                                                                                                                      <  1:2042632/150‑1 (MQ=255)
                 tACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAg                                                                                                                                      >  1:3796443/1‑150 (MQ=255)
                 tACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAg                                                                                                                                      >  1:2748125/1‑150 (MQ=255)
                  aCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAgg                                                                                                                                     <  2:4919593/150‑1 (MQ=255)
                    cGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCg                                                                                                                                   >  2:4224138/1‑150 (MQ=255)
                     ggATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCgg                                                                                                                                  >  2:3574692/1‑150 (MQ=255)
                     ggATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTAGCAGGCgg                                                                                                                                  >  2:1383736/1‑150 (MQ=255)
                     ggATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAAGTTCATGAGCGCAGCTTGCAGGCgg                                                                                                                                  <  2:1944803/150‑1 (MQ=255)
                      gATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGa                                                                                                                                 <  2:4999214/150‑1 (MQ=255)
                      gATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGa                                                                                                                                 >  2:3351485/1‑150 (MQ=255)
                       aTGCGTGCGCAACACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGAt                                                                                                                                <  2:3113438/150‑1 (MQ=255)
                       aTGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGAt                                                                                                                                >  1:1065451/1‑150 (MQ=255)
                        tgcgtgcgCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGAtt                                                                                                                               <  1:462694/150‑1 (MQ=255)
                           gtgcgCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttgg                                                                                                                            <  2:2642061/150‑1 (MQ=255)
                           gtgcgCAAAACTGCGTAGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttgg                                                                                                                            <  2:189760/150‑1 (MQ=255)
                            tgcgCAACACGGCGTCGGGCCGAAGACCCGCCTGGAAGTGGTTACCGACGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATttggc                                                                                                                           <  2:3852423/150‑1 (MQ=255)
                            tgcgCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggc                                                                                                                           <  2:4509510/150‑1 (MQ=255)
                            tgcgCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggc                                                                                                                           <  1:3183472/150‑1 (MQ=255)
                            tgcgCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggc                                                                                                                           <  2:497378/150‑1 (MQ=255)
                               gcAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgtt                                                                                                                        <  2:1641834/150‑1 (MQ=255)
                               gcAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgtt                                                                                                                        <  2:3023014/150‑1 (MQ=255)
                                 aaatCTGCGTAGGGCCGTCTACCAGCCTGGAAGGGGTTACCGAAGGCGTGCTGACTCGAATGATCCAGCGTGACCCGGAACTGATGGGTGTTGTTCTTGTGTTCCTTGATGTATTTAATGAGCGCAGCTTGCAGGCGGATttggcgttgg                                                                                                                      <  1:5008395/150‑1 (MQ=255)
                                 aaaaCTGCGTCGGGCCGAATACCCGCCTGTAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgttgg                                                                                                                      >  1:4546722/1‑150 (MQ=255)
                                  aaaCTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATttggcgttggc                                                                                                                     >  1:4445215/1‑150 (MQ=255)
                                     cTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCt                                                                                                                  <  2:2057664/150‑1 (MQ=255)
                                     cTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCt                                                                                                                  <  2:1394283/150‑1 (MQ=255)
                                           cGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACt                                                                                                            >  2:1461517/1‑150 (MQ=255)
                                            gggCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTc                                                                                                           >  1:1285094/1‑150 (MQ=255)
                                            gggCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTc                                                                                                           >  1:2552552/1‑150 (MQ=255)
                                               ccGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAt                                                                                                        >  2:2409977/1‑150 (MQ=255)
                                                cGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAtg                                                                                                       <  1:307517/150‑1 (MQ=255)
                                                 gAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAtgt                                                                                                      <  1:2277220/150‑1 (MQ=255)
                                                  aaTACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAtgtg                                                                                                     >  2:3460819/1‑150 (MQ=255)
                                                      cccGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaac                                                                                                 >  1:2844838/1‑150 (MQ=255)
                                                        cGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaacaa                                                                                               >  1:1599578/1‑150 (MQ=255)
                                                          ccTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAgg                                                                                             >  1:4371467/1‑150 (MQ=255)
                                                          ccTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAgg                                                                                             <  1:3450063/150‑1 (MQ=255)
                                                            tGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTc                                                                                           >  2:584633/1‑150 (MQ=255)
                                                                  tGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtg                                                                                     >  1:924921/1‑150 (MQ=255)
                                                                   ggTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtga                                                                                    <  2:3028209/150‑1 (MQ=255)
                                                                   ggTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTAGTCGATGTGCAACAAGGTCTGCGtga                                                                                    >  2:1091451/1‑150 (MQ=255)
                                                                     ttACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtgatg                                                                                  <  1:4316982/150‑1 (MQ=255)
                                                                     ttACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtgatg                                                                                  <  1:1539955/150‑1 (MQ=255)
                                                                       aCCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGAc                                                                                <  1:2403333/150‑1 (MQ=255)
                                                                       aCCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGAc                                                                                >  2:3617384/1‑150 (MQ=255)
                                                                       aCCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGAc                                                                                <  2:2758247/150‑1 (MQ=255)
                                                                        ccGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGAcc                                                                               >  2:2104771/1‑150 (MQ=255)
                                                                          gAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCtt                                                                             <  2:2420562/150‑1 (MQ=255)
                                                                             ggCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTaaa                                                                          <  1:209786/150‑1 (MQ=255)
                                                                             ggCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTaaa                                                                          >  1:4974163/1‑150 (MQ=255)
                                                                               cGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAAct                                                                        <  1:5188337/150‑1 (MQ=255)
                                                                                     gACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGAt                                                                  >  2:3250032/1‑150 (MQ=255)
                                                                                      aCGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGAtt                                                                 <  1:3939259/150‑1 (MQ=255)
                                                                                            aTGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCg                                                           >  2:3949683/1‑150 (MQ=255)
                                                                                            aTGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCg                                                           >  2:5222978/1‑150 (MQ=255)
                                                                                             tGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCgg                                                          >  2:4158352/1‑150 (MQ=255)
                                                                                              gATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGc                                                         >  2:8848/1‑150 (MQ=255)
                                                                                               aTCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCt                                                        <  2:2257174/150‑1 (MQ=255)
                                                                                                tCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTa                                                       >  2:687360/1‑150 (MQ=255)
                                                                                                  cAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTAcc                                                     <  1:1189096/150‑1 (MQ=255)
                                                                                                  cAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTAcc                                                     <  1:629870/150‑1 (MQ=255)
                                                                                                   aGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTAccc                                                    >  2:816873/1‑150 (MQ=255)
                                                                                                   aGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTAccc                                                    <  1:1946053/150‑1 (MQ=255)
                                                                                                    gCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCt                                                   <  2:2457816/150‑1 (MQ=255)
                                                                                                     cGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTg                                                  >  1:4759225/1‑150 (MQ=255)
                                                                                                         aCCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACa                                              >  2:5175707/1‑150 (MQ=255)
                                                                                                            cGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACg                                           >  1:168820/1‑150 (MQ=255)
                                                                                                            cGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACg                                           <  2:2691082/150‑1 (MQ=255)
                                                                                                              gAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGAc                                         >  1:3139675/1‑150 (MQ=255)
                                                                                                               aaCTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGAcc                                        <  1:764152/150‑1 (MQ=255)
                                                                                                                     gCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgc                                  <  2:1215367/150‑1 (MQ=255)
                                                                                                                     gCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgc                                  >  2:776049/1‑150 (MQ=255)
                                                                                                                      cGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgca                                 <  2:2936153/150‑1 (MQ=255)
                                                                                                                      cGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgca                                 <  2:4439360/150‑1 (MQ=255)
                                                                                                                       ggTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgcag                                <  2:4098908/150‑1 (MQ=255)
                                                                                                                        gtgtTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTgcagc                               >  1:2938328/1‑150 (MQ=255)
                                                                                                                           ttGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCaaa                            <  1:408259/150‑1 (MQ=255)
                                                                                                                                cTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAAtgctg                       >  2:4338847/1‑150 (MQ=255)
                                                                                                                                cTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAAtgctg                       >  1:2181499/1‑150 (MQ=255)
                                                                                                                                  ggTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGcc                     <  1:1475125/150‑1 (MQ=255)
                                                                                                                                     gATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGa                  <  2:3437260/150‑1 (MQ=255)
                                                                                                                                     gATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGa                  >  2:2786800/1‑150 (MQ=255)
                                                                                                                                      aTCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGaa                 >  1:3957547/1‑150 (MQ=255)
                                                                                                                                      aTCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGaa                 <  2:2673814/150‑1 (MQ=255)
                                                                                                                                       tCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAg                >  1:85111/1‑150 (MQ=255)
                                                                                                                                         cTTGAGGAATTTCATGAGCGCAGCTTGCAGGCGGATATGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTTATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAgcg              <  2:2814553/150‑1 (MQ=255)
                                                                                                                                         cTTGACGAATTTCTCGCGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACACCGACCGCTTGCAGCAAATGCTGCCAGAAgcg              <  2:3674526/150‑1 (MQ=255)
                                                                                                                                             atgaATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTg          <  1:594302/150‑1 (MQ=255)
                                                                                                                                                     cATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGTCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCAtc  <  2:4839810/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATC  >  CP009273/158818‑159116

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: