Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,429,253 G→T 8.1% pseudogene (13/51 nt) ttcC ← pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,429,2530GT 8.1% 568.6 / 3.4 213pseudogene (13/51 nt)ttcCpseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (105/79);  new base T (2/27);  total (107/106)
Fisher's exact test for biased strand distribution p-value = 1.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

GCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACATTTGC  >  CP009273/1429106‑1429399
                                                                                                                                                   |                                                                                                                                                  
gcaaaaaaAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTg                                                                                                                                                    <  2:3471321/148‑1 (MQ=255)
  aaaaaaaGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGAt                                                                                                                                                  <  2:1173150/148‑1 (MQ=255)
   aaaataGTGTGGCTTATCGTAGCTTTTTTTTGTCGATGATTTTATTTAACTCATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTTTTGTTCTTTCTTACTAATTTCTTGAtt                                                                                                                                                 <  1:1916786/148‑1 (MQ=255)
    aaaaaGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGAttt                                                                                                                                                <  2:1535220/148‑1 (MQ=255)
      aaaGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTg                                                                                                                                              <  1:904734/148‑1 (MQ=255)
       aaGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTATTACTAATTTATTGATTTTGc                                                                                                                                             >  1:1835274/1‑148 (MQ=255)
          agTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTTCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGAc                                                                                                                                          >  2:1442710/1‑148 (MQ=255)
            tGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACAt                                                                                                                                        >  2:900464/1‑148 (MQ=255)
             ggATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACAtt                                                                                                                                       <  2:494724/148‑1 (MQ=255)
             ggATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACAtt                                                                                                                                       <  2:1520948/148‑1 (MQ=255)
                 tATCGTAGCTTATTGTTGTCGCAGATGATAATTAACACATTGAAAAATAAGTAAAATAATTATGAGTAAGAGAGTTGTGATTTTTTCCCTTACTTGTTAAGGTTGTATTGTTCTTTCTTACTAATTTATTTATTTTTCGAAAATTaaa                                                                                                                                   >  1:515387/1‑148 (MQ=255)
                  aTCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTaaaa                                                                                                                                  <  1:2727291/148‑1 (MQ=255)
                     gTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCg                                                                                                                               <  1:2478832/148‑1 (MQ=255)
                      tAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGa                                                                                                                              <  1:2330508/148‑1 (MQ=255)
                         cTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTc                                                                                                                           >  2:1468471/1‑148 (MQ=255)
                            aTTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACGTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAAt                                                                                                                        >  2:1322853/1‑148 (MQ=255)
                              tgttgtCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTc                                                                                                                      >  1:3218897/1‑148 (MQ=255)
                                ttgtCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                                                                                                                    >  1:1108860/1‑148 (MQ=255)
                                ttgtCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                                                                                                                    >  1:2018831/1‑148 (MQ=255)
                                 tgtCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGtt                                                                                                                   >  2:111613/1‑148 (MQ=255)
                                   tCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTtat                                                                                                                 <  2:2788860/148‑1 (MQ=255)
                                    cGCTGATGATATTAAACAAATTCAAAAATAAGAAAAATACTTATGAGTCAGTGAGTTGTGATTTTTGCCCTTACTTGTTCTGTTTGTATTGTTCTTTCTTACTAATTTCTTGTTTTTGCGAAATTTAAAAGCGACTCAATTCGTtata                                                                                                                <  1:3075154/148‑1 (MQ=255)
                                      cTGATTATATTTAACAAATTCAAAAATAATTAAAATAATTATGAGTCAGTGTGTTGTTATTTTTGCCCTTAGTTGTTCTGTTTGTATTGTTATTTCTTACTAATTTCTTTTTTTTGCGACATTTAAACGCGACTCAATTCGTTATATg                                                                                                              <  1:1338724/148‑1 (MQ=255)
                                      cTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATg                                                                                                              >  2:2656693/1‑148 (MQ=255)
                                       ttattatatatAACAAATTCAAAAATAAGAAAAATAATTTTGTGTCAGATTGTTGTTATTTTTGCCATTAGTTGTTCTGTTTGTTTTTTTCTTTCTTACTAATTTCTTTTTTTTGCGAAATTAAAAAGCGACTCAATTCGTTATATgg                                                                                                             <  1:367073/143‑1 (MQ=255)
                                           gatATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATc                                                                                                         >  2:2013271/1‑148 (MQ=255)
                                               tttAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTGTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCagaa                                                                                                     >  1:898223/1‑148 (MQ=255)
                                                ttAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCagaag                                                                                                    >  1:350714/1‑148 (MQ=255)
                                                 tAACAATTTGAAAAATAATTAAAATATTTATGAGTCAGAGTGTTGTTATTTTTGCCCTTACTTTTTCTGGTTTTTTTTTTCTTTATTACTAATTTATTGATTTTGCGAAATTTAAAAGCGACTCAATTCGTTATATGGCATCagaaga                                                                                                   <  1:531834/148‑1 (MQ=255)
                                                   acacATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGGTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATACGTGATATGGCATCAGAAGAGt                                                                                                 >  1:1629657/1‑148 (MQ=255)
                                                       aTTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGc                                                                                             >  2:3056535/1‑148 (MQ=255)
                                                       aTTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGc                                                                                             >  2:1517875/1‑148 (MQ=255)
                                                       aTTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGc                                                                                             >  2:369322/1‑148 (MQ=255)
                                                            aaaaTAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCAt                                                                                        >  1:498543/1‑148 (MQ=255)
                                                             aaaTAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATg                                                                                       >  2:1185086/1‑148 (MQ=255)
                                                              aaTAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGc                                                                                      >  1:420377/1‑148 (MQ=255)
                                                                 aaGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCgg                                                                                   >  2:552616/1‑148 (MQ=255)
                                                                 aaGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCgg                                                                                   >  2:255994/1‑148 (MQ=255)
                                                                  agaaaaaTAATTATGAGTCAGAGTTTTTTTATTTTTGCCCTTACTTTTTCAGGTTTTTTTGTTCTTTATTACTAATTTGTTGATTTTGCCACATTTAAACGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGa                                                                                  <  1:1513437/145‑1 (MQ=255)
                                                                   gTAAAATACTTTTGAGTCAGAGAGTTGTGATTTTTGCCCTTACTGGTTCTGGTTGTTTTGTTCTTTCTTACTATTTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGaa                                                                                 <  1:3127129/148‑1 (MQ=255)
                                                                     aaaaTACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACg                                                                               >  1:1192355/1‑148 (MQ=255)
                                                                     aaaaTACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGAGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACg                                                                               >  1:1192010/1‑148 (MQ=255)
                                                                      aaaTACTTATGAGACAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTTTTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGc                                                                              <  1:1400869/148‑1 (MQ=255)
                                                                          aCTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAg                                                                          >  1:775844/1‑148 (MQ=255)
                                                                           cTTATGAGTCTGAGTGTTTTTTTTTTTGCCCTTACTTGTTCAGGTTTTTTTGTTCTTTGTTACTAATTTCTTGTTTTTGCGAAATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGc                                                                         <  1:3305859/148‑1 (MQ=255)
                                                                            ttATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCa                                                                        <  1:3499752/148‑1 (MQ=255)
                                                                            ttATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCa                                                                        <  2:2609746/148‑1 (MQ=255)
                                                                             tATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCAt                                                                       >  1:999421/1‑148 (MQ=255)
                                                                              ttGAGTCAGAGTTTTGTTATTTTTGCCCTTATTTGTTCTGTTTTTTTTTTTCTTTTTTACTAATTTCTTTTTTTTGATAAATTAAAACGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATa                                                                      <  1:900464/147‑1 (MQ=255)
                                                                               tGAGTCAGAGATTTGTGATTTTTGCCCTTACTTGTTCAGGTTTTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATaa                                                                     <  1:2861911/148‑1 (MQ=255)
                                                                               tGAGACAGTGAGTTGTGATTTTTGCCCTTACTTGTTCTGGTTGTATTTTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATaa                                                                     <  1:3447141/148‑1 (MQ=255)
                                                                                gAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAg                                                                    >  1:1732155/1‑148 (MQ=255)
                                                                                gAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAg                                                                    >  2:712751/1‑148 (MQ=255)
                                                                                 aGTCAGAGTGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGa                                                                   <  2:1127526/148‑1 (MQ=255)
                                                                                  gTCAGAGTGTTTTTATTTTTGCCCTTACTTGTTCTGGTTGTTTTGTTATTTCTTACTAATTTCTTTATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGaa                                                                  <  1:2014009/148‑1 (MQ=255)
                                                                                    cAGAGATTTGTTATTTTTGCCCTTACTTGTTCAGTTTTTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAAt                                                                <  2:1903750/148‑1 (MQ=255)
                                                                                    cAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAAt                                                                >  2:1275883/1‑148 (MQ=255)
                                                                                       agtgttgttTTTTTTGCCGTTACTTTTTCTGGTTGTTTTTTTCTTTCTTACTAATTTCTTTATTTTGCGAAATTAAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTg                                                             <  1:2681838/145‑1 (MQ=255)
                                                                                         tgTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGat                                                           <  2:1284325/147‑1 (MQ=255)
                                                                                           ttGTTATTTTTGCCCTTACTTGTTCAGTTTTTTTTTTTATTTCTAACTATTTTATTTTTTTTGCTAAATTTAAACGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGatat                                                         <  1:1797863/148‑1 (MQ=255)
                                                                                            tgtgATTTTTGCCCTTACTTGTTCTGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGAtata                                                        <  2:2216262/148‑1 (MQ=255)
                                                                                               gATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATaaaa                                                     >  2:1483853/1‑148 (MQ=255)
                                                                                                aTTTTTGCCATTACTTGTTCTGGTTGTTTTTTTCTTTCTTAATATTTTTTTTATTTTGCGACTTTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATaaaaa                                                    <  1:1268272/148‑1 (MQ=255)
                                                                                                 tttttGCCTTTTCTTTTTGTTGTTTTTTTGTTTTTTCTTAATAATTTCTTTTTTTTGCAACATTAAACCGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATaaaaaa                                                   <  1:15565/148‑1 (MQ=255)
                                                                                                 tttttGCCTTTAGTTTTTCTTTTTTTTTTTTTTTTTTTTACTATTTTTTTTTTTTTGCCAAATTAAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATaaaaaa                                                   <  1:1792888/148‑1 (MQ=255)
                                                                                                      gCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACtgtg                                              <  2:2908699/148‑1 (MQ=255)
                                                                                                         cTTACTTGTTCAGGTTGTTTTGTTATTTCTTACTAATTTCTTTATTTTGCGAAATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCg                                           <  1:490952/148‑1 (MQ=255)
                                                                                                         cTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCg                                           >  2:1961936/1‑148 (MQ=255)
                                                                                                             cTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTa                                       >  1:1443894/1‑148 (MQ=255)
                                                                                                               tGTTCAGGTTGTATTGTACTTTCTTACTAATTTCTTTATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATg                                     <  1:576877/148‑1 (MQ=255)
                                                                                                                 ttCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGt                                   <  1:2918075/148‑1 (MQ=255)
                                                                                                                  tCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGtc                                  >  1:2651232/1‑147 (MQ=255)
                                                                                                                  tCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTa                                  <  2:3218897/148‑1 (MQ=255)
                                                                                                                  tCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTa                                  >  1:573712/1‑148 (MQ=255)
                                                                                                                  tCAGGTTGTATTGTTCTTTCTTACTAATTTATTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTa                                  <  2:1237348/148‑1 (MQ=255)
                                                                                                                     ggTTTTATTTTTATTTATTATTAATTTCTTTATTTTGCCAATTTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACgg                               <  1:2800579/148‑1 (MQ=255)
                                                                                                                     ggTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACgg                               >  2:1665795/1‑148 (MQ=255)
                                                                                                                       ttGTTTTTTTCTTTTTAAATATTTTCTTTATTTTGCTAAATTAAAACGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGAt                             <  1:616645/148‑1 (MQ=255)
                                                                                                                        tGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGAtt                            >  2:414160/1‑148 (MQ=255)
                                                                                                                        tGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGAtt                            >  2:2472156/1‑148 (MQ=255)
                                                                                                                          tATTTTTCTTTCTTACTAATTTATTTATTTTGCTACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTag                          <  1:1902063/148‑1 (MQ=255)
                                                                                                                          tATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAACGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTag                          <  1:977444/148‑1 (MQ=255)
                                                                                                                          tATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTag                          <  2:88782/148‑1 (MQ=255)
                                                                                                                            tttttctttctTACTAATTTCTTGATTTTGCGAAATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTagag                        <  1:768403/145‑1 (MQ=255)
                                                                                                                             tGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAgg                       <  2:1701597/148‑1 (MQ=255)
                                                                                                                              gTTCTTTCTTACTAATTTCTTTATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAggg                      <  2:1744684/148‑1 (MQ=255)
                                                                                                                                tctttcttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGAGTGTATGGTACGGATTAGAGGGGa                    >  1:1027578/1‑148 (MQ=255)
                                                                                                                                   ttcttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGaaaa                 <  2:3170144/148‑1 (MQ=255)
                                                                                                                                   ttcttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGaaaa                 <  1:1442710/148‑1 (MQ=255)
                                                                                                                                   ttattAATAATTTCTTTATTTTGCGAAATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGaaaa                 <  1:1009615/145‑1 (MQ=255)
                                                                                                                                     cttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATg               <  2:2693121/148‑1 (MQ=255)
                                                                                                                                     cttACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATg               >  1:1145357/1‑148 (MQ=255)
                                                                                                                                       tACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTc             >  1:2658076/1‑148 (MQ=255)
                                                                                                                                           aaTTTCTTGATTTTGCGACATTAAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGca         <  2:325459/148‑1 (MQ=255)
                                                                                                                                            aTTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGcac        <  2:1027578/148‑1 (MQ=255)
                                                                                                                                             tttCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGcaca       >  2:2141323/1‑148 (MQ=255)
                                                                                                                                              ttCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACAt      <  1:216251/148‑1 (MQ=255)
                                                                                                                                               tCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACAtt     <  2:3354056/148‑1 (MQ=255)
                                                                                                                                                cTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACAttt    <  2:1919910/148‑1 (MQ=255)
                                                                                                                                                  tGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACATTTGc  >  2:3409203/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTGTATGGTACGGATTAGAGGGGAAAATGTCAGCACATTTGC  >  CP009273/1429106‑1429399

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: