Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,158,783 T→G I98M (ATT→ATG oxyR → oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,158,7830TG96.0% 259.5 / ‑3.3 76I98M (ATT→ATGoxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (35/38);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTG  >  NC_000913_3_bme_pgi/4158636‑4158921
                                                                                                                                                   |                                                                                                                                          
gCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGAGGATCAGGCACGACCAGGGATACGTGAGTTGAAGTTCCTAAATGAGATGACAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGcc                                                                                                                                          <  1:810132/150‑1 (MQ=255)
  ggACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCca                                                                                                                                        <  4:172966/150‑1 (MQ=255)
   gACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTGAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCcac                                                                                                                                       >  2:605601/1‑150 (MQ=255)
   gACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCcac                                                                                                                                       <  4:376339/150‑1 (MQ=255)
    aCCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCcaca                                                                                                                                      >  2:676954/1‑150 (MQ=255)
      cagcagTAAAGGGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGAAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGt                                                                                                                                    <  1:749364/150‑1 (MQ=255)
      cAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGt                                                                                                                                    >  1:404265/1‑150 (MQ=255)
       aGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGtt                                                                                                                                   <  4:406432/150‑1 (MQ=255)
             aaaGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCg                                                                                                                             >  3:265828/1‑150 (MQ=255)
                   ttgttCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTg                                                                                                                       <  3:187614/150‑1 (MQ=255)
                     gttCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCt                                                                                                                     >  1:264646/1‑150 (MQ=255)
                         aCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                >  2:415425/1‑135 (MQ=255)
                         aCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                <  1:415425/135‑1 (MQ=255)
                          cccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGc                                                                                                                <  3:365404/150‑1 (MQ=255)
                           ccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTg                                                                                                                                  <  2:638463/131‑1 (MQ=255)
                           ccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTg                                                                                                                                  >  1:638463/1‑131 (MQ=255)
                           ccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGt                                                                                                                            <  3:260235/137‑1 (MQ=255)
                           ccAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGt                                                                                                                            >  4:260235/1‑137 (MQ=255)
                            cAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCat                                                                                                              >  1:384213/1‑150 (MQ=255)
                            cAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGAAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCat                                                                                                              >  2:634202/1‑150 (MQ=255)
                                 gggAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCAtattat                                                                                                         >  2:823026/1‑150 (MQ=255)
                                          gctgGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtgct                                                                                                <  1:676954/150‑1 (MQ=255)
                                             ggtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCa                                                                                             <  2:404265/150‑1 (MQ=255)
                                              gtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGc                                                                                                                      <  2:684490/124‑1 (MQ=255)
                                              gtggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGc                                                                                                                      >  1:684490/1‑124 (MQ=255)
                                                ggATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCa                                                                                          >  2:117982/1‑150 (MQ=255)
                                                 gATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAg                                                                                         >  1:58074/1‑150 (MQ=255)
                                                  aTCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGa                                                                                        <  2:486869/150‑1 (MQ=255)
                                                     aGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCt                                                                                     >  2:476608/1‑150 (MQ=255)
                                                      ggCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCt                                                                                                                     <  3:348643/117‑1 (MQ=255)
                                                      ggCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCt                                                                                                                     >  4:348643/1‑117 (MQ=255)
                                                      ggCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCtt                                                                                    >  2:169888/1‑150 (MQ=255)
                                                       gcgcGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCttt                                                                                   >  1:6489/1‑150 (MQ=255)
                                                         gcgTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCa                                                                                                               <  1:271106/120‑1 (MQ=255)
                                                         gcgTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCa                                                                                                               >  2:271106/1‑120 (MQ=255)
                                                           gTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCaa                                                                               <  2:264646/150‑1 (MQ=255)
                                                                    tgcGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGaaa                                                                      <  1:634202/150‑1 (MQ=255)
                                                                       gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCt                                                                                                     >  3:77263/1‑116 (MQ=255)
                                                                       gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCt                                                                                                     <  4:77263/116‑1 (MQ=255)
                                                                       gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCt                                                                                     >  2:55053/1‑132 (MQ=255)
                                                                       gTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCt                                                                                     <  1:55053/132‑1 (MQ=255)
                                                                         gAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                >  2:849355/1‑87 (MQ=255)
                                                                         gAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                >  2:849290/1‑87 (MQ=255)
                                                                         gAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                <  1:849355/87‑1 (MQ=255)
                                                                         gAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGa                                                                                                                                <  1:849290/87‑1 (MQ=255)
                                                                             tGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGc                                                             >  2:594987/1‑150 (MQ=255)
                                                                              gAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCa                                                            >  4:490890/1‑150 (MQ=255)
                                                                              gAAAGTCCTTAAAGAGATGGCAAGCAAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACATTTCCAAAGCTGGAAATGTATCTGCa                                                            >  3:389884/1‑150 (MQ=255)
                                                                                   tCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAg                                                       <  1:823026/150‑1 (MQ=255)
                                                                                     cTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGca                                                     >  4:61817/1‑150 (MQ=255)
                                                                                         aaGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGa                                                 <  1:117982/150‑1 (MQ=255)
                                                                                               tGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                                                                       >  1:297530/1‑90 (MQ=255)
                                                                                               tGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                                                                       <  2:297530/90‑1 (MQ=255)
                                                                                               tGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacc                                           <  1:605601/150‑1 (MQ=255)
                                                                                                ggCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacca                                          <  1:169888/150‑1 (MQ=255)
                                                                                                ggCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccacca                                          <  3:61817/150‑1 (MQ=255)
                                                                                                      ccAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACt                                    <  2:58074/150‑1 (MQ=255)
                                                                                                       cagcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTg                                   <  1:594987/150‑1 (MQ=255)
                                                                                                         gcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCt                                                                                                     <  2:397038/82‑1 (MQ=255)
                                                                                                         gcagGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCt                                                                                                     >  1:397038/1‑82 (MQ=255)
                                                                                                                  gaCGATGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGa                        >  1:740486/1‑150 (MQ=255)
                                                                                                                       tGTCCGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCg                   <  2:6489/150‑1 (MQ=255)
                                                                                                                           cGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc                                                                                            <  1:714942/73‑1 (MQ=255)
                                                                                                                           cGGACCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCAc                                                                                            >  2:714942/1‑73 (MQ=255)
                                                                                                                              aCCGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACt            <  3:490890/150‑1 (MQ=255)
                                                                                                                                cGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                                                                       <  1:375363/57‑1 (MQ=255)
                                                                                                                                cGCTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                                                                       >  2:375363/1‑57 (MQ=255)
                                                                                                                                  cTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCa                                                                                             <  1:641375/65‑1 (MQ=255)
                                                                                                                                  cTGCACATTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCa                                                                                             >  2:641375/1‑65 (MQ=255)
                                                                                                                                     cacaTTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccac                                            <  2:313692/111‑1 (MQ=255)
                                                                                                                                     cacaTTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACccac                                            >  1:313692/1‑111 (MQ=255)
                                                                                                                                        aTTGGTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTg  >  3:394902/1‑150 (MQ=255)
                                                                                                                                            gTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAg                     <  4:153726/127‑1 (MQ=255)
                                                                                                                                            gTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAg                     >  3:153726/1‑127 (MQ=255)
                                                                                                                                            gTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCg    <  3:68920/144‑1 (MQ=255)
                                                                                                                                            gTTTGATGCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCg    >  4:68920/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTG  >  NC_000913_3_bme_pgi/4158636‑4158921

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