Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 198,073 G→A 5.8% R49H (CGC→CAC)  bamA → BamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi198,0730GA 5.8% 193.5 / 2.0 69R49H (CGC→CAC) bamABamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor
Reads supporting (aligned to +/- strand):  ref base G (27/38);  new base A (2/2);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

TGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCATTAC  >  NC_000913_3_bme_pgi/197929‑198220
                                                                                                                                                |                                                                                                                                                   
tGGCGATGAAAAAGTTGCCCATAGGGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGTTATTCATTTCGCAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGcaca                                                                                                                                                 <  2:224032/149‑1 (MQ=255)
gggCGATGAAAAAGTTGCTCATAGCGCCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGcaca                                                                                                                                                 <  1:525017/148‑1 (MQ=255)
 ggggATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGTTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCATCGTGTCGCCGTTTGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAg                                                                                                                                                <  2:73514/146‑1 (MQ=255)
 ggCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCTGCGCCTCCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAg                                                                                                                                                <  1:327881/149‑1 (MQ=255)
 ggCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATTCGGTGTTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAg                                                                                                                                                <  1:15679/149‑1 (MQ=255)
       gAAAAATTTGCGCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTCTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGaca                                                                                                                                          <  2:691171/149‑1 (MQ=255)
       gAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGaca                                                                                                                                          <  1:90968/149‑1 (MQ=255)
        aaaaaGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGacac                                                                                                                                         >  1:703740/1‑149 (MQ=255)
                 cTCATAGCGTTGCTGCTGTTTAGCAGCGCCACCGTATACGGTGGTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAatga                                                                                                                                <  1:144699/149‑1 (MQ=255)
                  tCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGCAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAatgat                                                                                                                               >  2:50505/1‑149 (MQ=255)
                  tCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCACACAGGCGACACGGTTAatgat                                                                                                                               >  1:123195/1‑149 (MQ=255)
                  tCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCACACAGGCGACACGGTTAatgat                                                                                                                               >  1:124804/1‑149 (MQ=255)
                     tAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGaa                                                                                                                            <  1:726779/149‑1 (MQ=255)
                     caGCGTCGGTGCTGTTTAGCAGCGCCACCGTATACGGTTCTGAAGGGTTCGTAGTGAATGATATTCATTTCGAAGGGCTTCAGAGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGaa                                                                                                                            <  1:344368/148‑1 (MQ=255)
                       gCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGGGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGa                                                                                                                          <  1:368551/149‑1 (MQ=255)
                         gTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGata                                                                                                                        <  1:653357/149‑1 (MQ=255)
                          tCGCTGCTGTTTAGCAGCGGCACCGGATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGatat                                                                                                                       <  2:790130/149‑1 (MQ=255)
                              tgctgTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGt                                                                                                                   >  1:138826/1‑149 (MQ=255)
                                ctgTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCTTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTaa                                                                                                                 <  1:775495/149‑1 (MQ=255)
                                ctgTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTaa                                                                                                                 <  1:739778/149‑1 (MQ=255)
                                ctgTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGACCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGATGATATCAGTaa                                                                                                                 >  1:95952/1‑149 (MQ=255)
                                 tgTTTAGCAGCGCCACCGTATACGGTGTTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGTCTTCAGCGTGTCGCCGTTGGTGCGGCTCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAAt                                                                                                                <  1:220905/149‑1 (MQ=255)
                                  gTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATa                                                                                                               <  2:882765/149‑1 (MQ=255)
                                  gTTTAGCAGCGCCACCGTATACGGTGCTGAAGGATTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTAGTGCGGCCCTCCTCAGTATGCCGTTGCGCACAGGCGACACGGTTAATTATGACAATATCAGTGATa                                                                                                               >  2:813067/1‑149 (MQ=255)
                                   tttAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATAc                                                                                                              <  2:249814/149‑1 (MQ=255)
                                   tttAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATAc                                                                                                              <  2:273491/149‑1 (MQ=255)
                                         acggcCACCGTATACGGTTCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGATGGCTTTCAGTGTGTCGTCGTTGGTGCGGTCCTTGTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTcg                                                                                                        <  1:310506/146‑1 (MQ=255)
                                                         ctGCAGAAGGGATCGTAGTCAAAGAAAATCATATCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTAcc                                                                                        >  2:757270/2‑149 (MQ=255)
                                                             tGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGca                                                                                                                                                   >  2:731339/1‑86 (MQ=255)
                                                             tGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGca                                                                                                                                                   <  1:731339/86‑1 (MQ=255)
                                                               aaGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAAc                                                                                  >  1:864385/1‑149 (MQ=255)
                                                                     tggtggTGAAAGATATTCATTTTCGAGGCCTTCAGCGGGTCGCTGTTGGTGCGGCTCTCCGCAGTATGCCGGTGCGCACAGGCGACATGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAg                                                                            <  2:62083/144‑1 (MQ=255)
                                                                       gTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGACGCTGATTTGGCCCTCCTCAGTATGCCGGTTCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTATGAGGa                                                                          >  2:148082/1‑149 (MQ=255)
                                                                             aaaGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTAc                                                                                         <  1:345945/128‑1 (MQ=255)
                                                                             aaaGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTAc                                                                                         >  2:345945/1‑128 (MQ=255)
                                                                             aaaGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCCGCACTCCTCAGTCTGCCGGAGCGCACAGGCGACACGGTTAACGATGAAGACATCAGTAATACCATTCGCGCTCTGCTTGCTACTGGCAACTTTGAGGATGTTcg                                                                    >  2:633911/1‑149 (MQ=255)
                                                                                 atcttCATTTCGAAGGCCTTCAGCGTGTCCGCGTTGGTGCGGACCTCCTCAGTATGCCGGTGCGCACAGGCGACACTGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTc                                                                <  1:774611/146‑1 (MQ=255)
                                                                                    ttCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTGCTCAGTATGCCGGAAAGAAAAAGAGGCGCTGTTTATAAAGAAGTGATACGTAACACCATAGGCGCTCGGTTTGATACCGGCAACTTTGAGGATGTTCGCGTCCtt                                                             >  2:223183/1‑149 (MQ=255)
                                                                                    ttCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCAtt                                                                                                          <  1:172464/104‑1 (MQ=255)
                                                                                    ttCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCAtt                                                                                                          >  2:172464/1‑104 (MQ=255)
                                                                                       aTTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACCTTGAGGATGTACGCGACCTTCGt                                                          >  1:872252/1‑149 (MQ=255)
                                                                                          tCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGAt                                                       <  2:740815/149‑1 (MQ=255)
                                                                                           cGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATg                                                      <  1:50505/149‑1 (MQ=255)
                                                                                           cGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTCGTCGATTTCAGCTTTGAGGATGTTCGCGTCCTTCGTGATg                                                      <  1:666554/149‑1 (MQ=255)
                                                                                           agAAGGCCTTCAGCGTGTCTCTGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGACGAAGATATCAGTAATACCATT‑‑CGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGt                                                    <  1:233139/148‑1 (MQ=255)
                                                                                             aaGGCCTGCAGCGTGTCGCCGTGGGTGCGCCCCTCCTCAGTATGCCGGTGCGCAAAGGCGACACGGTTAATGATGAAGATAGCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGt                                                    <  2:287394/149‑1 (MQ=255)
                                                                                               ggCCTTCAGCGTGTCGCCGTTGGTGCGGCGCTTCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGa                                                  <  1:787103/149‑1 (MQ=255)
                                                                                               ggCCTTCAGCGTGGCGCCGTTGGTGCGGCCCTCCTCAGTAGGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGa                                                  <  2:38586/149‑1 (MQ=255)
                                                                                                 ccTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATa                                                <  1:89534/149‑1 (MQ=255)
                                                                                                  cTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATAc                                               >  1:316633/1‑149 (MQ=255)
                                                                                                   ttCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATAcc                                              >  2:723347/1‑149 (MQ=255)
                                                                                                      aGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCtt                                           >  2:725771/1‑149 (MQ=255)
                                                                                                          tgtCCCCGGTGGTCCGGCCCTCCCCAGGATGCCGGTGCGCACAGGCGACAGGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTgg                                       <  2:134420/149‑1 (MQ=255)
                                                                                                                        cGGCCCTCCTCAGTATGCCGGTGCACACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGaa                         <  2:124804/149‑1 (MQ=255)
                                                                                                                        cGGCCCTCCTCAGTATGCCGGTGCACACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGaa                         <  2:123195/149‑1 (MQ=255)
                                                                                                                           cccTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGGAAGACCATGCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTAATACCCTTCTGGTTCAGGTAAAAGAAAGt                      >  1:172299/1‑149 (MQ=255)
                                                                                                                              tcctcAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCg                   >  1:175853/1‑149 (MQ=255)
                                                                                                                              tcctcAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCg                   >  1:5561/1‑149 (MQ=255)
                                                                                                                               cctcAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATATTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTCAAAGaagatccta                  >  1:65184/1‑142 (MQ=255)
                                                                                                                                ctcAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCACGGAAAAGAACGTCTGAc                 >  1:894552/1‑149 (MQ=255)
                                                                                                                                ctcAGTATGCCGGTCCGCACAAGCGACACGGTTAATGATGAAGCTATCAGTAATACCATTCGCGCTCTATTTGCTACCGGCAACTTAGAGGCTGTTTGAGTCCCTCCTGATCGCCATACTCATCTGACTAACGTAAAAGAACGTCCTAc                 >  2:797530/1‑149 (MQ=255)
                                                                                                                                 tcAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCAAACGCGCTCTGTTTGCTACCGGCAAATTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGAcc                >  1:253238/1‑149 (MQ=255)
                                                                                                                                    gTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCAtt             <  1:723347/149‑1 (MQ=255)
                                                                                                                                     tATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCAGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTg            >  2:368977/1‑149 (MQ=255)
                                                                                                                                       tGCAGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTAGCGCTCTGTTTGCTACCGGAAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCGTTCTGGATCAGGTAAAAGAACGTACAACCATTGcc          >  2:493897/1‑149 (MQ=255)
                                                                                                                                        gCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCa         <  2:175853/149‑1 (MQ=255)
                                                                                                                                           ggTGCGCACAGGCGACACGGTTAATGATGCAGATATCAGTAATACCATTCGCGTACTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCa      <  1:774081/149‑1 (MQ=255)
                                                                                                                                           ggTGCGCACAGGCGACACGGTTAAGTATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGTGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCa      <  2:509237/149‑1 (MQ=255)
                                                                                                                                             tGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCAtt    >  2:136062/1‑149 (MQ=255)
                                                                                                                                              gcgcACAGGCGACATGGTTAATGTTGAAGATATCAGTAATACCATTCGCTCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCATTa   <  2:200887/149‑1 (MQ=255)
                                                                                                                                               cgcACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCATTAc  >  1:101553/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGTTCGTAGTGAAAGATATTCATTTCGAAGGCCTTCAGCGTGTCGCCGTTGGTGCGGCCCTCCTCAGTATGCCGGTGCGCACAGGCGACACGGTTAATGATGAAGATATCAGTAATACCATTCGCGCTCTGTTTGCTACCGGCAACTTTGAGGATGTTCGCGTCCTTCGTGATGGTGATACCCTTCTGGTTCAGGTAAAAGAACGTCCGACCATTGCCAGCATTAC  >  NC_000913_3_bme_pgi/197929‑198220

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