Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,944 G→C 5.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 5.8% 213.1 / 4.1 69intergenic (+28/‑138)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTC  >  NC_002947/289799‑290090
                                                                                                                                                  |                                                                                                                                                  
cACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATg                                                                                                                                                 >  8:70446/1‑149 (MQ=255)
 aCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                >  7:259010/1‑149 (MQ=255)
 aCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                                >  8:327090/1‑149 (MQ=255)
       aaaGGCCTGGGCTTCAGCCTGCGTCACGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACTCTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAg                                                                                                                                             <  1:90221/146‑1 (MQ=255)
       aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAGGACCAGAGGCGGTTCATCGTGAAGTACCGTTACAAGTTCGTGTAGGCCGTAGGTGAAAAAGAAGCGTCGCGTGATGCGAGGCt                                                                                                                                          >  8:225912/1‑149 (MQ=255)
       aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCt                                                                                                                                          >  1:66923/1‑149 (MQ=255)
         aGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCttt                                                                                                                                        >  3:29695/1‑149 (MQ=255)
           gCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTTGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTCGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                       <  1:254590/149‑1 (MQ=255)
            ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCTCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                                     <  5:364525/149‑1 (MQ=255)
             gtGGGCTTCAGCCTGCGTCAAGGTGTCTACCTTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtg                                                                                                                                    <  8:283036/148‑1 (MQ=255)
              tGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAg                                                                                                                                             >  7:196604/1‑139 (MQ=255)
              tGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAg                                                                                                                                             <  8:196604/139‑1 (MQ=255)
                ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgt                                                                                                                                   <  7:327090/147‑1 (MQ=255)
                       gCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttct                                                                                                                                       >  3:272641/1‑134 (MQ=255)
                       gCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttct                                                                                                                                       <  4:272641/136‑3 (MQ=255)
                        ccTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAt                                                                                                                         >  8:197007/1‑149 (MQ=255)
                         cTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAtt                                                                                                                        >  4:249862/1‑149 (MQ=255)
                            cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                                     <  4:29695/149‑1 (MQ=255)
                             gTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCGTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtata                                                                                                                    >  6:370467/1‑149 (MQ=255)
                                 atgggtCTACCTGTGTATCGGCGAAACCGCCCTTGTCCAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAACCCGTAGCTTAAAAACAATCCTCGCCTGATGCGAGGCTTTCTTGTGCCTGCAGATTTATATGt                                                                                                                 <  7:225912/144‑1 (MQ=255)
                                   gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtt                                                                                                              >  7:35900/1‑149 (MQ=255)
                                      tCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTg                                                                                                           <  7:309970/149‑1 (MQ=255)
                                             tGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTataa                                                                                                    <  8:35900/149‑1 (MQ=255)
                                               gTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTataaat                                                                                                  <  8:259010/149‑1 (MQ=255)
                                                   cggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                              <  3:154488/149‑1 (MQ=255)
                                                    ggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAAtt                                                                                             >  4:68511/1‑149 (MQ=255)
                                                       gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtt                                                                                          >  6:268835/1‑149 (MQ=255)
                                                       gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGt                                                                                           >  8:246279/1‑148 (MQ=255)
                                                        aaaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                              >  5:20045/1‑144 (MQ=255)
                                                        aaaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                              <  6:20045/144‑1 (MQ=255)
                                                         aaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTAGACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                        >  5:251019/1‑149 (MQ=255)
                                                                   gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTGGCCTGATGCGAGGCTTTTTTGTGCCTGGAGATTTATATGTATTTTGGTTTTAAATAAATTGTTATTTATTCtttt                                                                              <  2:66923/149‑1 (MQ=255)
                                                                          aCCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                        >  2:42517/1‑148 (MQ=255)
                                                                            cAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTaa                                                                     <  5:268835/149‑1 (MQ=255)
                                                                             aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                                    <  3:68511/149‑1 (MQ=255)
                                                                                     ttCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttattta                                                                                     <  4:23718/124‑1 (MQ=255)
                                                                                       cATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAg                                                           <  7:70446/148‑1 (MQ=255)
                                                                                           gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGcaca                                                      >  3:108168/1‑149 (MQ=255)
                                                                                           gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGcaca                                                      >  7:288063/1‑149 (MQ=255)
                                                                                            tGATCTACACTTACAAGTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAt                                                     <  5:370467/149‑1 (MQ=255)
                                                                                              aaCTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                        >  7:64407/1‑112 (MQ=255)
                                                                                              aaCTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                        <  8:64407/112‑1 (MQ=255)
                                                                                              aaCTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTg                                                   <  7:197007/149‑1 (MQ=255)
                                                                                                   cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAcc                                               <  1:42517/148‑1 (MQ=255)
                                                                                                   cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                >  5:55925/1‑147 (MQ=255)
                                                                                                   cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                <  6:55925/147‑1 (MQ=255)
                                                                                                     cTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCa                                             >  8:9117/1‑148 (MQ=255)
                                                                                                      ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                        >  5:27113/1‑104 (MQ=255)
                                                                                                      ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                        <  6:27113/104‑1 (MQ=255)
                                                                                                      ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                           >  3:150877/1‑149 (MQ=255)
                                                                                                      ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                           >  6:355267/1‑149 (MQ=255)
                                                                                                      ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTATGTTTTTAATCGATTGCAGGCACATTGAAACCATa                                           >  8:133967/1‑149 (MQ=255)
                                                                                                             ttCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                                    >  7:196701/1‑117 (MQ=255)
                                                                                                             ttCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                                    <  8:196701/117‑1 (MQ=255)
                                                                                                                   gTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGt                                       <  7:10624/141‑1 (MQ=255)
                                                                                                                   gTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGt                                       >  8:10624/1‑141 (MQ=255)
                                                                                                                         cGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                          >  3:319989/1‑148 (MQ=255)
                                                                                                                                      gaagCCTCGCCTGATGCGAGGCTTTTTTGAGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTTTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGcc            <  5:95883/149‑1 (MQ=255)
                                                                                                                                       aagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagc                           <  4:263814/133‑1 (MQ=255)
                                                                                                                                       aagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagc                           >  3:263814/1‑133 (MQ=255)
                                                                                                                                           cTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGa       <  5:355267/149‑1 (MQ=255)
                                                                                                                                           cTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGa       <  7:246279/149‑1 (MQ=255)
                                                                                                                                              gCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGc     <  8:288063/148‑1 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAccca                             <  2:241698/122‑1 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAccca                             >  1:241698/1‑122 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                   <  2:303344/132‑1 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGc                   >  1:303344/1‑132 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCt   >  4:345514/1‑148 (MQ=255)
                                                                                                                                                cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTc  <  4:108168/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTC  >  NC_002947/289799‑290090

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: