Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,944 G→C 7.8% intergenic (+28/‑138) 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,9440GC 7.8% 260.8 / 11.8 90intergenic (+28/‑138)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base G (40/43);  new base C (4/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

AGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATG  >  NC_002947/289796‑290084
                                                                                                                                                    |                                                                                                                                            
aGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTg                                                                                                                                              <  2:303831/149‑1 (MQ=255)
  gCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGAt                                                                                                                                            >  1:42442/1‑149 (MQ=255)
    aCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTGCACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                          >  2:224288/1‑149 (MQ=255)
      tctcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGCCCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                          >  6:203401/1‑147 (MQ=255)
        tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAgg                                                                                                                                      >  7:337839/1‑149 (MQ=255)
        tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAgg                                                                                                                                      >  2:139609/1‑149 (MQ=255)
          aaaGGCGTGGTCTTCAGCCTGCGTCAAGGTTTCTACCGTGTTAGCGGCGAAAGCGCCCTGTACAAAGACCATACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAg                                                                                                                                       <  4:381210/146‑1 (MQ=255)
          aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCt                                                                                                                                    >  2:85555/1‑149 (MQ=255)
            aGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttc                                                                                                                                  <  5:181048/149‑2 (MQ=255)
             ggCCTGGGCTTCAGCCTGCGGCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttct                                                                                                                                 <  5:189327/149‑3 (MQ=255)
               ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               <  6:260136/149‑1 (MQ=255)
                   ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                           <  1:214971/149‑1 (MQ=255)
                   ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                           <  7:81062/149‑1 (MQ=255)
                              gCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTta                                                                                                                >  8:353080/1‑149 (MQ=255)
                               cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAATAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTta                                                                                                                <  1:81961/148‑1 (MQ=255)
                                      gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               <  7:153700/126‑1 (MQ=255)
                                      gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               >  8:153700/1‑126 (MQ=255)
                                      gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtt                                                                                                        <  7:164471/149‑1 (MQ=255)
                                           tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                   >  6:323819/1‑149 (MQ=255)
                                            aCCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTATCTGAAAAAGAATCCTCGACTGATCCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGttt                                                                                                  >  2:149692/1‑149 (MQ=255)
                                             ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTgg                                                                                                    <  1:5245/146‑1 (MQ=255)
                                              cGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTtat                                                                                                <  5:323856/149‑1 (MQ=255)
                                                     gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                               >  4:50946/1‑127 (MQ=255)
                                                     gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                               <  3:50946/127‑1 (MQ=255)
                                                           aaaGCGCCCTGGACAAGGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtta                                                                                   >  1:318894/1‑149 (MQ=255)
                                                            aaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGGAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                  <  1:70668/149‑1 (MQ=255)
                                                               cgcCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttattta                                                                               >  4:238942/1‑149 (MQ=255)
                                                                  ccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTc                                                                            >  5:212281/1‑149 (MQ=255)
                                                                  ccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTc                                                                            >  2:137673/1‑149 (MQ=255)
                                                                      gACAAAGACCAGACTCGCTTCATCGTGAACTGCACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTCGTTATAAATAAATTGTTATTTATTCttt                                                                         <  1:224288/148‑1 (MQ=255)
                                                                      gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCTAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                         <  5:203401/148‑1 (MQ=255)
                                                                      gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                         <  2:320144/148‑1 (MQ=255)
                                                                       aCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttttt                                                                       <  6:148163/149‑1 (MQ=255)
                                                                        cAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                        >  2:363037/1‑131 (MQ=255)
                                                                        cAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                        <  1:363037/131‑1 (MQ=255)
                                                                         aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  3:61358/1‑149 (MQ=255)
                                                                         aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  4:287357/1‑149 (MQ=255)
                                                                         aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  2:125205/1‑149 (MQ=255)
                                                                            gACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                  <  5:323819/149‑1 (MQ=255)
                                                                                aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                              <  3:238942/149‑1 (MQ=255)
                                                                                aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                              <  8:112445/149‑1 (MQ=255)
                                                                                 gACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                             <  1:85555/149‑1 (MQ=255)
                                                                                   cTCGCTTCATCGTGAACAACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTTTTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTGTTTAATcca                                                           >  4:22892/1‑147 (MQ=255)
                                                                                    tCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                          >  5:186632/1‑149 (MQ=255)
                                                                                      gcatCATCGTGACCTACACTTACCACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTg                                                        <  8:253552/146‑1 (MQ=255)
                                                                                        ttCATCGTGAACTACACTTAAAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                         <  1:4437/146‑1 (MQ=255)
                                                                                          cATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAg                                                     >  2:223182/1‑148 (MQ=255)
                                                                                          cATCGTGAACTACACTTACAACATCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCa                                                      <  8:335441/147‑1 (MQ=255)
                                                                                           aTCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAgg                                                    <  8:337839/148‑1 (MQ=255)
                                                                                              gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                 <  3:372086/68‑1 (MQ=255)
                                                                                              gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGcac                                                 <  3:287357/148‑1 (MQ=255)
                                                                                              gTGAACTAAACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                 >  4:372086/1‑68 (MQ=255)
                                                                                               tGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAt                                               >  1:29051/1‑149 (MQ=255)
                                                                                                 aaCTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAtt                                              <  2:29051/148‑1 (MQ=255)
                                                                                                   cTACACTTACCACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGaa                                           <  2:333873/149‑1 (MQ=255)
                                                                                                      cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGACTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATTTATTTTGGTTATAAATAAATTGTTATTTATTATTTTTTTTTTTAATAGATTGAAGGAACATTGaaac                                         >  7:221339/1‑146 (MQ=255)
                                                                                                      caattACAACTTACTGTAAGCCGTATCTGAAAAAGAAGCATCGCCTGATGAGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGTTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAcc                                         <  8:221339/145‑1 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  8:3212/1‑149 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  6:359009/1‑149 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  7:226024/1‑149 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  7:211476/1‑149 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     <  1:137673/149‑1 (MQ=255)
                                                                                                         ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  5:187967/1‑149 (MQ=255)
                                                                                                           aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACg                                   >  5:145312/1‑149 (MQ=255)
                                                                                                             aaCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACgg                                  >  7:105520/1‑148 (MQ=255)
                                                                                                              aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                <  4:61358/149‑1 (MQ=255)
                                                                                                                ttCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                >  1:312137/1‑147 (MQ=255)
                                                                                                                      tAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                         >  2:295484/1‑148 (MQ=255)
                                                                                                                            gTATCTGAAAAAGAAGCCTCGACTGATGCGATGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  5:351521/149‑1 (MQ=255)
                                                                                                                            gTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  8:105520/149‑1 (MQ=255)
                                                                                                                            gTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  2:318894/149‑1 (MQ=255)
                                                                                                                             tAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcccag                 <  2:42442/149‑1 (MQ=255)
                                                                                                                               gCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  6:187967/149‑1 (MQ=255)
                                                                                                                               gCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAAAAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  1:125205/149‑1 (MQ=255)
                                                                                                                               aggcacaaaaaaGCCTCGCATCAGGCGAGGCTTCTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcccag                 >  4:98545/7‑147 (MQ=38)
                                                                                                                                 tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGACTGCAGCTTTATATGTATTTTGGTTATAAATTAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAGCCCAGCCGAGGAGc             >  6:329742/1‑149 (MQ=255)
                                                                                                                                  gAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGcc            >  7:223977/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCACAGGAGCCg           >  5:352010/1‑147 (MQ=255)
                                                                                                                                     aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  2:272171/140‑1 (MQ=255)
                                                                                                                                     aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  >  1:272171/1‑140 (MQ=255)
                                                                                                                                      aagaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAg        <  1:223182/149‑1 (MQ=255)
                                                                                                                                       agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                    <  2:92682/136‑1 (MQ=255)
                                                                                                                                       agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                    >  1:92682/1‑136 (MQ=255)
                                                                                                                                       agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGc       >  7:281398/1‑149 (MQ=255)
                                                                                                                                       agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGACGCAGc       >  8:7711/1‑149 (MQ=255)
                                                                                                                                       aaaagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGc       >  4:188049/3‑149 (MQ=255)
                                                                                                                                          agCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGccc     <  7:353080/148‑1 (MQ=255)
                                                                                                                                           gCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCAt   <  6:186632/149‑1 (MQ=255)
                                                                                                                                           gCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  2:256333/137‑1 (MQ=255)
                                                                                                                                           gCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAAATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               >  1:256333/1‑137 (MQ=255)
                                                                                                                                             cTCGCCTGATGCGATGCTTTTTTGTGCATGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATg  <  1:149692/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATG  >  NC_002947/289796‑290084

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