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‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCC  >  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‑‑gggg                                                                                                                                                 <  1:586874/151‑1 (MQ=255)
      aGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggg                                                                                                                                            >  2:504293/1‑151 (MQ=255)
         aCTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCTTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggggggg                                                                                                                                         >  2:338805/1‑150 (MQ=255)
         aCTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑ggggggggggg                                                                                                                                         >  2:123360/1‑150 (MQ=255)
           tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑gggggggggggct                                                                                                                                       >  2:334495/1‑148 (MQ=255)
              gCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAAt                                                                                                                                    <  1:1059534/151‑1 (MQ=255)
               cATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATc                                                                                                                                   <  2:213470/151‑1 (MQ=255)
                aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCg                                                                                                                                  <  1:728502/151‑1 (MQ=255)
                   tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGc                                                                                                                                <  2:608508/150‑1 (MQ=255)
                   tCGCACCTCATCATCTGCATGATATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGc                                                                                                                                <  2:433439/150‑1 (MQ=255)
                     gCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAgg                                                                                                                             >  2:483990/1‑151 (MQ=255)
                     gCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCCAATGGGCAggg                                                                                                                            >  2:43820/1‑151 (MQ=255)
                       aCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAggg                                                                                                                            >  1:552735/1‑150 (MQ=255)
                         cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCAAATCGGCAGGGa                                                                                                                           >  2:217766/1‑149 (MQ=255)
                         cTCATCATCTGCATGCAATATACTATAGAGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGAGGGGGGGATAAATGGCAGGGa                                                                                                                           >  1:184812/1‑149 (MQ=255)
                              catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGgccgc                                                                                                                    >  2:301161/1‑151 (MQ=255)
                                     atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                             >  1:647231/1‑150 (MQ=255)
                                         aaTATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                                         >  2:523458/1‑151 (MQ=255)
                                            ataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg                                                                                                     >  1:832346/1‑151 (MQ=255)
                                                  taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGTCTAATCGGCAGGGAAGGCCCCCCCGGATAgcgggcgg                                                                                                    >  2:56020/1‑147 (MQ=255)
                                                   aGGGGGGTTTTCTATATGTCAATGCATACCCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                                                    <  1:477560/146‑1 (MQ=255)
                                                    ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAg                                                                                              <  2:2874/151‑1 (MQ=255)
                                                     gggggTATTCTATATGTCAATGCAAACCCCCCTATAGTATATAACACGGGAGAATAAAATTTATCCGGGGAATGTGGGCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                                             <  1:301386/151‑1 (MQ=255)
                                                       gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGaa                                                                                           <  1:1059645/151‑1 (MQ=255)
                                                              cTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                                                    >  1:202150/1‑151 (MQ=255)
                                                                 tctGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                                 >  1:5357/3‑151 (MQ=255)
                                                                 tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAttt                                                                                  <  2:1041727/150‑1 (MQ=255)
                                                                  atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                <  1:767575/151‑1 (MQ=255)
                                                                     tCAATGCATCCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                                                                              <  2:263227/150‑1 (MQ=255)
                                                                       aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                           <  2:227931/151‑1 (MQ=255)
                                                                            aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACgg                                                                     >  1:859366/1‑151 (MQ=255)
                                                                               ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                                  <  1:741452/151‑1 (MQ=255)
                                                                                      ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                                              <  1:692164/149‑1 (MQ=255)
                                                                                      ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                            >  1:218041/1‑151 (MQ=255)
                                                                                       taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                             <  1:1006009/149‑1 (MQ=255)
                                                                                         gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                         <  1:898619/151‑1 (MQ=255)
                                                                                           atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                            >  1:600016/1‑147 (MQ=255)
                                                                                           atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                                       <  1:164422/151‑1 (MQ=255)
                                                                                             aaaaCACTGGAAAAAAAAATTTATCCGGTAAATGTGGTCGAAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                      <  1:408687/149‑1 (MQ=255)
                                                                                               aacaggggAAAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGTGCGGCATAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                      <  2:736022/142‑1 (MQ=255)
                                                                                                  ccTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGt                                                  >  2:443617/2‑149 (MQ=255)
                                                                                                  acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                                                 >  2:191830/1‑150 (MQ=255)
                                                                                                     ggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                              <  1:359030/150‑1 (MQ=255)
                                                                                                        gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          <  2:237915/151‑1 (MQ=255)
                                                                                                              aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                                    >  2:1060750/1‑151 (MQ=255)
                                                                                                              aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                                    >  1:720758/1‑151 (MQ=255)
                                                                                                                     ccGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCCTAATCGGAGGGGAAGGCCGCCCGGGTAGGGGGGGGGAAGAGGGATTAAGAATTTCCAGGCCAGAGGGGTTGCAAGTTGCGACCCGTTAAAATCATGGGTTGGGgt                               >  2:760535/1‑149 (MQ=255)
                                                                                                                       ggTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTc                            <  2:174535/151‑1 (MQ=255)
                                                                                                                          gAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa                         <  1:88407/150‑1 (MQ=255)
                                                                                                                          gAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                        >  2:636657/1‑151 (MQ=255)
                                                                                                                                   cGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACttt                >  1:206771/1‑150 (MQ=255)
                                                                                                                                   cGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa               <  1:1027783/151‑1 (MQ=255)
                                                                                                                                      aaaaCAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAcc             >  2:117843/1‑151 (MQ=255)
                                                                                                                                       aaaCAAAGAGGAAAGGGGGGGGGGGCTAATCGGAAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTg           >  1:650699/1‑151 (MQ=255)
                                                                                                                                        aaaaaaGAGGAAA‑‑GGGGGGTGGGGCTAACCGTCAGGGCAGGCCGCCCCGAATCGCGGGCGGCAGAAGGAATCAGAATTTCTAGGTCAGACGGGCCGCAAGTTGCAGACCGATAAAATCATCGCTTGGGGGGTCGTCCCACACTTTAc              >  2:1030221/4‑147 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          >  2:10421/1‑118 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                                          <  1:10421/118‑1 (MQ=255)
                                                                                                                                         aCAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGGGTCGTACCACACTTTACCTGcc         >  2:619951/1‑151 (MQ=255)
                                                                                                                                            aaGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc      >  1:743654/1‑151 (MQ=255)
                                                                                                                                            aaGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCa     <  2:345623/151‑1 (MQ=255)
                                                                                                                                             agagGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc      >  1:118725/1‑150 (MQ=255)
                                                                                                                                                ggAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGCCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGcc  >  2:974722/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCC  >  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: