Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,904,418 C→A 50.4% intergenic (‑62/‑115) kdgR ← / → yebQ KDG regulon transcriptional repressor/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,904,4180CA50.4% 95.0 / inf 229intergenic (‑62/‑115)kdgR/yebQKDG regulon transcriptional repressor/putative transporter
Reads supporting (aligned to +/- strand):  ref base C (61/52);  new base A (63/52);  total (124/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 232 total aligned reads displayed.

TCACCCAGCGCCTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTG  >  CP009273/1904274‑1904565
                                                                                                                                                |                                                                                                                                                   
tCACCCAGCGCCTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAgg                                                                                                                                                  <  2:781124/148‑1 (MQ=255)
  aCCCAGCGCCTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGggt                                                                                                                                                <  2:95854/148‑1 (MQ=255)
     cAGCGCCTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATGTCGCACAAGggtggt                                                                                                                                             >  2:2757023/1‑148 (MQ=255)
           cTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAAt                                                                                                                                       >  1:61311/1‑148 (MQ=255)
                aaaTGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATc                                                                                                                                  >  1:338552/1‑148 (MQ=255)
                  aTGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCtt                                                                                                                                <  2:318924/148‑1 (MQ=255)
                  aTGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCtt                                                                                                                                >  2:1282623/1‑148 (MQ=255)
                  aTGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCtt                                                                                                                                <  1:3162981/148‑1 (MQ=255)
                   tGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCttt                                                                                                                               >  2:3110229/1‑148 (MQ=255)
                         aaaCTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCgg                                                                                                                         >  2:1898802/1‑148 (MQ=255)
                          aaCTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGgt                                                                                                                        <  1:459128/148‑1 (MQ=255)
                           aCTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGgtt                                                                                                                       <  1:1864256/148‑1 (MQ=255)
                             ttttAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGgttag                                                                                                                     <  1:1044506/148‑1 (MQ=255)
                                  gCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACa                                                                                                                >  2:1765005/1‑148 (MQ=255)
                                  gCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACa                                                                                                                >  2:2901144/1‑148 (MQ=255)
                                   cACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACaa                                                                                                               <  2:2039716/148‑1 (MQ=255)
                                      ggAAGATACAGAATCAGGCTGTTTATCCAGATCTGAGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATATAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTc                                                                                                            <  1:3536869/148‑1 (MQ=255)
                                        aaGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTg                                                                                                          >  1:1968449/1‑148 (MQ=255)
                                           aTACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGa                                                                                                       >  2:3099616/1‑148 (MQ=255)
                                           aTACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGa                                                                                                       >  2:2908356/1‑148 (MQ=255)
                                            tACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGAc                                                                                                      <  2:3501748/148‑1 (MQ=255)
                                             aCAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACa                                                                                                     >  1:1585577/1‑148 (MQ=255)
                                               aGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATc                                                                                                   <  2:2274286/148‑1 (MQ=255)
                                                gAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATcc                                                                                                  <  2:2365599/148‑1 (MQ=255)
                                                      ggCTGTTTATCCAGATCTGCGTTAGCCAATTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTg                                                                                            <  2:2473455/148‑1 (MQ=255)
                                                       gCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGa                                                                                           >  2:866606/1‑148 (MQ=255)
                                                        cTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCTGTTGGGTGGTCGCTTTGCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGaa                                                                                          <  1:3490101/148‑1 (MQ=255)
                                                        cTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGaa                                                                                          <  2:1966450/148‑1 (MQ=255)
                                                        cTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGaa                                                                                          <  2:1378149/148‑1 (MQ=255)
                                                        cTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGaa                                                                                          <  1:952621/148‑1 (MQ=255)
                                                         tGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAat                                                                                         <  1:2262982/148‑1 (MQ=255)
                                                         tGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAat                                                                                         <  1:1476578/148‑1 (MQ=255)
                                                          gtatATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTAGCACAAGGGTGGTCGCAATCCATCTTTAGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAata                                                                                        >  2:202162/4‑148 (MQ=255)
                                                          gTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAata                                                                                        <  2:1316164/148‑1 (MQ=255)
                                                           tttATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAatat                                                                                       >  1:2254480/1‑148 (MQ=255)
                                                            ttATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATc                                                                                      >  1:2767974/1‑148 (MQ=255)
                                                             tATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATANTTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCa                                                                                     <  1:1805974/148‑1 (MQ=255)
                                                             tATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCa                                                                                     <  1:2553/148‑1 (MQ=255)
                                                             tATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCa                                                                                     <  1:3136621/148‑1 (MQ=255)
                                                                ccAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTg                                                                                  <  2:1866359/148‑1 (MQ=255)
                                                                ccAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTg                                                                                  <  1:1288925/148‑1 (MQ=255)
                                                                ccAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTg                                                                                  <  1:2427851/148‑1 (MQ=255)
                                                                  aGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCt                                                                                >  2:1928015/1‑148 (MQ=255)
                                                                      cTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGAGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGAGAATATCAGTGCTAGaa                                                                            >  2:2375010/1‑148 (MQ=255)
                                                                        gCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATc                                                                          <  2:80411/148‑1 (MQ=255)
                                                                         cGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCa                                                                         >  1:2853495/1‑148 (MQ=255)
                                                                           ttAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATa                                                                       >  2:362486/1‑148 (MQ=255)
                                                                            tAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATAc                                                                      <  2:839246/148‑1 (MQ=255)
                                                                             aGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATAcc                                                                     >  2:1264792/1‑148 (MQ=255)
                                                                              gCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATAccc                                                                    >  1:2405154/1‑148 (MQ=255)
                                                                               ccATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTCGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATAcccc                                                                   >  1:3042775/1‑148 (MQ=255)
                                                                               ccATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATAcccc                                                                   >  1:1039172/1‑148 (MQ=255)
                                                                                 aTTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTg                                                                 <  1:1062827/148‑1 (MQ=255)
                                                                                    tATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTg                                                              <  2:300987/148‑1 (MQ=255)
                                                                                     aTCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGa                                                             <  1:1545439/148‑1 (MQ=255)
                                                                                         ccTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGattat                                                         <  1:3129120/148‑1 (MQ=255)
                                                                                           tCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTc                                                       >  1:2721011/1‑148 (MQ=255)
                                                                                            cATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCa                                                      >  1:2834465/1‑148 (MQ=255)
                                                                                             aTTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCAc                                                     >  1:1484720/1‑148 (MQ=255)
                                                                                               tGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCa                                                   >  2:829156/1‑148 (MQ=255)
                                                                                                   aGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGa                                               <  2:2994183/148‑1 (MQ=255)
                                                                                                      gtTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGatat                                            <  1:3349068/148‑1 (MQ=255)
                                                                                                       ttttATAAAAATAAGAACAGTTTTTTATTATAATTTAGCACAAGGGTGGTCGCAATCCATCTTTTTCCGGTTAGTTAAAAATATGCGACAACCACCGTGAAAATAAGTGCTAGAATAATAACCCAGTTGATTATTAACCAAAGataaa                                           >  2:50560/1‑146 (MQ=255)
                                                                                                          tataAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATaaaa                                        <  2:3285642/148‑1 (MQ=255)
                                                                                                          tataAAAATCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATaaaa                                        <  2:3228269/148‑1 (MQ=255)
                                                                                                              aaaaTCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTcc                                    >  1:337387/1‑148 (MQ=255)
                                                                                                               aaaTCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCt                                   >  1:920344/1‑148 (MQ=255)
                                                                                                                aaTCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTa                                  <  2:2392824/148‑1 (MQ=255)
                                                                                                                 aTCAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTAt                                 >  1:2053988/1‑148 (MQ=255)
                                                                                                                   cAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGc                               >  1:1338787/1‑148 (MQ=255)
                                                                                                                   cAGAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGc                               >  2:2810817/1‑148 (MQ=255)
                                                                                                                     gAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTTCTAGAATCATACCCATGTTGATTATTCACCAAAGATATAAAATTCCTATGACa                             >  2:253095/1‑148 (MQ=255)
                                                                                                                     gAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCa                             >  1:494397/1‑148 (MQ=255)
                                                                                                                     gAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCa                             >  2:2890341/1‑148 (MQ=255)
                                                                                                                     gAACTGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCa                             >  2:494196/1‑148 (MQ=255)
                                                                                                                      aaCTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTAGAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCaa                            <  1:3484003/148‑1 (MQ=255)
                                                                                                                      aaCTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCaa                            >  2:2176498/1‑148 (MQ=255)
                                                                                                                      aaCTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCaa                            >  2:3063157/1‑148 (MQ=255)
                                                                                                                       ccTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCaaa                           >  2:2249331/2‑148 (MQ=255)
                                                                                                                       aCTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCaaa                           >  2:224429/1‑148 (MQ=255)
                                                                                                                         tGTTTTTTATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAg                         >  1:1085124/1‑148 (MQ=255)
                                                                                                                           ttttttGTTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGtt                       >  2:1625948/1‑148 (MQ=255)
                                                                                                                           ttttttATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGtt                       <  1:2134659/148‑1 (MQ=255)
                                                                                                                            tttttATTATAATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTc                      <  2:1674028/148‑1 (MQ=255)
                                                                                                                             ttttATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCa                     <  2:442749/148‑1 (MQ=255)
                                                                                                                               ttattaTAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAgg                   >  2:1548394/1‑148 (MQ=255)
                                                                                                                                  ttatAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCg                >  2:1526681/1‑148 (MQ=255)
                                                                                                                                   tataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGa               >  2:2723596/1‑148 (MQ=255)
                                                                                                                                    ataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGAc              >  1:2619804/1‑148 (MQ=255)
                                                                                                                                     taATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACg             <  1:1654874/148‑1 (MQ=255)
                                                                                                                                     taATTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACg             >  1:1303739/1‑148 (MQ=255)
                                                                                                                                       aTTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc           <  1:2371102/148‑1 (MQ=255)
                                                                                                                                       aTTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc           >  1:3216134/1‑148 (MQ=255)
                                                                                                                                       aTTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc           >  2:3091118/1‑148 (MQ=255)
                                                                                                                                       aTTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc           <  1:547019/148‑1 (MQ=255)
                                                                                                                                       aTTTCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc           <  2:263838/148‑1 (MQ=255)
                                                                                                                                        tttCGCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGcc          >  1:2810779/1‑148 (MQ=255)
                                                                                                                                          tCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTg        <  1:3099616/148‑1 (MQ=255)
                                                                                                                                            gCACAAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGcc      >  2:409862/1‑148 (MQ=255)
                                                                                                                                             cacaAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCa     >  2:2289044/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCACCCAGCGCCTGCAAAATGCCAAAAACTTTTAGCACGGAAGATACAGAATCAGGCTGTTTATCCAGATCTGCGTTAGCCATTTATCACCTCATTGCGAGTGTTTTATAAAAATCAGAACTGTTTTTTATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTG  >  CP009273/1904274‑1904565

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: