Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,612,053 G→T 9.2% S28Y (TCC→TAC)  hda ← ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,612,0530GT 9.2% 558.1 / 22.4 158S28Y (TCC→TAC) hdaATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base G (52/91);  new base T (6/9);  total (58/100)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAA  >  CP009273/2611907‑2612200
                                                                                                                                                  |                                                                                                                                                   
  tAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGt                                                                                                                                                  <  1:307593/148‑1 (MQ=255)
    gCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTa                                                                                                                                                <  2:817743/148‑1 (MQ=255)
     cccACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTAt                                                                                                                                               <  2:1303141/148‑1 (MQ=255)
           gCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCg                                                                                                                                         >  2:293644/1‑148 (MQ=255)
            cATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCgg                                                                                                                                        <  1:19484/148‑1 (MQ=255)
            cATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCgg                                                                                                                                        >  1:1904236/1‑148 (MQ=255)
               cGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCa                                                                                                                                     >  1:2241799/1‑148 (MQ=255)
               cGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCa                                                                                                                                     >  1:3246656/1‑148 (MQ=255)
               cGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCa                                                                                                                                     >  2:1804819/1‑148 (MQ=255)
                  cACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGaa                                                                                                                                  <  1:534854/148‑1 (MQ=255)
                   aCGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGaaa                                                                                                                                 <  2:3230277/148‑1 (MQ=255)
                    cGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAAc                                                                                                                                >  2:2000018/1‑148 (MQ=255)
                    cGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAAc                                                                                                                                <  1:1568925/148‑1 (MQ=255)
                        gCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGc                                                                                                                            >  1:1086030/1‑148 (MQ=255)
                         cGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCa                                                                                                                           <  1:2838806/148‑1 (MQ=255)
                            cAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAg                                                                                                                        >  2:1629717/1‑148 (MQ=255)
                             aaTTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAgg                                                                                                                       >  2:3153338/1‑148 (MQ=255)
                             aaTTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAgg                                                                                                                       <  2:2752549/148‑1 (MQ=255)
                              aTTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGt                                                                                                                      <  1:1899028/148‑1 (MQ=255)
                               ttCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGtt                                                                                                                     <  1:387465/148‑1 (MQ=255)
                                 ccGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGACACTTGCAAAGGTTtc                                                                                                                   >  2:2028506/1‑148 (MQ=255)
                                   gcgcAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTtcgt                                                                                                                 <  1:3135809/148‑1 (MQ=255)
                                      cAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTcgtcgt                                                                                                              <  2:875264/148‑1 (MQ=255)
                                       aaGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCgtcgtc                                                                                                             <  2:1835982/148‑1 (MQ=255)
                                        aGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCa                                                                                                            >  2:78467/1‑148 (MQ=255)
                                          ccGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAgg                                                                                                          >  2:2748090/1‑148 (MQ=255)
                                           cgcgTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGa                                                                                                         >  1:364871/1‑148 (MQ=255)
                                           cgcgTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGa                                                                                                         >  1:424552/1‑148 (MQ=255)
                                           cgcgTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGa                                                                                                         >  1:1384726/1‑148 (MQ=255)
                                             cgTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAg                                                                                                       <  2:2257259/148‑1 (MQ=255)
                                              gTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGa                                                                                                      <  2:2682463/148‑1 (MQ=255)
                                               tGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGAt                                                                                                     >  1:2470104/1‑148 (MQ=255)
                                                gcagcaGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATa                                                                                                    <  1:3100666/148‑1 (MQ=255)
                                                gcagcaGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATa                                                                                                    <  1:1661280/148‑1 (MQ=255)
                                                   gcagATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAg                                                                                                 <  1:354572/148‑1 (MQ=255)
                                                    cagATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGt                                                                                                >  2:2592175/1‑148 (MQ=255)
                                                    cagATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGt                                                                                                <  1:2236869/148‑1 (MQ=255)
                                                     agATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTg                                                                                               <  2:1708364/148‑1 (MQ=255)
                                                     agATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTg                                                                                               <  1:990818/148‑1 (MQ=255)
                                                          gCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCaaa                                                                                          >  1:2705979/1‑148 (MQ=255)
                                                          gCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCaaa                                                                                          >  1:2532182/1‑148 (MQ=255)
                                                           cTGCGCCCCGCGCCTTCGCGTGCCCAGAAATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAag                                                                                         >  2:1031651/1‑148 (MQ=255)
                                                            tGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAaga                                                                                        <  1:322208/148‑1 (MQ=255)
                                                            tGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAaga                                                                                        <  1:1988372/148‑1 (MQ=255)
                                                            tGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAaga                                                                                        <  1:2598080/148‑1 (MQ=255)
                                                             gcgcCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAagag                                                                                       <  1:2537243/148‑1 (MQ=255)
                                                               gcCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAagagag                                                                                     <  1:2395934/148‑1 (MQ=255)
                                                                ccccGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGc                                                                                    <  1:2023572/148‑1 (MQ=255)
                                                                 cccGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCt                                                                                   <  2:1736429/148‑1 (MQ=255)
                                                                  ccGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCtg                                                                                  <  2:945008/148‑1 (MQ=255)
                                                                     cgcCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGc                                                                               >  2:3018901/1‑148 (MQ=255)
                                                                     cgcCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGc                                                                               >  2:1725012/1‑148 (MQ=255)
                                                                     cgcCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGGTCCTGACGCAGCACGATTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGc                                                                               <  2:2495688/148‑1 (MQ=255)
                                                                        cTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCgg                                                                            >  2:3307049/1‑148 (MQ=255)
                                                                          tCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGgtg                                                                          >  2:2320168/1‑148 (MQ=255)
                                                                            gcgTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGgtgtg                                                                        <  1:1000053/148‑1 (MQ=255)
                                                                            gcgTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGgtgtg                                                                        <  1:2839602/148‑1 (MQ=255)
                                                                              gTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGtt                                                                      <  2:4834/148‑1 (MQ=255)
                                                                                   cAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCagag                                                                 >  2:2871432/1‑148 (MQ=255)
                                                                                     gagaTAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGAt                                                               >  1:1475885/1‑148 (MQ=255)
                                                                                      agatagatGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATa                                                              >  2:1159934/1‑148 (MQ=255)
                                                                                           gatGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCa                                                         <  2:2470104/148‑1 (MQ=255)
                                                                                           gatGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCa                                                         >  1:1140877/1‑148 (MQ=255)
                                                                                              gTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAAcc                                                      >  1:481491/1‑148 (MQ=255)
                                                                                               tAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCa                                                     <  1:1048173/148‑1 (MQ=255)
                                                                                               tAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCa                                                     <  1:3157374/148‑1 (MQ=255)
                                                                                                aaCCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAg                                                    <  1:1521729/148‑1 (MQ=255)
                                                                                                aaCCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAg                                                    >  1:1723734/1‑148 (MQ=255)
                                                                                                 aCCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAgg                                                   <  1:1209796/148‑1 (MQ=255)
                                                                                                   cGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAATGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGAt                                                 <  2:1610766/148‑1 (MQ=255)
                                                                                                   cGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGAt                                                 <  2:3248951/148‑1 (MQ=255)
                                                                                                     cTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGAttt                                               <  1:1042068/148‑1 (MQ=255)
                                                                                                           tCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACaaa                                         <  1:2514213/148‑1 (MQ=255)
                                                                                                           tCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACaaa                                         <  1:667550/148‑1 (MQ=255)
                                                                                                            ccTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACaaaa                                        <  1:2319725/148‑1 (MQ=255)
                                                                                                             cTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAAt                                       <  2:2062321/148‑1 (MQ=255)
                                                                                                             cTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAAt                                       <  1:3272099/148‑1 (MQ=255)
                                                                                                                 cGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGTAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATcgcg                                   <  2:2856165/148‑1 (MQ=255)
                                                                                                                 cGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATcgcg                                   <  2:63951/148‑1 (MQ=255)
                                                                                                                    agcaCGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCgaga                                <  1:726657/148‑1 (MQ=255)
                                                                                                                     gcaCGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCTTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGaa                               <  2:2241138/148‑1 (MQ=255)
                                                                                                                     gcaCGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGaa                               <  2:2241799/148‑1 (MQ=255)
                                                                                                                      caCGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAg                              >  2:2347883/1‑148 (MQ=255)
                                                                                                                       aCGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGt                             <  2:1904236/148‑1 (MQ=255)
                                                                                                                        cGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGtt                            >  1:1118777/1‑148 (MQ=255)
                                                                                                                        cGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGtt                            >  1:1162927/1‑148 (MQ=255)
                                                                                                                           tttGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTAc                         >  1:2077625/1‑148 (MQ=255)
                                                                                                                           tttGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTAc                         >  1:1422006/1‑148 (MQ=255)
                                                                                                                                   gcgGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAAt                 <  1:960395/148‑1 (MQ=255)
                                                                                                                                      gCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATccc              <  2:866293/148‑1 (MQ=255)
                                                                                                                                       ccAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCa             <  2:1618669/148‑1 (MQ=255)
                                                                                                                                        cAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCat            >  2:3003212/1‑148 (MQ=255)
                                                                                                                                           tAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCatgat         <  1:3002878/148‑1 (MQ=255)
                                                                                                                                           tAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCatgat         <  2:364871/148‑1 (MQ=255)
                                                                                                                                           tAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCatgat         <  2:1261857/148‑1 (MQ=255)
                                                                                                                                            aaaGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCatgatg        <  2:128357/148‑1 (MQ=255)
                                                                                                                                             aaGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGt       >  1:35546/1‑148 (MQ=255)
                                                                                                                                                aggagTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTg    >  1:2923432/1‑148 (MQ=255)
                                                                                                                                                  gagTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGaa  >  1:2913470/1‑148 (MQ=255)
                                                                                                                                                  gagTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGaa  >  1:878923/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCGTGCCCAGAGATAGATGTAACCGCTATGTTCCTGACGCAGCACGTTTTGCAGCGCGGCCAGTAAAGAGGAGTTATCCCCCGGCCAGAAACTTGCAAAGGTTTCGTCGTCAGGAAGATAAAGTGGCAAAGAGAGCTGTGCCGGTGTGTTCAGAGATACCTCAACCAGGATTTCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAA  >  CP009273/2611907‑2612200

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: