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.G82.5% 168.3 / 34.9 57intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (8/2);  new base G (22/25);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 8.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCC  >  NC_000913/379864‑380157
                                                                                                                                                     |                                                                                                                                                     
ttcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gg                                                                                                                                                    >  1:302022/1‑151 (MQ=255)
 tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                                        <  2:799958/148‑1 (MQ=255)
 tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa                                                                                                                                                        <  1:235711/148‑1 (MQ=255)
  ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg‑g                                                                                                                                                 <  1:586874/151‑1 (MQ=255)
      aGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg‑ggg‑gg                                                                                                                                            >  2:504293/1‑151 (MQ=255)
         aCTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCTTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg‑ggg‑ggggg                                                                                                                                         >  2:338805/1‑150 (MQ=255)
         aCTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg‑ggg‑ggggg                                                                                                                                         >  2:123360/1‑150 (MQ=255)
           tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggg‑ggg‑gggggct                                                                                                                                       >  2:334495/1‑148 (MQ=255)
              gCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAAt                                                                                                                                    <  1:1059534/151‑1 (MQ=255)
               cATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATc                                                                                                                                   <  2:213470/151‑1 (MQ=255)
                aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCg                                                                                                                                  <  1:728502/151‑1 (MQ=255)
                   tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGc                                                                                                                                <  2:608508/150‑1 (MQ=255)
                   tCGCACCTCATCATCTGCATGATATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGc                                                                                                                                <  2:433439/150‑1 (MQ=255)
                     gCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAgg                                                                                                                             >  2:483990/1‑151 (MQ=255)
                     gCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CCAATGGGCAggg                                                                                                                            >  2:43820/1‑151 (MQ=255)
                       aCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAggg                                                                                                                            >  1:552735/1‑150 (MQ=255)
                         cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CAAATCGGCAGGGa                                                                                                                           >  2:217766/1‑149 (MQ=255)
                         cTCATCATCTGCATGCAATATACTATAGAGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGAGGG‑GGGG‑ATAAATGGCAGGGa                                                                                                                           >  1:184812/1‑149 (MQ=255)
                              catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGgccgc                                                                                                                    >  2:301161/1‑151 (MQ=255)
                                     atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                             >  1:647231/1‑150 (MQ=255)
                                         aaTATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                                         >  2:523458/1‑151 (MQ=255)
                                            ataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg                                                                                                     >  1:832346/1‑151 (MQ=255)
                                                  taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGGTCTAATCGGCAGGGAAGGCCCCCCCGGATAgcgggcgg                                                                                                    >  2:56020/1‑147 (MQ=255)
                                                   aGGGGGGTTTTCTATATGTCAATGCATACCCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                                                    <  1:477560/146‑1 (MQ=255)
                                                    ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAg                                                                                              <  2:2874/151‑1 (MQ=255)
                                                     gggggTATTCTATATGTCAATGCAAACCCCCCTATAGTATATAACACGGGAGAATAAAATTTATCCGGGGAATGTGGGCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                                             <  1:301386/151‑1 (MQ=255)
                                                       gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGaa                                                                                           <  1:1059645/151‑1 (MQ=255)
                                                              cTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                                                    >  1:202150/1‑151 (MQ=255)
                                                                 tctGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                                 >  1:5357/3‑151 (MQ=255)
                                                                 tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAttt                                                                                  <  2:1041727/150‑1 (MQ=255)
                                                                  atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                <  1:767575/151‑1 (MQ=255)
                                                                     tCAATGCATCCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                                                                              <  2:263227/150‑1 (MQ=255)
                                                                       aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                           <  2:227931/151‑1 (MQ=255)
                                                                            aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACgg                                                                     >  1:859366/1‑151 (MQ=255)
                                                                               ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                                  <  1:741452/151‑1 (MQ=255)
                                                                                      ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                                              <  1:692164/149‑1 (MQ=255)
                                                                                      ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                            >  1:218041/1‑151 (MQ=255)
                                                                                       taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                             <  1:1006009/149‑1 (MQ=255)
                                                                                         gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                         <  1:898619/151‑1 (MQ=255)
                                                                                           atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                            >  1:600016/1‑147 (MQ=255)
                                                                                           atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                                       <  1:164422/151‑1 (MQ=255)
                                                                                             aaaaCACTGGAAAAAAAAATTTATCCGGTAAATGTGGTCGAAAAACAAAGAGGAAAGGGGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                      <  1:408687/149‑1 (MQ=255)
                                                                                               aacaggggAAAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGTGCGGCATAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                      <  2:736022/142‑1 (MQ=255)
                                                                                                  ccTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGt                                                  >  2:443617/2‑149 (MQ=255)
                                                                                                  acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                                                 >  2:191830/1‑150 (MQ=255)
                                                                                                     ggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                              <  1:359030/150‑1 (MQ=255)
                                                                                                        gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          <  2:237915/151‑1 (MQ=255)
                                                                                                              aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                                    >  2:1060750/1‑151 (MQ=255)
                                                                                                              aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                                    >  1:720758/1‑151 (MQ=255)
                                                                                                                     ccGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGGCCTAATCGGAGGGGAAGGCCGCCCGGGTAGGGGGGGGGAAGAGGGATTAAGAATTTCCAGGCCAGAGGGGTTGCAAGTTGCGACCCGTTAAAATCATGGGTTGGGgt                               >  2:760535/1‑149 (MQ=255)
                                                                                                                       ggTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc                            <  2:174535/151‑1 (MQ=255)
                                                                                                                          gAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa                         <  1:88407/150‑1 (MQ=255)
                                                                                                                          gAATGTGGTCGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                        >  2:636657/1‑151 (MQ=255)
                                                                                                                                   cGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACttt                >  1:206771/1‑150 (MQ=255)
                                                                                                                                   cGGAAAACAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa               <  1:1027783/151‑1 (MQ=255)
                                                                                                                                      aaaaCAAAGAGGAAAGGGGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAcc             >  2:117843/1‑151 (MQ=255)
                                                                                                                                       aaaCAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGAAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTg           >  1:650699/1‑151 (MQ=255)
                                                                                                                                        aaaaaaGAGGAAA‑‑GGG‑GGGTGGGG‑CTAACCGTCAGGGCAGGCCGCCCCGAATCGCGGGCGGCAGAAGGAATCAGAATTTCTAGGTCAGACGGGCCGCAAGTTGCAGACCGATAAAATCATCGCTTGGGGGGTCGTCCCACACTTTAc              >  2:1030221/4‑147 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          >  2:10421/1‑118 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          <  1:10421/118‑1 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGGGTCGTACCACACTTTACCTGcc         >  2:619951/1‑151 (MQ=255)
                                                                                                                                            aaGAGGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc      >  1:743654/1‑151 (MQ=255)
                                                                                                                                            aaGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCa     <  2:345623/151‑1 (MQ=255)
                                                                                                                                             agagGAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc      >  1:118725/1‑150 (MQ=255)
                                                                                                                                                ggAAAG‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGCCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGcc  >  2:974722/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGG‑GGG‑GGGG‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCC  >  NC_000913/379864‑380157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑