Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 338,121 T→A 10.6% H753Q (CAT→CAA)  essC → virulence protein EssC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730338,1210TA10.6% 171.5 / 8.0 68H753Q (CAT→CAA) essCvirulence protein EssC
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base A (4/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGTGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTTGTATC  >  CP000730/337981‑338261
                                                                                                                                            |                                                                                                                                            
ggTGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCat                                                                                                                                              >  1:254931/1‑141 (MQ=255)
ggTGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCat                                                                                                                                              >  1:186545/1‑141 (MQ=255)
ggTGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACGTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCat                                                                                                                                              >  2:248969/1‑141 (MQ=255)
 gTGGGGGTATGGCGAAGTTAATTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCTAATTGATAAATCGTCAACGTTTATTCGGAGAGCatg                                                                                                                                             <  2:643251/141‑1 (MQ=255)
  tGGGTGTTTGGGGAACTTATTTAAAGATTTAGTCCATTTAGTTGGGACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCatga                                                                                                                                            <  1:215954/141‑1 (MQ=255)
  tGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCatga                                                                                                                                            <  1:427196/141‑1 (MQ=255)
  tGGGGGTATGGCGAACTTATTTAAAGATGTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCatga                                                                                                                                            <  1:464244/141‑1 (MQ=255)
     gggTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGt                                                                                                                                         <  2:304249/141‑1 (MQ=255)
     gggTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGt                                                                                                                                         <  1:337242/141‑1 (MQ=255)
       gTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTa                                                                                                                                       >  2:423573/1‑141 (MQ=255)
       gTATGGCGAACTTATTTAAAGATCTAGTCCATTTAGTAGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAACGCCGAATTCAGAAAACGCCAACGTTTATTCGAGGAGCAAGAGGTTa                                                                                                                                       >  2:386077/1‑141 (MQ=255)
         aTGTCGAACTTATTTAAAGATTTAGTCCATTTAGTTTGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAAc                                                                                                                                     <  2:195341/141‑1 (MQ=255)
         aTGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAAc                                                                                                                                     >  1:275124/1‑141 (MQ=255)
          tGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAAcc                                                                                                                                    <  1:117564/141‑1 (MQ=255)
           ggcgaaacaaTTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCa                                                                                                                                   <  1:337602/132‑1 (MQ=255)
            gCGAACTTAGTTAAAGATTTAGTTCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCat                                                                                                                                  <  1:576276/141‑1 (MQ=255)
               aaCTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGACGCGATGCGTGCCTTAACATCAATCAAAAGGCAGGCGACAAAACGACAAAGTTTAGTCAGGTATCATGATGTTAACCATAtt                                                                                                                               >  1:315385/1‑141 (MQ=255)
               aaCTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATAtt                                                                                                                               <  1:68456/141‑1 (MQ=255)
                    aTTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAAGTGAGAAACCGGCAACGTTTATTCGGAGCGCATGAGGTTAACCATATTAATCa                                                                                                                          >  2:382657/1‑141 (MQ=255)
                     tttAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCaa                                                                                                                         >  1:456597/1‑141 (MQ=255)
                             tttAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATaa                                                                                                                 >  1:608620/1‑141 (MQ=255)
                             tttAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATaa                                                                                                                 <  1:317212/141‑1 (MQ=255)
                                  tCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCGGAATTGAGAAAACGACACAGTTTATTCGGAGAGCACGATGTTAAGCATATTAATCAATACCATAAGTTAt                                                                                                            >  2:30873/1‑141 (MQ=255)
                                  tCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTAt                                                                                                            >  1:291893/1‑141 (MQ=255)
                                   ccATTTAGTTGGTACGATTACACACTTAGATGGCGATGAAGCGATGCGTGACTTAACCTCAAACAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTGAACCATATTAATCAATACCATAAGTTAtt                                                                                                           <  2:659246/141‑1 (MQ=255)
                                         aGTTGGTACGATTGCAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAaaga                                                                                                     <  1:551701/141‑1 (MQ=255)
                                         aGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAaaga                                                                                                     >  1:501879/1‑141 (MQ=255)
                                         aGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAaaga                                                                                                     >  1:501905/1‑141 (MQ=255)
                                         aGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAaaga                                                                                                     >  1:17179/1‑141 (MQ=255)
                                          gTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAaagaa                                                                                                    <  1:385104/141‑1 (MQ=255)
                                             ggTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGt                                                                                                 >  1:40018/1‑141 (MQ=255)
                                               tACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTAt                                                                                               >  2:13836/1‑141 (MQ=255)
                                                     tACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGAc                                                                                         <  2:430297/141‑1 (MQ=255)
                                                     tACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGAc                                                                                         <  2:430300/141‑1 (MQ=255)
                                                       cAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAg                                                                                       >  1:59148/1‑141 (MQ=255)
                                                         aaCTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGaa                                                                                     <  1:212977/141‑1 (MQ=255)
                                                               gATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATg                                                                               >  1:546506/1‑141 (MQ=255)
                                                               gATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATg                                                                               >  1:528041/1‑141 (MQ=255)
                                                                   gCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCcac                                                                           >  2:106618/1‑141 (MQ=255)
                                                                      aTGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACAtt                                                                        <  2:212561/141‑1 (MQ=255)
                                                                       tGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCCCAttt                                                                       >  2:668172/1‑141 (MQ=255)
                                                                        gAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACCTTTa                                                                      >  2:500240/1‑141 (MQ=255)
                                                                          aGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTAtt                                                                    <  2:9721/141‑1 (MQ=255)
                                                                           gCGGTGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTc                                                                   <  1:312169/141‑1 (MQ=255)
                                                                           gCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGTGCATGATGTTAACCATATTAATCAATACCATAATTTACTTGGATAATGTATTGCGACAGAACCAATGCCACATTTATTc                                                                   <  2:275124/141‑1 (MQ=255)
                                                                           gCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATGAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCcacaagtatcg                                                                   >  2:300527/1‑134 (MQ=255)
                                                                           gCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTc                                                                   >  1:645891/1‑141 (MQ=255)
                                                                           gCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATAACATAAGTTATTTAAAGAAGGTATTGCGACAGAACACATGACACATCTATTc                                                                   >  2:464321/1‑141 (MQ=255)
                                                                                    gCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTcc                                                          >  2:88561/1‑141 (MQ=255)
                                                                                       ttAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCCTTATTTCCGAt                                                       >  1:124276/1‑141 (MQ=255)
                                                                                        tAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATg                                                      <  2:648852/141‑1 (MQ=255)
                                                                                        tAACATAAATCAAAGCCGAATTGAGAAGACGGCAACTTTTATCCGTAGAGCATGATGTTAACCATATGAAGCAAGACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATg                                                      <  2:580901/141‑1 (MQ=255)
                                                                                        tAACACCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTGAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATg                                                      <  2:211354/141‑1 (MQ=255)
                                                                                           cATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATCCCATAAGTTATTTGAAGAAGTTATTTCGACAGCACCGATGCCACATTTATTCATTATTTCCGATGAGt                                                   <  2:140974/141‑1 (MQ=255)
                                                                                             tcaatcaaAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGttt                                                 >  2:535165/1‑141 (MQ=255)
                                                                                                atcaaAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGcc                                              <  2:608620/141‑1 (MQ=255)
                                                                                                 tcaaAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCg                                             <  1:401264/141‑1 (MQ=255)
                                                                                                        cGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTaaa                                      >  1:423732/1‑141 (MQ=255)
                                                                                                            ttGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCa                                  <  1:467267/141‑1 (MQ=255)
                                                                                                                gaAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGaac                              >  1:331011/1‑141 (MQ=255)
                                                                                                                    aCGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAAcc                          <  1:560941/141‑1 (MQ=255)
                                                                                                                    aCGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACCTTTATTCATGATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAAcc                          <  2:501879/141‑1 (MQ=255)
                                                                                                                    aCGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAAcc                          <  2:17179/141‑1 (MQ=255)
                                                                                                                    aCGTCAACGTTTATTCGGAGAGCAAGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAAcc                          <  2:501905/141‑1 (MQ=255)
                                                                                                                     cGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCt                         <  1:382657/141‑1 (MQ=255)
                                                                                                                        cAACGTTTATTCGGAGAGCATGATGCTAACCATATGAATCAAATCCATAAGTTATTTAAAGAATGTATTGCGATAGAACCAATGCCACATTCATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGAt                      <  1:302372/141‑1 (MQ=255)
                                                                                                                         aaCGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGGACAACCTGAtt                     >  1:8373/1‑141 (MQ=255)
                                                                                                                          aCGTTTATTGGGGGAGCAGGAGGTTAACCATATTAATCAATACCATAAGTTATTTAAGGAAGGTACTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGAttt                    >  2:505160/1‑141 (MQ=255)
                                                                                                                          aCGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGAttt                    >  1:255349/1‑141 (MQ=255)
                                                                                                                            tgttATTCGGAAGGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTa                  <  2:399432/139‑1 (MQ=255)
                                                                                                                               tctacGGACAGCAGGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATGGCGACAGAACCAATGCCACAGTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATAGGACCAGACTGAttgcagga               >  2:576868/5‑135 (MQ=255)
                                                                                                                                  tCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAg            <  1:70042/141‑1 (MQ=255)
                                                                                                                                   cGGAGAGCAGGATGTTAACCATATTAATCAATACCATAAGTTATTTAGAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGa           >  2:430062/1‑141 (MQ=255)
                                                                                                                                          cATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTTGTa    >  1:410684/1‑141 (MQ=255)
                                                                                                                                          cATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTTGTa    >  2:606139/1‑141 (MQ=255)
                                                                                                                                           atgatgTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTAGtcc   >  2:328537/1‑139 (MQ=255)
                                                                                                                                            tgatgTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTTGTATc  >  1:360942/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GGTGGGGGTATGGCGAACTTATTTAAAGATTTAGTCCATTTAGTTGGTACGATTACAAACTTAGATGGCGATGAAGCGATGCGTGCCTTAACATCAATCAAAGCCGAATTGAGAAAACGTCAACGTTTATTCGGAGAGCATGATGTTAACCATATTAATCAATACCATAAGTTATTTAAAGAAGGTATTGCGACAGAACCAATGCCACATTTATTCATTATTTCCGATGAGTTTGCCGAATTAAAATCAGAACAACCTGATTTTATGAAAGAACTTGTATC  >  CP000730/337981‑338261

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑