Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,946 T→G 9.5% 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 9.5% 229.1 / 8.5 80intergenic (+30/‑136)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base G (4/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCAAA  >  NC_002947/289798‑290093
                                                                                                                                                     |                                                                                                                                                    
gCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGAt                                                                                                                                                      <  8:271255/149‑1 (MQ=255)
gCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGCCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGAt                                                                                                                                                      >  6:314494/1‑149 (MQ=255)
        aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGc                                                                                                                                               <  8:93697/148‑1 (MQ=255)
         aaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtt                                                                                                                                             >  2:32791/1‑149 (MQ=255)
             ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGACTGAGGCGAGGCtttttt                                                                                                                                         >  4:114055/1‑149 (MQ=255)
              cTGGGCTTCAGCATGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttctttt                                                                                                                                        <  6:233137/149‑2 (MQ=255)
               tGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCTCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttc                                                                                                                                            <  1:60351/144‑2 (MQ=255)
               tGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgt                                                                                                                                       >  8:303287/1‑149 (MQ=255)
               tGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttc                                                                                                                                            >  2:60351/1‑143 (MQ=255)
                  gCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAATGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                                     <  1:209458/149‑1 (MQ=255)
                    ttCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTg                                                                                                                                  >  3:304151/1‑149 (MQ=255)
                       aGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgtg                                                                                                                                      >  4:269490/1‑142 (MQ=255)
                       aGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTtgtg                                                                                                                                      <  3:269490/142‑1 (MQ=255)
                       aGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAAAAGCCTCGCATCAGGCGAGGCttcttt                                                                                                                                         <  4:261009/139‑1 (MQ=255)
                       aGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAAAAGCCTCGCATCAGGCGAGGCttcttt                                                                                                                                         >  3:261009/1‑139 (MQ=255)
                          cTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                           <  5:51899/134‑1 (MQ=255)
                          cTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                           >  6:51899/1‑134 (MQ=255)
                             cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGAttt                                                                                                                           <  4:304151/147‑1 (MQ=255)
                              gTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGCACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtata                                                                                                                        <  3:61108/149‑1 (MQ=255)
                               tCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTatat                                                                                                                       <  5:212046/149‑1 (MQ=255)
                                 aaGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGCCTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTCATGCGAGGCTTTTTTGTGCCTGCAGATTTATATg                                                                                                                      >  1:274684/1‑148 (MQ=255)
                                    gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtt                                                                                                                  >  1:2120/1‑149 (MQ=255)
                                         tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                             >  3:309852/1‑149 (MQ=255)
                                          aCCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGtt                                                                                                            >  4:347941/1‑149 (MQ=255)
                                           ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                             <  2:2120/147‑1 (MQ=255)
                                           ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTta                                                                                                           >  8:158342/1‑149 (MQ=255)
                                              tGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                                         >  8:276014/1‑132 (MQ=255)
                                              tGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                                         <  7:276014/132‑1 (MQ=255)
                                                  aGCGGCGAAAGCGCCCTGGACACTGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGt                                                                                                                     >  7:55252/1‑132 (MQ=255)
                                                  aGCGGCGAAAGCGCCCTGGACAAAGACCAGACTAGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGt                                                                                                                     <  8:55252/132‑1 (MQ=255)
                                                   gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTAtaaataaa                                                                                                   <  6:320552/149‑1 (MQ=255)
                                                     ggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGAGAGGCTTTTTTGTGACTGCAGATTTATATGTATTTTGGTTATaaataaata                                                                                                 >  4:20158/1‑148 (MQ=255)
                                                        gAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtt                                                                                              <  8:119453/149‑1 (MQ=255)
                                                          aaGCTCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACTCTTACAACTTCCTGTAATCCGTATCTTAATAAGAAGCCTCGCCTGATTGGAGGCTTTTTTGTGCGTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                            <  5:93670/149‑1 (MQ=255)
                                                               cccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttatttat                                                                                        >  3:9676/1‑148 (MQ=255)
                                                                 cTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCt                                                                                     >  2:203191/1‑149 (MQ=255)
                                                                    gagaaaGACCAGACTCGCTTCATCGTGAACTACAATTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGACTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAATTTGTTATTTATTCttt                                                                                   <  2:274684/145‑1 (MQ=255)
                                                                     aCAAAGACCAGACTCGCTTCATCGTGAATTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttttt                                                                                 <  2:39838/149‑1 (MQ=255)
                                                                       aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                               >  3:36283/1‑149 (MQ=255)
                                                                       aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                               >  2:308101/1‑149 (MQ=255)
                                                                          gACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                            <  7:303287/149‑1 (MQ=255)
                                                                             cAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTaa                                                                         <  4:36283/149‑1 (MQ=255)
                                                                               gACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                      >  6:91828/1‑149 (MQ=255)
                                                                                  tCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                                   >  7:161010/1‑149 (MQ=255)
                                                                                     cTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTg                                                                 >  3:351996/1‑148 (MQ=255)
                                                                                     cTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGc                                                                <  5:71295/149‑1 (MQ=255)
                                                                                      ttCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                                   >  3:135890/1‑145 (MQ=255)
                                                                                      ttCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                                   <  4:135890/145‑1 (MQ=255)
                                                                                        cATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAg                                                              <  7:158342/148‑1 (MQ=255)
                                                                                                 cTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                      <  4:52842/131‑1 (MQ=255)
                                                                                                 cTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                                      >  3:52842/1‑131 (MQ=255)
                                                                                                 cTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGACTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTAGTCTTTTTGTTTTAAATCGATTGCAGGCAAATTGaa                                                    >  8:69556/1‑149 (MQ=255)
                                                                                                  tACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                   <  1:308101/149‑1 (MQ=255)
                                                                                                  cgatcTTACAACTTCCTTTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                   <  7:224212/145‑1 (MQ=255)
                                                                                                    cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAccc                                                 <  8:161010/149‑1 (MQ=255)
                                                                                                       ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                              >  1:53204/1‑149 (MQ=255)
                                                                                                       ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGACTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTATTTTTGTTTTTAATCGATTGCAGGAACATTGAACCAATa                                              >  3:296200/1‑149 (MQ=255)
                                                                                                         aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                   >  8:224212/1‑142 (MQ=255)
                                                                                                         aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACg                                            >  8:97529/1‑149 (MQ=255)
                                                                                                          cAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                   <  6:229779/141‑1 (MQ=255)
                                                                                                          cAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAc                                                   >  5:229779/1‑141 (MQ=255)
                                                                                                          cAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACgg                                           <  4:309852/149‑1 (MQ=255)
                                                                                                           aaCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACgg                                           >  1:357483/1‑148 (MQ=255)
                                                                                                            aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                         <  2:357483/149‑1 (MQ=255)
                                                                                                            aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                         >  4:24320/1‑149 (MQ=255)
                                                                                                            aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                         <  5:314494/149‑1 (MQ=255)
                                                                                                            aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                         <  4:351996/149‑1 (MQ=255)
                                                                                                            aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTAAATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                         <  3:347941/149‑1 (MQ=255)
                                                                                                             cTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATTCGATTGCAGGCACATTGAACCCATACGGTa                                         <  1:32791/149‑1 (MQ=255)
                                                                                                                    ttaggcGTAGCTTAAAAATAAGCCTCGCCTGATGCGAGGATTTTTTTTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                                  <  3:20158/143‑1 (MQ=255)
                                                                                                                    tAAGCTGTAGGTGAATATGAAGCCTCGCCTTATGCTAGGCTTTTTTGTTCATGCAGATTTATTTGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTATTCGATTGAAGGCACATTGAACCCATACGGTAAAGAAcc                                  <  4:296200/148‑1 (MQ=255)
                                                                                                                     aaGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATGTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                                  >  6:364469/1‑147 (MQ=255)
                                                                                                                      aGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccag                               <  3:24320/149‑1 (MQ=255)
                                                                                                                        ccGTAGCTGAAAAAGAACCCTCGCCTGCTGCGAGGCTTTTTTGTCCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTAAACTATTTGTTTTTAAACGATTGCAGGCACATTGAACCCAGACGGTACAGACcccagcc                             >  2:195127/1‑149 (MQ=255)
                                                                                                                                gAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGcc                     >  4:85790/1‑149 (MQ=255)
                                                                                                                                 aaaaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                                  <  7:27706/103‑1 (MQ=255)
                                                                                                                                 aaaaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                                  >  8:27706/1‑103 (MQ=255)
                                                                                                                                   aaaaaaGCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGGTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGACCCCATACGGGACAGCACCCAGCCCAGGAGCCg                    >  5:157919/1‑147 (MQ=255)
                                                                                                                                    aagaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAg                 <  2:53204/149‑1 (MQ=255)
                                                                                                                                            cTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACACAACCCAGCCCAGGAGCCGCAGCCCATGa          <  1:203191/149‑1 (MQ=255)
                                                                                                                                                 cTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTc     <  4:9676/149‑1 (MQ=255)
                                                                                                                                                   gATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCa    >  3:280907/1‑148 (MQ=255)
                                                                                                                                                     tGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCaaa  <  5:91828/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATGAGCCTCAAA  >  NC_002947/289798‑290093

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: