Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,194,160 T→A 5.9% Y135N (TAT→AAT)  yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,194,1600TA 5.9% 689.2 / 26.3 220Y135N (TAT→AAT) yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (92/115);  new base A (5/8);  total (97/123)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 223 total aligned reads displayed.

AAACCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCAGCGTTGCGCC  >  CP009273/2194013‑2194306
                                                                                                                                                   |                                                                                                                                                  
aaaCCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGtt                                                                                                                                                    <  2:2869908/148‑1 (MQ=255)
   cccTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTGTCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATg                                                                                                                                                 >  2:2218863/1‑148 (MQ=255)
   cccTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATg                                                                                                                                                 >  2:1925284/1‑148 (MQ=255)
   cccTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATg                                                                                                                                                 <  2:854796/148‑1 (MQ=255)
    ccTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATga                                                                                                                                                >  1:2522576/1‑148 (MQ=255)
    ccTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTAGTGCGCGGTGTTATga                                                                                                                                                >  2:569882/1‑148 (MQ=255)
     cTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATgag                                                                                                                                               >  2:1599758/1‑148 (MQ=255)
      tCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATgaga                                                                                                                                              >  2:2047215/1‑148 (MQ=255)
       cAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGGGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGagag                                                                                                                                             >  1:1916092/1‑148 (MQ=255)
       cAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTAATGagag                                                                                                                                             <  1:299460/148‑1 (MQ=255)
        aaTCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGt                                                                                                                                            <  1:1947485/148‑1 (MQ=255)
          tCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGa                                                                                                                                          >  1:2394281/1‑148 (MQ=255)
               gcccggACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCctc                                                                                                                                     <  2:478554/148‑1 (MQ=255)
                cccggACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGACAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCctct                                                                                                                                    >  1:2633562/1‑148 (MQ=255)
                 ccggACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTg                                                                                                                                   <  2:830939/148‑1 (MQ=255)
                    gACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTa                                                                                                                                >  1:1296521/1‑148 (MQ=255)
                     aCAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAg                                                                                                                               >  2:1326039/1‑148 (MQ=255)
                      cactGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGt                                                                                                                              <  2:2320170/145‑1 (MQ=255)
                      cAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTAATGAGAGTGATCCTCTGGTAGt                                                                                                                              <  2:1844024/148‑1 (MQ=255)
                       aaTGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGTAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTg                                                                                                                             <  1:132241/148‑1 (MQ=255)
                         tGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGAc                                                                                                                           <  1:1432320/148‑1 (MQ=255)
                          ggcggcGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTAATGAGAGTGATCCTCTGGTAGTGACa                                                                                                                          <  1:3272431/148‑1 (MQ=255)
                               cgcTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCa                                                                                                                     <  1:3235250/148‑1 (MQ=255)
                                  tAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTa                                                                                                                  <  2:2614465/148‑1 (MQ=255)
                                  tAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTa                                                                                                                  <  2:1688703/148‑1 (MQ=255)
                                    gCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATc                                                                                                                >  1:1577575/1‑148 (MQ=255)
                                       tGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAgg                                                                                                             >  2:2498278/1‑148 (MQ=255)
                                         tATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGAc                                                                                                           >  2:1951612/1‑148 (MQ=255)
                                          aTCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACg                                                                                                          <  2:925389/148‑1 (MQ=255)
                                          aTCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACg                                                                                                          <  2:2941786/148‑1 (MQ=255)
                                              cgcAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGa                                                                                                      <  2:1154260/148‑1 (MQ=255)
                                              cgcAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGa                                                                                                      <  2:1154073/148‑1 (MQ=255)
                                                 aaaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTa                                                                                                   <  2:2168225/148‑1 (MQ=255)
                                                 aaaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTa                                                                                                   <  2:3010498/148‑1 (MQ=255)
                                                 aaaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTa                                                                                                   <  1:1124824/148‑1 (MQ=255)
                                                 aaaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTa                                                                                                   <  1:2877739/148‑1 (MQ=255)
                                                 aaaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTa                                                                                                   <  1:207393/148‑1 (MQ=255)
                                                  aaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTAt                                                                                                  <  1:683054/148‑1 (MQ=255)
                                                  aaaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTAt                                                                                                  <  1:800996/148‑1 (MQ=255)
                                                   aaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATg                                                                                                 <  2:1459078/148‑1 (MQ=255)
                                                   aaTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATg                                                                                                 >  2:3310176/1‑148 (MQ=255)
                                                    aTGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATgg                                                                                                >  2:2560226/1‑148 (MQ=255)
                                                     tGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGc                                                                                               <  1:2429679/148‑1 (MQ=255)
                                                      gACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCt                                                                                              >  1:1766641/1‑148 (MQ=255)
                                                       aCAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTa                                                                                             <  1:741865/148‑1 (MQ=255)
                                                       aCAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTa                                                                                             >  2:1937212/1‑148 (MQ=255)
                                                          gCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATgg                                                                                          <  1:675798/148‑1 (MQ=255)
                                                          gCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATgg                                                                                          <  1:2916545/148‑1 (MQ=255)
                                                           ccAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGt                                                                                         >  1:3100161/1‑148 (MQ=255)
                                                           ccAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGAGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGt                                                                                         >  2:2658505/1‑148 (MQ=255)
                                                            cAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGtt                                                                                        <  1:1590414/148‑1 (MQ=255)
                                                            cAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGtt                                                                                        <  2:2079792/148‑1 (MQ=255)
                                                            cAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTAATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGtt                                                                                        <  2:2484036/148‑1 (MQ=255)
                                                               cAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATc                                                                                     >  1:2198289/1‑148 (MQ=255)
                                                                  ttGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCt                                                                                  >  1:1758434/1‑148 (MQ=255)
                                                                       cgacgaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGaa                                                                             >  2:1094453/1‑148 (MQ=255)
                                                                         acgaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTAt                                                                                                  >  1:1059614/1‑125 (MQ=255)
                                                                         acgaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTAt                                                                                                  <  2:1059614/125‑1 (MQ=255)
                                                                         acgaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGGTCAGAAAc                                                                           >  1:905915/1‑148 (MQ=255)
                                                                         acgaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAAc                                                                           <  1:2844601/148‑1 (MQ=255)
                                                                           gaAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCt                                                                         <  2:1715020/148‑1 (MQ=255)
                                                                             aaTCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCtat                                                                       >  1:26793/1‑148 (MQ=255)
                                                                              aTCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCtata                                                                      <  2:1220719/148‑1 (MQ=255)
                                                                               tCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTAATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATaa                                                                     >  2:65489/1‑148 (MQ=255)
                                                                                  tCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGa                                                                  >  2:2357945/1‑148 (MQ=255)
                                                                                  tCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGa                                                                  <  2:3059925/148‑1 (MQ=255)
                                                                                   ccAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGaa                                                                 >  1:1874271/1‑148 (MQ=255)
                                                                                   ccAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGaa                                                                 >  2:3114559/1‑148 (MQ=255)
                                                                                    caacaaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAAt                                                                <  2:1700338/148‑1 (MQ=255)
                                                                                    caacaaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAAt                                                                <  2:1560842/148‑1 (MQ=255)
                                                                                     aacaaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAAtt                                                               <  1:213941/148‑1 (MQ=255)
                                                                                     aacaaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAAtt                                                               >  2:1400732/1‑148 (MQ=255)
                                                                                        aaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGa                                                            <  2:2927774/148‑1 (MQ=255)
                                                                                        aaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGa                                                            <  2:1278475/148‑1 (MQ=255)
                                                                                        aaGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGa                                                            <  1:2659107/148‑1 (MQ=255)
                                                                                          gagGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATTGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGAtc                                                          <  2:2605219/148‑1 (MQ=255)
                                                                                                  ggAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTg                                                  >  1:1466411/1‑148 (MQ=255)
                                                                                                  ggAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTg                                                  >  1:545241/1‑148 (MQ=255)
                                                                                                     atatGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGt                                               >  1:1341200/1‑148 (MQ=255)
                                                                                                      tatGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTa                                              <  1:225713/148‑1 (MQ=255)
                                                                                                       atGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTaa                                             <  2:1804742/148‑1 (MQ=255)
                                                                                                            aCTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCAtc                                        <  2:620794/148‑1 (MQ=255)
                                                                                                               gAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCAtctct                                     >  2:1814361/1‑148 (MQ=255)
                                                                                                                aGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCAtctctc                                    <  1:91136/148‑1 (MQ=255)
                                                                                                                aGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCAtctctc                                    <  1:2554888/148‑1 (MQ=255)
                                                                                                                  gtgttgATTGCACGGCAGTTTATTGCGCGGTTTTATGAGAGTGATCCTCTGGTTGTGAAATTGAAGTATCAGGACGAGGTTTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCAtctctctc                                  <  2:228867/148‑1 (MQ=255)
                                                                                                                   tggtgATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATctctctct                                 >  2:2397265/1‑148 (MQ=255)
                                                                                                                    ggtgATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTg                                >  2:2679063/1‑148 (MQ=255)
                                                                                                                         ttGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTCTAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCaga                           >  2:2115239/1‑148 (MQ=255)
                                                                                                                          tGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCagag                          >  2:114565/1‑148 (MQ=255)
                                                                                                                          tGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCagag                          >  2:2428201/1‑148 (MQ=255)
                                                                                                                                ccAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAgc                    >  1:2571947/1‑148 (MQ=255)
                                                                                                                                   gTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAgctgc                 >  2:943854/1‑148 (MQ=255)
                                                                                                                                     ttATTGCGCGGTGTAATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGctgctg               >  2:1544625/1‑148 (MQ=255)
                                                                                                                                       aTTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGc             <  1:2465533/148‑1 (MQ=255)
                                                                                                                                        ttGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCa            >  2:3231643/1‑148 (MQ=255)
                                                                                                                                        ttGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCa            >  2:1898657/1‑148 (MQ=255)
                                                                                                                                         tGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCAg           >  2:3230658/1‑148 (MQ=255)
                                                                                                                                            gcgGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATGACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCAGCGt        <  1:56379/148‑1 (MQ=255)
                                                                                                                                              ggTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCAGCGTTg      <  2:1297409/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAACCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCTGTATCGCGCAAAATGACAGCCAACAGTTGCTCGACGAAATCGTCCAACAAGAGGGGCTGGAATATGCGACTGAGGTGGTGATTGCACGCCAGTTTATTGCGCGGTGTTATGAGAGTGATCCTCTGGTAGTGACATTGCAGTATCAGGACGAGGATTATGGCTATGGTTATCGCTCAGAAACCTATAACGAATTCGATCTCCGACTGCGTAAGCATCTCTCTCTGGCAGAGGAAAGCTGCTGGCAGCGTTGCGCC  >  CP009273/2194013‑2194306

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: