Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,946 T→G 7.8% intergenic (+30/‑136) oprE → / → PP_0235 outer‑membrane porin E/peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947289,9460TG 7.8% 217.7 / 7.3 74intergenic (+30/‑136)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base G (4/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCAAA  >  NC_002947/289798‑290093
                                                                                                                                                    |                                                                                                                                                    
gCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGAt                                                                                                                                                      >  6:124487/1‑149 (MQ=255)
  aCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                    >  6:324116/1‑149 (MQ=255)
      tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGACCTGTACATAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGa                                                                                                                                                  <  8:39568/147‑1 (MQ=255)
           ggCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                           <  4:112174/149‑1 (MQ=255)
             ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                                         <  8:189313/149‑1 (MQ=255)
              cTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtg                                                                                                                                        >  3:201658/1‑149 (MQ=255)
                 ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                           >  1:134191/1‑143 (MQ=255)
                 ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                           <  2:134191/143‑1 (MQ=255)
                 ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                                     >  2:15727/1‑149 (MQ=255)
                    ttCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTg                                                                                                                                  >  1:266245/1‑149 (MQ=255)
                        gCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGa                                                                                                                              <  2:210465/149‑1 (MQ=255)
                           tGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAtt                                                                                                                            >  5:1999/1‑148 (MQ=255)
                             cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTTGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAttt                                                                                                                           <  4:186236/147‑1 (MQ=255)
                             cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTta                                                                                                                          >  4:27717/1‑148 (MQ=255)
                              gTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCttttt                                                                                                                                          <  6:239931/131‑1 (MQ=255)
                              gTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCttttt                                                                                                                                          >  5:239931/1‑131 (MQ=255)
                                 aaGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                    >  2:340736/1‑118 (MQ=255)
                                 aaGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                    <  1:340736/118‑1 (MQ=255)
                                    gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgt                                                                                                                                       <  7:16352/128‑1 (MQ=255)
                                    gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgt                                                                                                                                       >  8:16352/1‑128 (MQ=255)
                                         tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                    <  3:264743/110‑1 (MQ=255)
                                         tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                    >  4:264743/1‑110 (MQ=255)
                                         tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAAAAGCCTCGCATCAGGCGAGGCttctttttcag                                                                                                                                    >  5:227585/1‑121 (MQ=38)
                                         tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAAAAGCCTCGCATCAGGCGAGGCttctttttcag                                                                                                                                    <  6:227585/126‑6 (MQ=38)
                                          aCCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                             >  8:187425/1‑148 (MQ=255)
                                           ccGTGGTAGAGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGtt                                                                                                            <  5:282545/148‑1 (MQ=255)
                                            cGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGCATTTTGGTtat                                                                                                          <  6:1999/149‑1 (MQ=255)
                                                gTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTataa‑a                                                                                                      <  1:15727/148‑1 (MQ=255)
                                                    cggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaa‑ataaat                                                                                                 <  1:16251/149‑1 (MQ=255)
                                                      gcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTg                                                                                               >  6:120093/1‑149 (MQ=255)
                                                        gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGtt                                                                                             <  3:166669/149‑1 (MQ=255)
                                                          aaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGttat                                                                                           <  7:187425/149‑1 (MQ=255)
                                                                ccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTc                                                                                     >  1:90045/1‑149 (MQ=255)
                                                                    taCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAATCCTTAGCTGAAAAAGAAGCCTCGCCTGATGCGATGCTTTTTTTTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCttt                                                                                  <  1:91169/147‑1 (MQ=255)
                                                                    gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCttt                                                                                  <  5:272291/148‑1 (MQ=255)
                                                                     aCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCtttt                                                                                 <  2:90045/148‑1 (MQ=255)
                                                                      cAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCATGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTg                                                                               >  4:197817/1‑149 (MQ=255)
                                                                       aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGt                                                                              <  2:266245/149‑1 (MQ=255)
                                                                         aGACCAGACTCGCTTCATCGTCAACTACACTTACAACTTCCTGTAATCCTTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCATATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGttt                                                                            <  3:229411/149‑1 (MQ=255)
                                                                              aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                                       <  5:124487/149‑1 (MQ=255)
                                                                                  tCGGTTCATCGTGAACTACACGTACAACTTGCTGTAAGCCGGAGCTGAAAAAGAAGCCTCGCCGGATGGGAGGCGTTTTTGTGGCTGCAGATTTTTATGTTTTTTTGTTTTTA‑ATTAATTGTTTTTTTTTGTTTTTGTTTTTTATGGAt                                                                   >  8:134258/1‑149 (MQ=255)
                                                                                  tCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                                   >  1:232051/1‑149 (MQ=255)
                                                                                       tCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTACAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCttttt                                                                                >  6:27405/1‑131 (MQ=255)
                                                                                       tCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTACAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCttttt                                                                                <  5:27405/131‑1 (MQ=255)
                                                                                        cATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCa                                                               <  2:244054/147‑1 (MQ=255)
                                                                                            gTGAAGTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGAAGGcaca                                                         >  4:156561/1‑149 (MQ=255)
                                                                                                 cTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGttatt                                                                                          <  5:286532/111‑1 (MQ=255)
                                                                                                 cTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGttatt                                                                                          >  6:286532/1‑111 (MQ=255)
                                                                                                       ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                              <  7:134258/149‑1 (MQ=255)
                                                                                                         aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACg                                            >  2:237395/1‑149 (MQ=255)
                                                                                                             cTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTAc                                        <  5:120093/149‑1 (MQ=255)
                                                                                                                 cTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGaa                                    <  3:27717/149‑1 (MQ=255)
                                                                                                                  tGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAc                                   >  5:69141/1‑149 (MQ=255)
                                                                                                                    tAAGCCTTAGCTTAAAAAGAAGCCTCTCCTGATGCGATGCTTTTTTTTGACTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCGGGCACATTGAACCCATACGGTACAGAccc                                  <  3:201039/148‑1 (MQ=255)
                                                                                                                        ccGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                             <  4:222948/149‑1 (MQ=255)
                                                                                                                             aggcacaaaaaaGCCTCGCATCAGGCGAGGCTTCTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa                        >  4:253989/7‑149 (MQ=38)
                                                                                                                               ttaaaaaCAATCCTCGCCTGATGCGAGGCTTTTTTGTTCCTGCAGATTTATATGTATTTTGGATATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                      <  8:272270/147‑1 (MQ=255)
                                                                                                                               tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAGTCGATTG‑‑GGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCg                    >  7:180621/1‑149 (MQ=255)
                                                                                                                               tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                      >  5:287822/1‑149 (MQ=255)
                                                                                                                               tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                      >  8:185122/1‑149 (MQ=255)
                                                                                                                                 aaaaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                      <  6:69141/147‑1 (MQ=255)
                                                                                                                                    aagaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACa                                       >  7:134745/1‑127 (MQ=255)
                                                                                                                                    aagaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACa                                       <  8:134745/127‑1 (MQ=255)
                                                                                                                                      gaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGcc               >  7:246344/1‑149 (MQ=255)
                                                                                                                                       aagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTGTTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                                  >  5:71486/1‑129 (MQ=255)
                                                                                                                                       aagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTGTTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                                  <  6:71486/129‑1 (MQ=255)
                                                                                                                                       aagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGcc               >  4:308802/1‑148 (MQ=255)
                                                                                                                                          ccTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                      >  3:310258/1‑90 (MQ=255)
                                                                                                                                          ccTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                      <  4:310258/90‑1 (MQ=255)
                                                                                                                                             cGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGc        >  7:253843/1‑149 (MQ=255)
                                                                                                                                              gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGcc       >  2:47852/1‑149 (MQ=255)
                                                                                                                                               ccTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAAGTGTTATTTATTCTTTTTGTTTTTAATAGATTGCAGGCACATTGAACGCATAAGGTAGAGAACCCAGCCGAGGCGCCGCCGCCCATGAGcgt      >  8:132992/1‑147 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAGATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCt      >  1:258579/1‑149 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTc     <  6:287822/149‑1 (MQ=255)
                                                                                                                                                   aTGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCaaa  <  8:253843/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAA‑ATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCAAA  >  NC_002947/289798‑290093

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