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.G83.6% 174.9 / 27.6 61intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (7/3);  new base G (34/17);  total (41/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

GGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379877‑380135
                                                                                                                                         |                                                                                                                              
ggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg                                                                                                                            >  2:739071/1‑139 (MQ=255)
                    aaaagacagatatACAAAAGGGGGGTATTCTATATGTCAATGCATGCGCCCCTATAATAAAGAATGGTGGAAAAGAAAATTTATCCGGTGAATGTGGAAGGAAAAAAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAgg                                                                                                       <  2:927111/130‑1 (MQ=255)
                       catgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCGGCTAATCGGCAGGGaa                                                                                                    >  2:175402/1‑139 (MQ=255)
                        atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAg                                                                                                   >  2:765269/1‑139 (MQ=255)
                                 aCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCgg                                                                                         >  1:214679/1‑139 (MQ=255)
                                      aGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa                                                                                      >  2:743475/1‑139 (MQ=255)
                                      aGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg                                                                                     >  2:330205/1‑139 (MQ=255)
                                        gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  <  1:5836/124‑1 (MQ=255)
                                        gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                  >  2:5836/1‑124 (MQ=255)
                                        gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                   >  1:276076/1‑139 (MQ=255)
                                        gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccctgttcttt                                                                                   >  1:927111/1‑131 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                 >  2:1014513/1‑139 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                 >  2:512181/1‑139 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                 >  2:550516/1‑139 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCAAATGGGAAGGGAAGGCCGCCCCGGATAggggg                                                                                 >  1:442609/1‑139 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggcta                                                                                                               <  2:881709/108‑5 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggcta                                                                                                               >  1:881709/1‑104 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggg                                                                                                                        <  2:563338/100‑1 (MQ=255)
                                          gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggg                                                                                                                        >  1:563338/1‑100 (MQ=255)
                                             tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                             >  2:921645/1‑139 (MQ=255)
                                               ttCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAACGGGGGGCa                                                                            >  1:1010819/1‑139 (MQ=255)
                                                    tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                       >  1:885935/1‑139 (MQ=255)
                                                       gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                                                    >  2:979924/1‑139 (MQ=255)
                                                          aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAACCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                                                 >  1:743864/1‑139 (MQ=255)
                                                           aTGCATACCCCCCTATAGAAAAAAACCCTGGGGAAAAAAATTTTTCCGGGGAAGGGGGTCGGAAAAAAAAGGGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                <  2:1010819/139‑1 (MQ=255)
                                                            tGCTTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                               >  1:717316/1‑139 (MQ=255)
                                                            tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                               >  2:444360/1‑139 (MQ=255)
                                                            tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCAAACGGGAAGGGAAGGCCCCCCCGGATGGGGGGCGGCAGAAGGAATCAGaa                                                               >  2:1033538/1‑139 (MQ=255)
                                                            tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggg                                                                                                                        >  1:1001535/1‑82 (MQ=255)
                                                            tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggg                                                                                                                        <  2:1001535/82‑1 (MQ=255)
                                                            tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                              >  1:568721/1‑139 (MQ=255)
                                                                tACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                           <  1:739071/139‑1 (MQ=255)
                                                                          taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                                                 <  2:628346/139‑1 (MQ=255)
                                                                              atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            <  2:214679/139‑1 (MQ=255)
                                                                               tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  1:873065/1‑139 (MQ=255)
                                                                               tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  1:76900/1‑139 (MQ=255)
                                                                               tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  1:343459/1‑139 (MQ=255)
                                                                               tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                            >  2:12370/1‑139 (MQ=255)
                                                                               tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCAGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATAAGAATTTCAGGGTAAGACGGCCt                                            >  2:787311/1‑139 (MQ=255)
                                                                                 aaaCACTGGAGAATAAAATTTAACCCGTGGATGGGGGTAAAAAACAAAAAGGGAAG‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                          <  1:787311/138‑1 (MQ=255)
                                                                                  aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                         >  2:818029/1‑139 (MQ=255)
                                                                                  aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                        >  2:106965/1‑139 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                        >  1:202488/1‑139 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                       >  2:356264/1‑139 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                       >  2:131567/1‑139 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggctaa                                                                                                              >  1:834553/1‑62 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggctaa                                                                                                              <  2:834553/67‑6 (MQ=255)
                                                                                    cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAGGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                      >  2:957161/1‑139 (MQ=255)
                                                                                     acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGCCGGGCTGCAAGt                                      >  1:773999/1‑139 (MQ=255)
                                                                                       tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                    >  2:312139/1‑139 (MQ=255)
                                                                                       tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggg                                                                                                                       >  2:149321/1‑56 (MQ=255)
                                                                                       tGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggg                                                                                                                       <  1:149321/56‑1 (MQ=255)
                                                                                         gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCGGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  >  2:365950/1‑139 (MQ=255)
                                                                                         gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  >  1:1026000/1‑139 (MQ=255)
                                                                                         gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                  >  1:198232/1‑139 (MQ=255)
                                                                                           gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGGGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                <  1:765269/139‑1 (MQ=255)
                                                                                           gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                <  1:550516/139‑1 (MQ=255)
                                                                                              tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAAt                                                                                                              >  2:872759/1‑58 (MQ=255)
                                                                                              tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAAt                                                                                                              <  1:872759/58‑1 (MQ=255)
                                                                                              tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                                                    >  2:1008209/1‑100 (MQ=255)
                                                                                              tAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                                                    <  1:1008209/100‑1 (MQ=255)
                                                                                                   tttATCCGGTGAATGTGGTCGGAAAACAACGAGGAAAGGGGGGGGGGGCTCATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAACCGGGCTGCAAGTTGCAGGCCGTTaaa                        >  1:1027922/1‑139 (MQ=255)
                                                                                                      aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  >  1:271820/1‑110 (MQ=255)
                                                                                                      aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                  <  2:271820/110‑1 (MQ=255)
                                                                                                      aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatc                     <  1:12370/139‑1 (MQ=255)
                                                                                                                tgtgGGCGGAAAACAAAGGGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg           <  1:1033538/139‑1 (MQ=255)
                                                                                                                 gtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg          <  1:131567/139‑1 (MQ=255)
                                                                                                                 gtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg          <  1:979924/139‑1 (MQ=255)
                                                                                                                   ggTCCGAAAAGATAGAGGGAG‑‑GGGGGTTGGTGCTAATCGGCAGGTGAGGCCACCACAGATATCCGGCGGCATAAGTAAACATAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCTTTAAAAGCATCGGTTGGGgtg        <  2:1027922/139‑1 (MQ=255)
                                                                                                                     tCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                          >  2:589130/1‑86 (MQ=255)
                                                                                                                     tCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                          <  1:589130/86‑1 (MQ=255)
                                                                                                                       ggAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg     <  1:743475/139‑1 (MQ=255)
                                                                                                                        gAAAACGAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCACGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc  >  2:648412/1‑139 (MQ=255)
                                                                                                                         aaaaCAAAGGGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc  <  1:356264/139‑1 (MQ=255)
                                                                                                                         aaaaCAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGCCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc  <  2:743864/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                              
GGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAC  >  NC_000913/379877‑380135

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