Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,775,841 T→C 85.3% H228R (CAT→CGT)  clpX ← S14 family endopeptidase ClpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,775,8410TC85.3% 178.0 / 15.9 77H228R (CAT→CGT) clpXS14 family endopeptidase ClpX
Reads supporting (aligned to +/- strand):  ref base T (5/7);  new base C (33/31);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CTTCATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAATTT  >  CP000730/1775710‑1775976
                                                                                                                                   |                                                                                                                                       
cTTCATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTc                                                                                                                                >  1:594019/1‑141 (MQ=255)
cTTCATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTc                                                                                                                                >  1:485732/1‑141 (MQ=255)
cTTCATTGCTTCAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTc                                                                                                                                <  1:531012/141‑1 (MQ=255)
   cATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTACTCACTTCTTCCATACGAACAAAGGCACCACCAAGAATAAATAAAAGATTTGTGGTATCAATTTGAATCATCTAATGGTTTTGATGTTTgtgtaccc                                                                                                                             >  2:206628/1‑137 (MQ=255)
   cATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCCCCAAGAATAAATCAAATATTTGTTGTATCAATTTGAATCAGTTTTTGGTTTGAATGTTTGCGttgaa                                                                                                                             >  2:660812/1‑137 (MQ=255)
   cATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGCATCATTTCTTGGTTTGGATGTTTGCGTCTAc                                                                                                                             >  2:231103/1‑141 (MQ=255)
     ttGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTCGCGTCCACCt                                                                                                                           <  1:502646/141‑1 (MQ=255)
     ttGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCt                                                                                                                           >  2:548297/1‑141 (MQ=255)
     ttGCTTGAGAAACCAATAACTTTATCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCt                                                                                                                           >  2:496688/1‑141 (MQ=255)
     ttGCTTGAGAAACCAATAACTTGTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCt                                                                                                                           >  2:507270/1‑141 (MQ=255)
      tGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCAACtt                                                                                                                          >  1:596162/1‑141 (MQ=255)
        cTTGAGAAACCAATGACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGCTATTTCGTGTATCAAGTTGCAGCATTTCGTGGTTTGGACGTTTGCGTCCACCTtgt                                                                                                                        <  2:665920/141‑1 (MQ=255)
                aCCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAGTCAATTCTTAGTTTGGACGTTTGGGTCCACATATTGGCAGaaa                                                                                                                >  2:201524/1‑140 (MQ=255)
                aCCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAac                                                                                                                >  1:253380/1‑141 (MQ=255)
                aCCAATAACTTTTTCACCAAGACGGCGCTTAATCACCTTTTCAATAACATCAAAGGAACCACAAAGAATAAATAAGATATTTGATGTATACACTTGCATCATGTCTTGCTTTAAACGTTAGCGTCCACCTTGTGGCGGacc                                                                                                                >  2:403786/1‑139 (MQ=255)
                  cAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAacac                                                                                                              <  2:85267/141‑1 (MQ=255)
                    aTAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTAGTAACAATTTGAATCATTTCTTGGTTTGGACGTTAGCGTCCACCTTGTGGCGGAACACtg                                                                                                            >  1:92872/1‑140 (MQ=255)
                        cTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAg                                                                                                        <  1:554917/141‑1 (MQ=255)
                         tttttCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGt                                                                                                       <  2:439428/141‑1 (MQ=255)
                          tttcCACCAAGCCATCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTc                                                                                                      <  2:34030/141‑1 (MQ=255)
                            ttCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGt                                                                                                    >  1:345089/1‑141 (MQ=255)
                            ttCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGt                                                                                                    >  1:40011/1‑141 (MQ=255)
                            ttCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTCGAATCGCTTCTTGGATTGGACGTTTGCGTCCACCTTGTGGCGGAACTGTTGCAGTCGt                                                                                                    >  1:648552/1‑141 (MQ=255)
                             tCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTCGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTa                                                                                                   >  2:576372/1‑141 (MQ=255)
                             tCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTCGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTAGCGGGACACTTGCAGTCTTa                                                                                                   >  2:266626/1‑141 (MQ=255)
                              cACCAAGACGGCGCTTAATCACTTCTTCAATACTATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTTAATCATTTGTTGGTTTGGAAGTTTGAGTTCAACTTGTAGCGGTACAAtttgagttgtta                                                                                                  >  2:450147/1‑131 (MQ=255)
                                ccAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCt                                                                                                <  2:335919/141‑1 (MQ=255)
                                 cAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCtt                                                                                               <  2:462810/141‑1 (MQ=255)
                                        gcgcTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGAt                                                                                        <  2:631326/141‑1 (MQ=255)
                                         cgcTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCAATTCTTTATTTGAACGTTTGTGTCCACCTTGTGGCGGAACACTTCCAGTTGTACCTTCTatgttt                                                                                       >  2:656283/1‑136 (MQ=255)
                                          gcTTAATCACTTCTTCAATACCATCCAAGGCACGCCCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAACAttt                                                                                      <  1:473204/141‑1 (MQ=255)
                                              aaTCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAg                                                                                  >  1:594040/1‑141 (MQ=255)
                                              aaTCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAg                                                                                  >  1:594188/1‑141 (MQ=255)
                                                       ttcAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCa                                                                                <  1:490419/134‑1 (MQ=255)
                                                       ttcAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCa                                                                                >  2:490419/1‑134 (MQ=255)
                                                        tcAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCttgt                                                                        >  2:126714/1‑141 (MQ=255)
                                                        tcAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCttgt                                                                        >  1:602584/1‑141 (MQ=255)
                                                           aTACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGa                                                                     <  1:95657/141‑1 (MQ=255)
                                                           aTACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGa                                                                     <  1:660554/141‑1 (MQ=255)
                                                            tACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGaa                                                                    <  2:573920/141‑1 (MQ=255)
                                                             aCCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTCCCACTTTGAAACATCTCTTGGTTTGGACGTTTGCGTCCACCTGCTGGCGGAATACTTGCAGTCGTACCATCTAAGACTTTAAGCAATACTTGTTGacc                                                                   >  2:96118/1‑139 (MQ=255)
                                                                  cAAAGGCACCACCAAGAATAAATAAGATTTTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCAACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTACGCAATGCTTGTTGAACACCtt                                                              <  1:482353/141‑1 (MQ=255)
                                                                  cAAAGGCACCACCAAAAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCTGAACACTTGCAGTCGTACCTTCTAATATTTTAATCAATGCTTGTTGAACACCtt                                                              <  1:450002/141‑1 (MQ=255)
                                                                   aagggCACCACCGAGGATAAATAAGGGATTTGTTGTATCAATTTGAATCCTTTCTGGGTTTGGACGTGTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTTTGTAATGCTTGTTGAATACCTTc                                                             <  2:301272/138‑1 (MQ=255)
                                                                      ggCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCAcc                                                          >  2:54877/1‑141 (MQ=255)
                                                                      ggCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCAcc                                                          <  1:285199/141‑1 (MQ=255)
                                                                      ggCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCAcc                                                          <  2:330018/141‑1 (MQ=255)
                                                                        caccacCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTg                                                        >  1:290489/1‑141 (MQ=255)
                                                                        caccacCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATCCTTGCGTAACACCTTCGCCTg                                                        >  1:262212/1‑141 (MQ=255)
                                                                          ccaccaAGAATAAATAAGGTATTTCTTGTATCAATTTGAATCATTTCTTTGTTTGGACGTTTGCGGCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGaa                                                      <  1:660850/141‑1 (MQ=255)
                                                                           caccaAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACCCTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGaaa                                                     <  1:507270/141‑1 (MQ=255)
                                                                           caccaAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCGAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCCGaaa                                                     <  1:496688/141‑1 (MQ=255)
                                                                           caccaAGAATAAATAAGAGATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGaaa                                                     <  2:311649/141‑1 (MQ=255)
                                                                               aaGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTc                                                 <  1:54877/141‑1 (MQ=255)
                                                                               aaGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTc                                                 <  1:586516/141‑1 (MQ=255)
                                                                                        taaGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGATCACCTTCACCTGAAACGTCACGTGTtat                                        <  1:441231/141‑1 (MQ=255)
                                                                                           gATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTCGACGTTTGGTTCCAGCGTGTGGCGGAACACTTGCAGTCGTACCTTCTCAGATTTTAAGCCATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGa                                     <  2:72440/141‑1 (MQ=255)
                                                                                               tttGTTGTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGAtgtg                                 <  1:243467/141‑1 (MQ=255)
                                                                                                ttgttgTATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATgtgt                                <  1:548297/141‑1 (MQ=255)
                                                                                                   ttgtATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGtttt                             <  1:161847/141‑1 (MQ=255)
                                                                                                    tgtATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCTCCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTc                            >  2:348887/1‑141 (MQ=255)
                                                                                                    tgtATCAATTTGAATCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAGCGTCACGTGTTATAGATGTGTTTTc                            >  1:167956/1‑141 (MQ=255)
                                                                                                           aTTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCCCCTGAAACGTCACGGGTTATAGAAGTGTTTTCAGATTTa                     >  1:488759/1‑141 (MQ=255)
                                                                                                            tttGAATCATTTCTTGGTTTGGCCGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGCAACGTCAAGTGTTATAGATGTGTTTTCAGATTTAc                    >  2:291018/1‑141 (MQ=255)
                                                                                                                aaTCGTTTCTTGATTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAACCGTCACGTGTTATAGGTGTGTCTTCAGATTTACGTGc                >  1:490460/1‑141 (MQ=255)
                                                                                                                aaTCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGc                >  2:41404/1‑141 (MQ=255)
                                                                                                                aaTCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGc                <  1:532443/141‑1 (MQ=255)
                                                                                                                aaTCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCAACTGAAACGTCACGTGTTATAGATGTGTTTTCAGCTGTACGTGc                >  2:554077/1‑141 (MQ=255)
                                                                                                                  tCATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCaa              >  1:325949/1‑141 (MQ=255)
                                                                                                                   cATTTCTTGGTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTGATAGATGTGTTTTCAGATTTACGTGCatt             >  1:40555/1‑139 (MQ=255)
                                                                                                                   cATTTCTTGGTTTGGACGTTGGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAAt             <  1:295177/141‑1 (MQ=255)
                                                                                                                          tGGTTTGGACGGTTGCGTCCACCTTGTGGCGGACCACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTTAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCa      <  1:35521/141‑1 (MQ=255)
                                                                                                                          tGGTGTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCa      <  1:574784/141‑1 (MQ=255)
                                                                                                                            gTTTGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAAt    >  2:567014/1‑141 (MQ=255)
                                                                                                                             gttGGACGTTTGCGTCCACCTTGTGGCGGAACCCTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAAtt   <  1:150923/140‑1 (MQ=255)
                                                                                                                              ttGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAAttt  <  1:373844/141‑1 (MQ=255)
                                                                                                                              ttGGACGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAAttt  <  2:412261/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CTTCATTGCTTGAGAAACCAATAACTTTTTCACCAAGACGGCGCTTAATCACTTCTTCAATACCATCAAAGGCACCACCAAGAATAAATAAGATATTTGTTGTATCAATTTGAATCATTTCTTGGTTTGGATGTTTGCGTCCACCTTGTGGCGGAACACTTGCAGTCGTACCTTCTAAGATTTTAAGCAATGCTTGTTGAACACCTTCACCTGAAACGTCACGTGTTATAGATGTGTTTTCAGATTTACGTGCAATTTTATCAATTT  >  CP000730/1775710‑1775976

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