Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,942 C→A 5.5% intergenic (+26/‑140) 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,9420CA 5.5% 265.6 / 4.0 82intergenic (+26/‑140)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base A (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

AAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCT  >  NC_002947/289795‑290089
                                                                                                                                                   |                                                                                                                                                   
aaGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCt                                                                                                                                                    <  5:8847/149‑1 (MQ=255)
       tctcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                               >  3:177794/1‑147 (MQ=255)
         tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGa                                                                                                                                             <  6:60841/147‑1 (MQ=255)
         tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGa                                                                                                                                             >  5:60841/1‑147 (MQ=255)
           aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCt                                                                                                                                         >  2:111917/1‑149 (MQ=255)
           aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCt                                                                                                                                         >  6:312637/1‑149 (MQ=255)
             aGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCttt                                                                                                                                       >  3:305743/1‑149 (MQ=255)
              ggCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                      <  5:312637/149‑1 (MQ=255)
              ggCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                      <  1:111917/149‑1 (MQ=255)
              ggCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAAAAGCCTCGCATCAGGCGAGGCttct                                                                                                                                      <  3:299096/149‑3 (MQ=255)
                ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                                    <  4:305743/149‑1 (MQ=255)
                      cTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCt                                                                                                                              >  3:275790/1‑149 (MQ=255)
                            ccTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAt                                                                                                                        <  6:315660/149‑1 (MQ=255)
                                 gTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtata                                                                                                                   <  6:333741/149‑1 (MQ=255)
                                      ggTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                                    <  1:201716/143‑1 (MQ=255)
                                      ggTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                                    >  2:201716/1‑143 (MQ=255)
                                         gtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtttt                                                                                                           >  8:127948/1‑149 (MQ=255)
                                         gtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtttt                                                                                                           >  3:213335/1‑149 (MQ=255)
                                            tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                        >  5:319253/1‑149 (MQ=255)
                                              ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGtt                                                                                                       <  4:213335/148‑1 (MQ=255)
                                              ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTta                                                                                                      >  6:277323/1‑149 (MQ=255)
                                              ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTAAAAAAAAAGCCTCGCATCAGGCGAGGCttctttttca                                                                                                                                >  1:95812/1‑119 (MQ=38)
                                              ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTAAAAAAAAAGCCTCGCATCAGGCGAGGCttctttttca                                                                                                                                <  2:95812/123‑5 (MQ=38)
                                                    tAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTataaat                                                                                                 >  6:213905/1‑148 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTAtaaataaa                                                                                              <  5:222490/149‑1 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTAtaaataaa                                                                                              >  3:248725/1‑149 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTAtaaataaa                                                                                              >  3:321314/1‑149 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTAtaaataaa                                                                                              <  7:127948/149‑1 (MQ=255)
                                                         gcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTg                                                                                           >  4:152393/1‑149 (MQ=255)
                                                         gcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTg                                                                                           >  6:280614/1‑149 (MQ=255)
                                                           gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtt                                                                                         <  4:321314/149‑1 (MQ=255)
                                                           gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtt                                                                                         <  4:248725/149‑1 (MQ=255)
                                                               gcgcCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttattt                                                                                     <  4:177794/149‑1 (MQ=255)
                                                                    cTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCt                                                                                >  3:6000/1‑149 (MQ=255)
                                                                       gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGa                                                                                                                                             >  5:196074/1‑85 (MQ=255)
                                                                       gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGa                                                                                                                                             <  6:196074/85‑1 (MQ=255)
                                                                        aCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttttt                                                                            <  3:265105/149‑1 (MQ=255)
                                                                        aCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttttt                                                                            <  4:275790/149‑1 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtta                                                                                        >  4:129942/1‑135 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtta                                                                                        <  3:129942/135‑1 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                              <  5:73869/145‑1 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                              >  6:73869/1‑145 (MQ=255)
                                                                                 aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                                   <  5:277323/149‑1 (MQ=255)
                                                                                    cTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGa                                                                <  6:319253/149‑1 (MQ=255)
                                                                                       gCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                       <  1:185604/139‑1 (MQ=255)
                                                                                       gCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                       >  2:185604/1‑139 (MQ=255)
                                                                                            aTCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGc                                                        <  3:152393/149‑1 (MQ=255)
                                                                                            aTCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGc                                                        <  5:280614/149‑1 (MQ=255)
                                                                                              cGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGca                                                       >  1:348207/1‑148 (MQ=255)
                                                                                                tGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAt                                                    >  4:15977/1‑149 (MQ=255)
                                                                                                        acTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                              >  7:122076/1‑131 (MQ=255)
                                                                                                        acTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                              <  8:122076/131‑1 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGTTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                  >  4:318012/1‑125 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                  <  3:318012/125‑1 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCAt                                           >  4:1964/1‑148 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                          >  6:301812/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                          >  4:145665/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCCGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                          >  7:83736/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCCGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                          >  5:336968/1‑149 (MQ=255)
                                                                                                            aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACg                                        >  6:227503/1‑149 (MQ=255)
                                                                                                             cAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACgg                                       <  5:227503/149‑1 (MQ=255)
                                                                                                              aaCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                              >  4:337000/1‑125 (MQ=255)
                                                                                                              aaCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                              <  3:337000/125‑1 (MQ=255)
                                                                                                               aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                     <  4:6000/149‑1 (MQ=255)
                                                                                                                cTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTAc                                    <  5:213905/149‑1 (MQ=255)
                                                                                                                 ttCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACa                                   <  3:145665/149‑1 (MQ=255)
                                                                                                                    cTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCCGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGaa                                <  6:336968/149‑1 (MQ=255)
                                                                                                                    cTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCCGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGaa                                <  8:83736/149‑1 (MQ=255)
                                                                                                                        aaGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                              >  6:285363/1‑147 (MQ=255)
                                                                                                                                  tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                  >  1:144260/1‑149 (MQ=255)
                                                                                                                                      aaaaaaGCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAg                   >  2:305283/1‑144 (MQ=255)
                                                                                                                                      aaaaaaGCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAg                   <  1:305283/144‑1 (MQ=255)
                                                                                                                                      aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCAAGGAGCCgtg              >  1:296358/1‑147 (MQ=255)
                                                                                                                                      aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACACATAATGTACAGAACCAAGCCCAGGAGCCg                >  6:169487/1‑147 (MQ=255)
                                                                                                                                        agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGc            <  2:348207/149‑1 (MQ=255)
                                                                                                                                            gCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCAt        >  2:51683/1‑149 (MQ=255)
                                                                                                                                               tCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAg     >  2:144456/1‑149 (MQ=255)
                                                                                                                                                 gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAg     <  3:15977/147‑1 (MQ=255)
                                                                                                                                                 gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGcc   >  2:170209/1‑149 (MQ=255)
                                                                                                                                                 gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGc    <  6:242712/148‑1 (MQ=255)
                                                                                                                                                 gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGc    >  5:242712/1‑148 (MQ=255)
                                                                                                                                                 gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGc    <  1:144456/148‑1 (MQ=255)
                                                                                                                                                   cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCAGCAGCCCATGAGCCt  >  2:6024/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCT  >  NC_002947/289795‑290089

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