Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 380,022 (G)10→11 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,0121.G77.0% 163.8 / 41.4 61intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (7/7);  new base G (31/16);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379874‑380135
                                                                                                                                              |                                                                                                                               
cTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                  >  2:161007/1‑139 (MQ=255)
                        tgcatgcaATATACTATGGGGGGGTATTCTATATGTCAATGAATACCCCCCTATAGTATATAACACTGGAGAAAAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAGG‑‑‑GGGGGGGGGG‑‑CTAATCGGCAggg                                                                                                      <  1:161007/139‑1 (MQ=255)
                            tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  >  1:673618/1‑139 (MQ=255)
                            tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGgc                                                                                                 >  1:1296477/1‑139 (MQ=255)
                              acagataCTATAGGGGGGTTTTCTATATGTCAAAGCATACCCCCCTATAGTATAAAACACTGGGGAATAAAATTTATCCGGTGAAGGTGGTCGGAAAACAAAGGGGAAGG‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGgcc                                                                                                <  1:26462/135‑1 (MQ=255)
                                 tataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGgccgcc                                                                                             >  1:1342099/1‑139 (MQ=255)
                                         aGGGGGGTATTTTATATGTCAAAGCATACCCCCCTATAGTAAATAACACTGGAGAATAAAAGTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAg                                                                                     <  2:1342099/139‑1 (MQ=255)
                                           acaggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACTCTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                   <  1:398805/136‑1 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CAAATCGGAAGGGAAGGCCGCCCCGGATAgctgt                                                                                 >  2:26462/1‑136 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggggg                                                                                                                        >  2:827891/1‑101 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggggg                                                                                                                        <  1:827891/101‑1 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggg                                                                                                                          >  2:31316/1‑99 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAAAAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggg                                                                                                                          <  1:31316/99‑1 (MQ=255)
                                                tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                              >  1:734154/1‑139 (MQ=255)
                                                          gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGCGGAAAG‑‑GGGGGGGGGG‑‑CGAATTGGCGGGGAAGGCGGCCCCGGATAGGGGGGGGCAGAAgctgt                                                                    >  2:398805/1‑135 (MQ=255)
                                                               tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAAGGGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                               >  2:1400157/1‑139 (MQ=255)
                                                               tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CAAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCAGAAGGACTAAGaa                                                               >  2:382619/1‑139 (MQ=255)
                                                               tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                              >  2:1421041/1‑139 (MQ=255)
                                                                    aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑gggggg                                                                                                                         <  2:105430/77‑1 (MQ=255)
                                                                    aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑gggggg                                                                                                                         >  1:105430/1‑77 (MQ=255)
                                                                     ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCa                                                         >  2:310211/1‑139 (MQ=255)
                                                                         ccTATAGTATATAACACTGGGGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                  <  2:945764/66‑1 (MQ=255)
                                                                         ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                  >  1:945764/1‑66 (MQ=255)
                                                                         ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGCGGGCGGCAGAAGGAATCAAAATTTCCAGGTc                                                     >  1:1027897/1‑139 (MQ=255)
                                                                                tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGG‑‑CTAATCGGCGGGGAAGGCCCCCCCGGATAGGGGGCGGCAAAAGGAAACAGAATTTCCAGGTAAGACg                                                >  2:787722/1‑139 (MQ=255)
                                                                                tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                                             >  1:192871/1‑139 (MQ=255)
                                                                                  tttaACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  <  1:1564855/83‑1 (MQ=255)
                                                                                  tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  >  2:1564855/1‑85 (MQ=255)
                                                                                  tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGTTAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  2:375825/1‑139 (MQ=255)
                                                                                  tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            <  1:310211/139‑1 (MQ=255)
                                                                                  tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  1:81921/1‑139 (MQ=255)
                                                                                  tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCt                                            >  2:960325/1‑139 (MQ=255)
                                                                                   ataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  2:1574222/139‑1 (MQ=255)
                                                                                    taACACTGGAGAATAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  2:81921/139‑1 (MQ=255)
                                                                                    taACACAGGAGAATAAAATTTATCCGGGGAATGTGGGCGGAAAATAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  2:673618/139‑1 (MQ=255)
                                                                                     aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CAAATGGGCGGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAGGGAATAAGAATTTCCAGGTCAGCCGGGCTGCa                                         >  1:1653936/1‑139 (MQ=255)
                                                                                       cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  <  2:1232134/80‑1 (MQ=255)
                                                                                       cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  >  1:1232134/1‑80 (MQ=255)
                                                                                       cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                       >  2:167073/1‑139 (MQ=255)
                                                                                          tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACgg                                               <  2:950343/128‑1 (MQ=255)
                                                                                          tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATAAGAATTTCCGGGTCAGACgg                                               >  1:950343/1‑128 (MQ=255)
                                                                                            gggaAAAAAATTTATCCGGGCAATGTGTTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  <  1:382619/137‑1 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGG‑‑CTAATGGGCAGGGAAGgccgcc                                                                                             >  2:1617006/1‑81 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGCATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  >  1:1090478/1‑139 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                >  2:1532982/1‑109 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  >  1:821723/1‑139 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGGGAATGTGGTGGGAAAACAAAGAGGAGAG‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                <  1:1532982/109‑1 (MQ=255)
                                                                                            gagaAAAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGG‑‑CTAATCGGCAGGGAAGgccgcc                                                                                             <  1:1617006/81‑1 (MQ=255)
                                                                                                 tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑gggggggggggcta                                                                                                                 >  1:995952/1‑52 (MQ=255)
                                                                                                 tAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAGAGGGAA‑‑‑gggggggggggcta                                                                                                                 <  2:995952/56‑5 (MQ=255)
                                                                                                  aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                           >  1:1249698/1‑139 (MQ=255)
                                                                                                         aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  <  2:388687/110‑1 (MQ=255)
                                                                                                         aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  >  1:388687/1‑110 (MQ=255)
                                                                                                         aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          >  2:1475752/1‑118 (MQ=255)
                                                                                                         aTCCGGTGAAAGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  1:1475752/118‑1 (MQ=255)
                                                                                                          tCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                   <  2:1296477/139‑1 (MQ=255)
                                                                                                               tGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGGGCCTAAACGGCAGGGAAGGCCGCCCCGGAAGCGGGCGGCAGAAGGAATCAAAATTTTCAGGTCCAGCGGGCTGCAAGTTGCCGACCGTTAGAataatcagc               >  2:233141/1‑136 (MQ=255)
                                                                                                                  aTGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg           <  2:1249698/139‑1 (MQ=255)
                                                                                                                   tgtgGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg           >  1:50219/1‑139 (MQ=255)
                                                                                                                      ggTCGGAAAACAAAGGGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  <  2:907948/97‑1 (MQ=255)
                                                                                                                      ggTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  >  1:907948/1‑97 (MQ=255)
                                                                                                                       ggcGGAAAACAAAGAGGGAAG‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  <  1:204386/46‑1 (MQ=37)
                                                                                                                       gTCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAgg                                                                                                  >  2:204386/1‑48 (MQ=255)
                                                                                                                        tCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                >  1:1213944/1‑65 (MQ=255)
                                                                                                                        tCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                <  2:1213944/65‑1 (MQ=255)
                                                                                                                        tCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc      >  1:549874/1‑139 (MQ=255)
                                                                                                                        tCGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc      >  2:1545121/1‑139 (MQ=255)
                                                                                                                        tCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGGCCAAATCGGCAGGGAGGGCCGCCCAGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGCCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc      >  2:1565061/1‑139 (MQ=255)
                                                                                                                         cGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                <  2:649735/80‑1 (MQ=255)
                                                                                                                         cGGAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                >  1:649735/1‑80 (MQ=255)
                                                                                                                           gAAAACAAAGAGGAAAG‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa   <  2:50219/139‑1 (MQ=255)
                                                                                                                           gAAAACAAAGAGGAAA‑‑‑GGGGGGGGGC‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATGGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCGGCCGGGCTGCAAGTTGCAGCCGGTTAAATTAATCGGTTGGGGTGCCGTAc  >  2:1476710/1‑139 (MQ=255)
                                                                                                                                 aaaGAGGAAAGGGGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  2:908486/95‑1 (MQ=255)
                                                                                                                                 aaaGAGGAAAGGGGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCCCCCCGGATAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          >  1:908486/1‑95 (MQ=255)
                                                                                                                                              |                                                                                                                               
CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379874‑380135

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