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.G87.7% 172.6 / 20.7 57intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base G (33/17);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC  >  NC_000913/379865‑380153
                                                                                                                                                     |                                                                                                                                                
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                                   <  1:47782/148‑1 (MQ=255)
    gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggg                                                                                                                                         >  1:334786/1‑151 (MQ=255)
      gTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAATACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggggg                                                                                                                                       >  1:173689/1‑151 (MQ=255)
        aCTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGCGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggg                                                                                                                                      >  1:350006/1‑150 (MQ=255)
         cTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggc                                                                                                                                    >  1:9777/1‑149 (MQ=255)
         cTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggggggg                                                                                                                                     <  1:404896/150‑2 (MQ=255)
            ggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACGAAGAGGAAA‑‑GGGGGGGGGCTAATc                                                                                                                               <  2:373127/151‑1 (MQ=255)
               aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCg                                                                                                                              >  2:382939/1‑151 (MQ=255)
                     ccccTCATCATCTGCATGCAATATACTATAGGGGGGTATTTTATATGTCAATGCATACCCCCCTATAGTATAAAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAAAGGAAAGGGGGGGGGGGCTAATCGGCAgg                                                                                                                         <  2:267308/148‑1 (MQ=255)
                      aCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAggg                                                                                                                        >  1:705484/1‑150 (MQ=255)
                       ccTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGa                                                                                                                       >  1:483754/1‑150 (MQ=255)
                                cTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                                                             >  1:458675/1‑151 (MQ=255)
                                  gcatgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                                                               <  2:705435/147‑1 (MQ=255)
                                     tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                         >  1:286014/1‑149 (MQ=255)
                                      gcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                         >  1:412081/1‑149 (MQ=255)
                                       caATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc                                                                                                      >  1:347750/1‑151 (MQ=255)
                                         atatatTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                                    <  1:831403/151‑1 (MQ=255)
                                           ataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                   >  1:402936/1‑151 (MQ=255)
                                            taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg                                                                                                 >  1:451361/1‑151 (MQ=255)
                                            taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg                                                                                                 >  2:685971/1‑151 (MQ=255)
                                               tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                               >  2:866275/1‑151 (MQ=255)
                                                   ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATGGGAAGGGAAGGCCCCCCCGGATAGCGGGCGGAAGa                                                                                            >  1:131611/1‑150 (MQ=255)
                                                    gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                                                        <  2:285988/151‑1 (MQ=255)
                                                           ttCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                                                                     >  2:308022/1‑148 (MQ=255)
                                                                   gTCAATGCATCCCCCCCTATAGTATATAACACTGGAAAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                                                                          <  1:336024/151‑1 (MQ=255)
                                                                   gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                             >  1:84890/1‑148 (MQ=255)
                                                                      aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                       >  1:477653/1‑151 (MQ=255)
                                                                        tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                       >  2:149606/1‑149 (MQ=255)
                                                                            tACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACgg                                                                 >  2:689015/1‑151 (MQ=255)
                                                                             aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                                >  1:72859/1‑151 (MQ=255)
                                                                               cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGAAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                                >  1:540331/1‑149 (MQ=255)
                                                                                ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                             <  2:281707/151‑1 (MQ=255)
                                                                                ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                             <  2:597378/151‑1 (MQ=255)
                                                                                  ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                            >  1:548674/1‑150 (MQ=255)
                                                                                    tataGTATATAACCCTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGAGGGGGGCTAATCGGCAGGGAAGGCCGCCCTGGATAGCGGGCGGCAGAAGGAATAAGAATTTCCAGGACCGACGGGCTGc                                                            >  2:810932/1‑148 (MQ=255)
                                                                                    tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                            >  1:41664/1‑116 (MQ=255)
                                                                                    taaaGTATATAACACTGGAAAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                            <  2:41642/113‑1 (MQ=255)
                                                                                     ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCGGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                           >  2:261015/1‑148 (MQ=255)
                                                                                      taGTATATAACACTGGAGAATAAAATTTATCCGGTTAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGGTAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                         <  2:543216/149‑1 (MQ=255)
                                                                                          atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                     >  2:311286/1‑149 (MQ=255)
                                                                                          atataACACTGGAGAATAAAAGTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                                   <  2:365042/151‑1 (MQ=255)
                                                                                                    ggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                          <  1:79377/150‑1 (MQ=255)
                                                                                                           aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCGGGTCGGACGGGCTGAAGGTTGCAGGCCGTTAAAATCATCgg                                  >  1:483458/1‑151 (MQ=255)
                                                                                                            aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGt                                 <  1:103951/151‑1 (MQ=255)
                                                                                                            aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGAAGGGAAGGCCGCCCCGGAAAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGt                                 >  2:787044/1‑151 (MQ=255)
                                                                                                              aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAGTCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                               <  1:224993/151‑1 (MQ=255)
                                                                                                                ttATCCGGTGAATGTGGACGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                              <  1:347412/150‑1 (MQ=255)
                                                                                                                 tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGTCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGTCAGAAGGAATAAGAATTTCCAGGACAGACGGGCTGCAAGTTGCAGACCTTTAAAATCATCGGTTgg                              >  2:587296/1‑149 (MQ=255)
                                                                                                                   tCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg                          >  1:471506/1‑151 (MQ=255)
                                                                                                                    ccGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc                        >  1:588310/1‑151 (MQ=255)
                                                                                                                     cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc                        >  1:297054/1‑151 (MQ=255)
                                                                                                                     cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc                        >  1:174175/1‑151 (MQ=255)
                                                                                                                      ggTGAATGTGGTCGGAAATCAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGTAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg                       <  2:460039/151‑1 (MQ=255)
                                                                                                                      ggTGAATGTGGCCGGAAAACAAAGAGGAAAGGGGGGGGGGGGTAATCGGCGGGGAAGGGCGCCCCGGATAACGGGGTGGACAAGGAATAAGATTTTCAGGGTCAGACGGGCTGAAAGTGGCAGACAGAAAAACTCATCGGTTggggaggcg                       >  2:585637/1‑146 (MQ=255)
                                                                                                                           aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACca                  <  1:452841/151‑1 (MQ=255)
                                                                                                                               gggcGGAAAAAAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACt              <  1:689064/148‑1 (MQ=255)
                                                                                                                                  cGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACttt            >  2:287529/1‑150 (MQ=255)
                                                                                                                                   ggAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGACGCCCGGGAAAGCGGGCGGTAGAAGGAATCAGAATTTCGAGGTCAGAGGGGCTGCAAGTTGCAGCCCGTTAAAATCATGGGTTGGGGTGTGGTACCACACt              >  1:461913/1‑147 (MQ=255)
                                                                                                                                          aaaGAGGAAGG‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGAAGAAGGAATAAGAATTCCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTGGGGGTGTCGTACCACACTTTACCTGCCGt   <  2:458628/151‑1 (MQ=255)
                                                                                                                                          aaaGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg    <  1:561569/150‑1 (MQ=255)
                                                                                                                                          aaaGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt   >  1:676653/1‑151 (MQ=255)
                                                                                                                                           aaGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc  >  1:485/1‑151 (MQ=255)
                                                                                                                                             gagGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc  >  1:666398/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC  >  NC_000913/379865‑380153

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