Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 330,730 T→C 5.2% intergenic (+21/‑63) esxA → / → esaA ESAT‑6 family virulence protein/virulence protein EsaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730330,7300TC 5.2% 223.9 / 3.0 78intergenic (+21/‑63)esxA/esaAESAT‑6 family virulence protein/virulence protein EsaA
Reads supporting (aligned to +/- strand):  ref base T (46/27);  new base C (2/2);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

CTAAAGTAGAAAAATTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGATATTTTTT  >  CP000730/330594‑330870
                                                                                                                                        |                                                                                                                                            
cTAAAGTAGAAAAATTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTcaca                                                                                                                                          >  2:355110/1‑141 (MQ=255)
      tAGAAAAATTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGCCACAttagaa                                                                                                                                    <  2:370697/141‑5 (MQ=255)
          aaaaTTTGCACAATTATTAGAAGAAATTCAACAACAAATGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAAtg                                                                                                                                >  1:417367/1‑141 (MQ=255)
          aaaaTTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAAtg                                                                                                                                >  1:325556/1‑141 (MQ=255)
          aaaaTTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGCCACAttagaaaatg                                                                                                                                <  1:509593/141‑9 (MQ=255)
            aaTTTGCACAATTATAAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAAtgtg                                                                                                                              <  1:538171/141‑1 (MQ=255)
               ttGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACCACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCt                                                                                                                           <  1:254626/141‑1 (MQ=255)
               ttGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCt                                                                                                                           >  1:765926/1‑141 (MQ=255)
                  tacaATTATTAGAAGAANTTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTg                                                                                                                        >  1:492290/2‑141 (MQ=255)
                  cacaATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTg                                                                                                                        <  2:779012/141‑1 (MQ=255)
                  cacaATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTg                                                                                                                        >  1:448774/1‑141 (MQ=255)
                    caATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCt                                                                                                                      >  2:597260/1‑141 (MQ=255)
                    caATTATNAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGGAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCt                                                                                                                      <  2:589764/141‑1 (MQ=255)
                        tattaGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATc                                                                                                                  <  1:101253/141‑1 (MQ=255)
                           taGAAGAAATTAAACAACAATTGAATAGCACTGCTGAAGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCAtt                                                                                                               >  1:567016/1‑141 (MQ=255)
                            agaagaAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCAttt                                                                                                              <  1:778239/141‑1 (MQ=255)
                              aagaaATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCAttttt                                                                                                            >  1:105606/1‑141 (MQ=255)
                              aagaaATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCAttttt                                                                                                            >  1:473915/1‑141 (MQ=255)
                               agaaATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCAtttttt                                                                                                           <  1:429936/141‑1 (MQ=255)
                                 aaaTTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTa                                                                                                         <  1:99867/141‑1 (MQ=255)
                                   aTTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAg                                                                                                       >  1:84518/1‑141 (MQ=255)
                                      aaaCAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGaaa                                                                                                    >  1:212512/1‑141 (MQ=255)
                                       aacaacAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGaaaa                                                                                                   >  1:393241/1‑141 (MQ=255)
                                       aacaacAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGaaaa                                                                                                   >  1:494235/1‑141 (MQ=255)
                                        acaacaATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAaac                                                                                                  <  2:494235/141‑1 (MQ=255)
                                            caaTTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTGATCATTTTTTTAAGAAAACAACt                                                                                              >  2:118213/1‑141 (MQ=255)
                                             aaTTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTg                                                                                             <  1:741918/141‑1 (MQ=255)
                                               ttGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGaa                                                                                           >  1:758284/1‑141 (MQ=255)
                                                tGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGaaa                                                                                          >  2:256335/1‑141 (MQ=255)
                                                 gAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAg                                                                                         <  1:585104/141‑1 (MQ=255)
                                                 gAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGACACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAg                                                                                         <  2:260422/141‑1 (MQ=255)
                                                   aTAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGa                                                                                       >  2:82823/1‑141 (MQ=255)
                                                     aGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGaaa                                                                                     >  2:571074/1‑141 (MQ=255)
                                                       cACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATa                                                                                   >  2:55565/1‑141 (MQ=255)
                                                          tgctgATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGc                                                                                >  1:87396/1‑141 (MQ=255)
                                                          tgctgATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGc                                                                                >  2:597480/1‑141 (MQ=255)
                                                          tgctgATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGc                                                                                >  2:763291/1‑141 (MQ=255)
                                                                  ccGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGAAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAATAAAACAACTGAAAGGAAATAAGCATgaaaaa                                                                        >  2:687851/1‑141 (MQ=255)
                                                                     ttCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATgaaaaagaa                                                                     >  2:395750/1‑141 (MQ=255)
                                                                      tCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATgaaaaagaaa                                                                    >  2:666992/1‑141 (MQ=255)
                                                                       cAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAAAAACTGAAAGGAAATAAGCATgaaaaagaaaa                                                                   >  1:234342/1‑141 (MQ=255)
                                                                         aGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGAAAAAAAAAAAAAAAt                                                                 <  2:10661/141‑1 (MQ=255)
                                                                          gAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAAtt                                                                >  2:721280/1‑141 (MQ=255)
                                                                           aaCAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTg                                                               >  1:175739/1‑141 (MQ=255)
                                                                             cAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGa                                                             >  1:758324/1‑141 (MQ=255)
                                                                              aaGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGAt                                                            >  2:378718/1‑141 (MQ=255)
                                                                                gACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGAttt                                                          >  2:455737/1‑141 (MQ=255)
                                                                                 aCCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTa                                                         >  1:59233/1‑141 (MQ=255)
                                                                                 aCCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTa                                                         >  2:524997/1‑141 (MQ=255)
                                                                                  ccAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTAt                                                        <  1:607172/141‑1 (MQ=255)
                                                                                    aacaacTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGc                                                      >  1:216644/1‑141 (MQ=255)
                                                                                      caacTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCAt                                                    <  1:230476/141‑1 (MQ=255)
                                                                                      caacTTTATAATAATAACGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGATTATCATTTTTTTAAGAAAACAGCTGAATGGAAATAAGAAAAAAAAAAAAAAATTGGATTTATGCAt                                                    <  2:20481/141‑1 (MQ=255)
                                                                                         cTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTaa                                                 >  1:358327/1‑141 (MQ=255)
                                                                                         cTTTCTAATAATGTCGGTTTTTAATAAGGAGTCTGAAAGTGGGATAGTGACTGTGTCTAATGTGGCTGTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTaa                                                 <  1:662097/141‑1 (MQ=255)
                                                                                              taataaTTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGAAAAAAAAAAAAAAATTGGATTTATGCATTAATTGTc                                            <  2:124620/141‑1 (MQ=255)
                                                                                               aataatTTCGGTTTGCAAAAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAAGTGTCa                                           <  2:676212/141‑1 (MQ=255)
                                                                                                 taatTTCGGTTTGAAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAAAAAAATTGGATTTATGCATTAATTGTCACt                                         <  2:693937/141‑1 (MQ=255)
                                                                                                         gTTTGTAATAAGCATTCTGAAATTGGAAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGTTTTATGCATTAATTGTCACTTTAattat                                 <  2:375381/141‑1 (MQ=255)
                                                                                                         gTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAattat                                 >  1:535523/1‑141 (MQ=255)
                                                                                                         gTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAattat                                 <  2:393241/141‑1 (MQ=255)
                                                                                                          tttGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAattatt                                <  1:140427/141‑1 (MQ=255)
                                                                                                             gCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATtata                             >  2:701444/1‑141 (MQ=255)
                                                                                                                 tAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTg                         <  1:414358/141‑1 (MQ=255)
                                                                                                                  aaGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTTATTTTACTACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGc                        >  2:554621/1‑141 (MQ=255)
                                                                                                                     cATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCAAAAAAAAAAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCAt                     <  2:87396/141‑1 (MQ=255)
                                                                                                                       ttCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCAtag                   >  2:571989/1‑141 (MQ=255)
                                                                                                                         cTGAAATTGGCAAAGTCAAATTTTCTAATGTGGCTTTGCTTATCATTTTTTTACGAACACAACTGAAAGGAAATAAGCAAGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCAtagtt                 >  2:180309/1‑141 (MQ=255)
                                                                                                                          tGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCAtagtta                >  2:456113/1‑141 (MQ=255)
                                                                                                                              ataagCAAAGCCACATTAGAAAATGTGACTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCAGAGTTAGTAt            >  1:517451/5‑141 (MQ=38)
                                                                                                                              aTTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTAt            <  2:280146/141‑1 (MQ=255)
                                                                                                                                tGGCAAAGTAAAATTTTATAATGTGGCTTTGCTTATCATTTTTTTAATAAAACAATTTAAATGAAATAAGCATAAAAAAGAAAAATTGGAGTTATGCATTAATTATGACTGTAATTATAAAAATTGCCATGGTTAGTATGa          >  2:90932/1‑141 (MQ=255)
                                                                                                                                  gCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATAAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGata        <  2:373454/141‑1 (MQ=255)
                                                                                                                                   cAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAGACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATATTTAGTATGatat       >  2:423526/1‑141 (MQ=255)
                                                                                                                                   cAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGatat       >  1:40695/1‑141 (MQ=255)
                                                                                                                                   cAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGatat       <  2:175739/141‑1 (MQ=255)
                                                                                                                                   cAAAGCCACATTAGAAAATGTGACTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGctat       >  2:761993/1‑141 (MQ=255)
                                                                                                                                      aGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGATAtttt    >  1:596432/1‑141 (MQ=255)
                                                                                                                                        tCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGATAtttttt  <  1:529171/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
CTAAAGTAGAAAAATTTGCACAATTATTAGAAGAAATTAAACAACAATTGAATAGCACTGCTGATGCCGTTCAAGAACAAGACCAACAACTTTCTAATAATTTCGGTTTGCAATAAGCATTCTGAAATTGGCAAAGTCACATTTTCTAATGTGGCTTTGCTTATCATTTTTTTAAGAAAACAACTGAAAGGAAATAAGCATGAAAAAGAAAAATTGGATTTATGCATTAATTGTCACTTTAATTATTATAATTGCCATAGTTAGTATGATATTTTTT  >  CP000730/330594‑330870

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