Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 380,022 (G)10→9 intergenic (‑141/+47) frmR ← / ← yaiO regulator protein that represses frmRAB operon/outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0130G.97.1% 269.5 / 1.9 69intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (38/28);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379902‑380135
                                                                                                               |                                                                                                                          
tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgcc                                                                                                >  1:576229/1‑139 (MQ=255)
     tataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                                           >  2:1343754/1‑139 (MQ=255)
          tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                                                    >  1:532466/1‑125 (MQ=255)
          tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                                                    <  2:820149/125‑1 (MQ=255)
          tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                                                    >  1:820149/1‑125 (MQ=255)
          tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                                                    <  2:532466/125‑1 (MQ=255)
            taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  >  1:369707/1‑125 (MQ=255)
            taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  <  2:369707/125‑1 (MQ=255)
            taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAACTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                             <  1:959175/130‑1 (MQ=255)
            gaGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGAAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                             >  2:959175/2‑130 (MQ=255)
             aGGGGGGTCTTCTATATGTCAATGCATAACTCCCAATAGTATATAACACCGGAGAATAAAAAGTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGTCTGATGGGCAGGAAGGGCCGACCCGGATAgcg                                                                                   >  2:992010/1‑139 (MQ=255)
               gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                 >  1:370672/1‑139 (MQ=255)
                 gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg                                                                               >  1:854214/1‑139 (MQ=255)
                    tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                                                 <  2:571933/103‑4 (MQ=255)
                    tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                                                 >  1:571933/1‑100 (MQ=255)
                    tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                                                            >  2:1272679/1‑139 (MQ=255)
                              gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  >  2:632606/1‑107 (MQ=255)
                              gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  <  1:632606/107‑1 (MQ=255)
                                cAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                               >  1:1291785/1‑139 (MQ=255)
                                cAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAGAGGAAAGGGGGGGAGCTAATCGGCAGGGAATGCTGGCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                >  2:626348/1‑139 (MQ=255)
                                 aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                                                      >  2:683457/1‑85 (MQ=255)
                                 aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                                                      <  1:683457/85‑1 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCg                                                                                                            <  2:31267/92‑1 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCg                                                                                                            >  1:31267/1‑92 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                             >  2:736331/1‑139 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                             >  2:1147200/1‑139 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCGGGACAGACggccggc                                          >  2:172228/1‑134 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          >  2:1365191/1‑139 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          >  1:523371/1‑139 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          >  1:424784/1‑139 (MQ=255)
                                                         aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  >  2:1484624/1‑80 (MQ=255)
                                                         aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  <  1:1484624/80‑1 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                             <  1:827002/83‑1 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGTGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                     >  1:916649/1‑139 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGAAGGGAAGgccgcc                                                                                             >  2:827002/1‑83 (MQ=255)
                                                              tGGCGAATAAAATTTATCCTGCGAATTTGGGCGGGAAACTAAGGGGACAGGGCGTGGGCCAATCGGCTGGGAATGCCGCCCCGGATAGCGGGCGGCAGAAGGAATGAGTATTGCCAGGCCAGTCTGGCCACAAGATGct                                  >  2:1081611/1‑138 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  <  1:1365191/139‑1 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                 >  2:134061/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATGGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                >  2:746567/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAgg                                                                                                       <  1:228860/68‑1 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAgg                                                                                                       >  2:228860/1‑68 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                >  1:1158456/1‑107 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                <  2:1488768/107‑1 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                >  1:1488768/1‑107 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                <  2:1158456/107‑1 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCACATTGCaca                                >  2:131702/1‑137 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGAAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                >  2:706306/1‑139 (MQ=255)
                                                                     tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                  >  1:542623/1‑100 (MQ=255)
                                                                     tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                  <  2:542623/100‑1 (MQ=255)
                                                                      aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                             <  1:1176896/72‑1 (MQ=255)
                                                                      aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                             >  2:1176896/1‑72 (MQ=255)
                                                                      aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAGACCGTTa                          <  2:916649/139‑1 (MQ=255)
                                                                        aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                          >  1:685430/1‑89 (MQ=255)
                                                                        aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                          <  2:685430/89‑1 (MQ=255)
                                                                             aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                   <  2:424784/139‑1 (MQ=255)
                                                                               ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                                                 >  1:441668/1‑59 (MQ=255)
                                                                               ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                                                 <  2:441668/59‑1 (MQ=255)
                                                                                      aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg          <  1:746567/139‑1 (MQ=255)
                                                                                       tgtgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg        <  2:1291785/139‑1 (MQ=255)
                                                                                        gtgGGCGGAAAAAAAAAGGGAAAGGGGGGGGGCTAAATGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATTAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg        <  1:172228/139‑1 (MQ=255)
                                                                                         tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt       <  2:576229/139‑1 (MQ=255)
                                                                                           gTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                  >  1:224269/1‑126 (MQ=255)
                                                                                           gTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                  <  2:224269/126‑1 (MQ=255)
                                                                                            tCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  >  2:1118471/1‑45 (MQ=255)
                                                                                            tCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  <  1:1118471/45‑1 (MQ=255)
                                                                                            tCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                   <  1:130319/92‑1 (MQ=255)
                                                                                            tCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                   >  2:130319/1‑92 (MQ=255)
                                                                                             cGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa   <  1:706306/139‑1 (MQ=255)
                                                                                             cGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc  <  1:134061/139‑1 (MQ=255)
                                                                                               gAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                          >  1:1378440/1‑66 (MQ=255)
                                                                                               gAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                          <  2:1378440/66‑1 (MQ=255)
                                                                                                     aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa   >  2:506094/1‑132 (MQ=255)
                                                                                                     aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa   <  1:506094/132‑1 (MQ=255)
                                                                                                     aaaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                              <  2:241848/88‑1 (MQ=255)
                                                                                                     aaaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                              >  1:241848/1‑88 (MQ=255)
                                                                                                               |                                                                                                                          
TGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379902‑380135

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