Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,158 C→A 23.2% A503A (GCG→GCT)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,038,1580CA23.2% 456.3 / 92.8 185A503A (GCG→GCT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base C (63/79);  new base A (11/32);  total (74/111)
Fisher's exact test for biased strand distribution p-value = 3.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 185 total aligned reads displayed.

AGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGG  >  CP009273/1038012‑1038303
                                                                                                                                                  |                                                                                                                                                 
agagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCg                                                                                                                                                  >  1:473270/1‑148 (MQ=255)
  agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTGTTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                                                                >  2:1938191/1‑148 (MQ=255)
   gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCa                                                                                                                                               <  2:1654453/148‑1 (MQ=255)
     tGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAgg                                                                                                                                             >  1:1665951/1‑148 (MQ=255)
       aaCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGaa                                                                                                                                           >  2:2892281/1‑148 (MQ=255)
       aaCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGaa                                                                                                                                           >  2:2437962/1‑148 (MQ=255)
        aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                                                                          >  1:2740766/1‑148 (MQ=255)
        aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                                                                          >  1:2800082/1‑148 (MQ=255)
            aCAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGa                                                                                                                                      <  2:1754737/148‑1 (MQ=255)
                 cGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttt                                                                                                                                 <  1:740568/148‑1 (MQ=255)
                   ttttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttttt                                                                                                                               <  1:3378315/148‑1 (MQ=255)
                   ttttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttttt                                                                                                                               <  1:3378923/148‑1 (MQ=255)
                    tttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTa                                                                                                                              <  2:1589933/148‑1 (MQ=255)
                     ttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAg                                                                                                                             >  2:590511/1‑148 (MQ=255)
                     ttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAg                                                                                                                             <  2:1487094/148‑1 (MQ=255)
                      tACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGt                                                                                                                            <  1:3213665/148‑1 (MQ=255)
                      tACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGt                                                                                                                            >  2:3323649/1‑148 (MQ=255)
                       aCCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTa                                                                                                                           <  1:1228073/148‑1 (MQ=255)
                         cACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACg                                                                                                                         <  2:3052429/148‑1 (MQ=255)
                            tGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATg                                                                                                                      >  1:1281800/1‑148 (MQ=255)
                             gTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATgg                                                                                                                     <  1:2013542/148‑1 (MQ=255)
                                tGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGgcgc                                                                                                                  <  1:25149/148‑1 (MQ=255)
                                  gCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctg                                                                                                                <  1:2912803/148‑1 (MQ=255)
                                   cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                               <  1:3458972/148‑1 (MQ=255)
                                   cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                               <  1:1235253/148‑1 (MQ=255)
                                      cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                            <  2:1102889/148‑1 (MQ=255)
                                      cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                            <  2:1102938/148‑1 (MQ=255)
                                      cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                            >  2:721598/1‑148 (MQ=255)
                                      cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                            <  2:2648335/148‑1 (MQ=255)
                                       gCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAt                                                                                                           <  2:856864/148‑1 (MQ=255)
                                        cACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAtt                                                                                                          <  1:3033899/148‑1 (MQ=255)
                                        cACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTCTTAGTACGATGGCGCTGCTGAtt                                                                                                          <  1:1262416/148‑1 (MQ=255)
                                              tGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaaa                                                                                                    <  2:3523502/148‑1 (MQ=255)
                                               gATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaat                                                                                                   >  2:2532153/1‑148 (MQ=255)
                                                aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaata                                                                                                  >  2:3291186/1‑148 (MQ=255)
                                                 tCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataa                                                                                                 >  1:980621/1‑148 (MQ=255)
                                                 tCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataa                                                                                                 >  1:1559534/1‑148 (MQ=255)
                                                  cATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaa                                                                                                <  1:1911429/148‑1 (MQ=255)
                                                    ttAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtt                                                                                              <  1:1743397/148‑1 (MQ=255)
                                                     tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                             >  2:2999171/1‑148 (MQ=255)
                                                     tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                             <  1:2892281/148‑1 (MQ=255)
                                                       gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttttt                                                                                           >  2:655814/1‑148 (MQ=255)
                                                        aaTGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTc                                                                                          >  2:838871/1‑148 (MQ=255)
                                                        aaTGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTc                                                                                          <  1:2977765/148‑1 (MQ=255)
                                                          tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                        >  1:2210824/1‑148 (MQ=255)
                                                          tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                        >  2:1637341/1‑148 (MQ=255)
                                                          tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                        <  1:1819490/148‑1 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       >  1:3028502/1‑148 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       >  1:1551063/1‑148 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       <  1:2611310/148‑1 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       >  1:2539356/1‑148 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       <  2:2753085/148‑1 (MQ=255)
                                                           gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                       <  2:63810/148‑1 (MQ=255)
                                                                tctcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                                                                  >  2:2651831/1‑148 (MQ=255)
                                                                    cGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACgg                                                                              <  1:1445290/148‑1 (MQ=255)
                                                                     gTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACggg                                                                             <  2:1664714/148‑1 (MQ=255)
                                                                      tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGt                                                                            <  2:1379134/148‑1 (MQ=255)
                                                                      tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGt                                                                            <  2:3148595/148‑1 (MQ=255)
                                                                      tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGt                                                                            >  1:3484278/1‑148 (MQ=255)
                                                                      tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCTGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGt                                                                            >  2:2854521/1‑148 (MQ=255)
                                                                       ctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTg                                                                           >  2:1413756/1‑148 (MQ=255)
                                                                       ctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTg                                                                           >  2:2379714/1‑148 (MQ=255)
                                                                       ctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTAACGCAGACGGGTg                                                                           >  1:111149/1‑148 (MQ=255)
                                                                          catcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                                                        <  1:908721/148‑1 (MQ=255)
                                                                          catcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                                                        <  1:1819275/148‑1 (MQ=255)
                                                                           atcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGtt                                                                       >  1:2009959/1‑148 (MQ=255)
                                                                               tAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATc                                                                   <  1:143292/148‑1 (MQ=255)
                                                                               tAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATc                                                                   <  1:1938191/148‑1 (MQ=255)
                                                                                 gCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGTAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCa                                                                 >  1:372602/1‑148 (MQ=255)
                                                                                 gCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCa                                                                 >  1:2111530/1‑148 (MQ=255)
                                                                                    aaaTGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGcc                                                              >  2:690944/1‑148 (MQ=255)
                                                                                      aTGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCAc                                                            <  2:2472515/148‑1 (MQ=255)
                                                                                         cAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCg                                                         >  2:3204325/1‑148 (MQ=255)
                                                                                            aCTGGTTCGTGGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                                                      >  2:2110112/1‑148 (MQ=255)
                                                                                             cTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACa                                                     <  2:3018100/148‑1 (MQ=255)
                                                                                                 ttCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTgg                                                 <  2:2009959/148‑1 (MQ=255)
                                                                                                  tCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTggg                                                >  1:637329/1‑148 (MQ=255)
                                                                                                  tCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTggg                                                >  1:2250375/1‑148 (MQ=255)
                                                                                                  tCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTggg                                                <  1:1723603/148‑1 (MQ=255)
                                                                                                    gTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAt                                              <  2:255828/148‑1 (MQ=255)
                                                                                                        cgcACGTACGGCTTCCACAGAAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTg                                          <  1:203579/148‑1 (MQ=255)
                                                                                                         gcACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg                                         <  1:3408321/148‑1 (MQ=255)
                                                                                                          cACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGc                                        <  2:21537/148‑1 (MQ=255)
                                                                                                            cGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCAt                                      <  1:1460714/148‑1 (MQ=255)
                                                                                                               aCGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATaaa                                   <  1:976324/148‑1 (MQ=255)
                                                                                                                cGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAAc                                  >  2:1135952/1‑148 (MQ=255)
                                                                                                                  gCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                <  2:1551102/148‑1 (MQ=255)
                                                                                                                  gCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                <  2:1551063/148‑1 (MQ=255)
                                                                                                                    ttCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTg                              <  1:761291/148‑1 (MQ=255)
                                                                                                                       cacaGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATg                           <  1:173879/148‑1 (MQ=255)
                                                                                                                            cagAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTg                      >  1:1588956/1‑148 (MQ=255)
                                                                                                                              gAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGct                    <  2:473270/148‑1 (MQ=255)
                                                                                                                               aaTCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGctc                   <  1:590511/148‑1 (MQ=255)
                                                                                                                                aTCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGctct                  >  2:1234127/1‑148 (MQ=255)
                                                                                                                                     ccgGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAg             >  2:3354638/1‑148 (MQ=255)
                                                                                                                                     ccgGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAg             <  1:3169145/148‑1 (MQ=255)
                                                                                                                                      cgGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgt            <  2:3513270/148‑1 (MQ=255)
                                                                                                                                       gggTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgtt           <  1:2220543/148‑1 (MQ=255)
                                                                                                                                          ttATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgttgtt        <  2:1705791/148‑1 (MQ=255)
                                                                                                                                                aCCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCgggg  <  1:1244613/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGG  >  CP009273/1038012‑1038303

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: ‑