Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,923 A→G 18.1% intergenic (+87/‑99) PP_1214 → / → ruvC transcriptional regulator/crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,392,9230AG18.1% 251.4 / 40.5 113intergenic (+87/‑99)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (37/50);  new base G (13/6);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 114 total aligned reads displayed.

GACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTA  >  NC_002947/1392774‑1393068
                                                                                                                                                     |                                                                                                                                                    
gacgTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGGAAACCCGCTCCGGCAGGTATTa                                                                                                                                                      >  2:3085708/1‑150 (MQ=255)
 acgTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAc                                                                                                                                                     <  2:4302806/150‑1 (MQ=255)
      cAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAGACCCGCGCCTACAGGTCTTAAGcccc                                                                                                                                               >  1:5080852/1‑149 (MQ=255)
         aaCGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCAttt                                                                                                                                            <  2:2290043/150‑1 (MQ=255)
          aCGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTc                                                                                                                                          <  2:1637901/150‑1 (MQ=255)
            gTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaa                                                                                                                                        >  1:5094534/1‑150 (MQ=255)
             tCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaaa                                                                                                                                       <  2:1271927/150‑1 (MQ=255)
             tCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaaa                                                                                                                                       <  2:4152926/150‑1 (MQ=255)
              cTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAt                                                                                                                                      <  1:1749282/150‑1 (MQ=255)
               tATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAtg                                                                                                                                     <  1:92568/150‑1 (MQ=255)
                 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCGCCGCCACGTGTTTCGCCCAGTTCAAAtggg                                                                                                                                   >  1:1220796/1‑148 (MQ=255)
                 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAtgtg                                                                                                                                   <  2:3676732/150‑1 (MQ=255)
                 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAtgtg                                                                                                                                   <  1:3812175/150‑1 (MQ=255)
                 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtgg                                                                                                                                    <  1:1488707/150‑2 (MQ=37)
                 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtgg                                                                                                                                    <  2:4837828/150‑2 (MQ=37)
                   ccAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggt                                                                                                                                 <  1:4876873/150‑1 (MQ=255)
                   ccAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtgggc                                                                                                                                  <  1:1957576/150‑4 (MQ=37)
                   ccAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtgggc                                                                                                                                  <  1:3706819/150‑4 (MQ=25)
                   ccAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtgggc                                                                                                                                  <  1:4860019/150‑4 (MQ=25)
                    cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAAACCCCGCCAACAGGTAGTAAGCCCAGTTCACATggggtg                                                                                                                                >  1:2104096/1‑150 (MQ=255)
                    cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggtg                                                                                                                                <  1:1032157/150‑1 (MQ=255)
                    cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCGTTCGCGGGGTGAACCGCTCCGACAGGTTTTTCGCCCAGTTGAAATgtggtg                                                                                                                                >  1:2104052/1‑150 (MQ=255)
                     aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCGCCAACAGGTATTACGCCCAAATAAAATGTGGTGc                                                                                                                               >  2:785955/1‑150 (MQ=39)
                     aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTAGCGGGTTAAGCCGCGCCAAGAGGTATTTAGCCCATTTGAAATgtggtgg                                                                                                                               >  2:482099/1‑149 (MQ=255)
                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                          <  2:1875345/150‑1 (MQ=255)
                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                          <  1:55804/150‑1 (MQ=255)
                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                          <  1:1186005/150‑1 (MQ=255)
                            aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCt                                                                                                                        <  2:2719953/150‑1 (MQ=255)
                              aTTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGt                                                                                                                      <  2:4194381/150‑1 (MQ=255)
                               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTa                                                                                                                     <  2:4524026/150‑1 (MQ=255)
                               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTa                                                                                                                     <  1:4928497/150‑1 (MQ=255)
                                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                              <  2:4752715/150‑1 (MQ=255)
                                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                              <  1:3834460/150‑1 (MQ=255)
                                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                              <  2:4752747/150‑1 (MQ=255)
                                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                              <  2:3867287/150‑1 (MQ=255)
                                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                              <  1:124856/150‑1 (MQ=255)
                                        agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGt                                                                                                            <  1:4141786/150‑1 (MQ=255)
                                        agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGCGCAGTACCTGTAGGAGCGGGt                                                                                                            <  1:4772719/150‑1 (MQ=255)
                                          aTCATTGAAAACCTCCGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCTATTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGCTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGatt                                                                                                          <  2:369646/150‑1 (MQ=38)
                                          aTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGttt                                                                                                          <  1:953420/150‑1 (MQ=255)
                                          aTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGttt                                                                                                          <  1:667113/150‑1 (MQ=255)
                                           tCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCGCCAACAGGTATTAAGCCCAGTTCAAATGTGGTGCAGTACCGGTAGGGGCGGGGTTa                                                                                                         >  2:3225359/1‑150 (MQ=38)
                                            cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCTACTAAAGGTATTAAGCCCAATTAAAAAGTTGTGCAGTAACAGTTGGAGCGGGGTTAc                                                                                                       >  2:1570490/1‑150 (MQ=255)
                                            cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAc                                                                                                       <  2:4220597/150‑1 (MQ=255)
                                            cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGATTAc                                                                                                       <  1:3245966/150‑1 (MQ=255)
                                               ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAccc                                                                                                     >  2:315923/1‑149 (MQ=255)
                                               ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAccc                                                                                                     <  1:315923/149‑1 (MQ=255)
                                                gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                   <  2:4425017/150‑1 (MQ=255)
                                                       tCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCACTACCAGTATGAGCCGGTTTAACcgcaccgaca                                                                                            >  1:3790326/1‑143 (MQ=38)
                                                         ggCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCTTTTGATAGCGCCTGGGCGGGCCCTTTCGCGGGTAAACCCGCTCCTCCAGGTATTCCGCCCATTTCAAATGTGGTGCAGTACCTGTAGGGGCGGGTTTACCCGCGAAAGGGcc                                                                                          <  2:1147061/150‑1 (MQ=37)
                                                                       ggCTTTTGGGGCCACCTGTGTGCGTTGTAGCGCGCGGGTGCCGCCCGTTCCGCGGGTAACCCCGCTCCACCCGTTATGCGGCCCTGTTCAACGGTGGTGCCGTCCGTGTTGGCGCGGTTTTCCCCGCGAATGGGCCGGCGCAGTTTGTag                                                                            <  2:3674697/150‑1 (MQ=255)
                                                                                  ccAACTCTGAGCGTTTTATTGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGTATTACGCCCATTTCAAATGTGGTGCTGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGcatcat                                                                 <  2:5078988/150‑1 (MQ=255)
                                                                                            gCGTTTGATAGCGCCTGTGCCGGCCCTTTCCCGGGTAAAACCCCACCAAAACGGCAAAAGCCCAATTAAAAATTGGTGCACTAACTGGAGGAGCGGGGTTAACCCCGAAACGGCCCGCGCCCTTTTTGGAGAAGCAAAAATTTTTAGACg                                                       >  2:3981021/1‑150 (MQ=255)
                                                                                                       cgcCTGTGCCGGCCCTTTGGGGGTTAAACCCGCTCCTACAGTTATTACGCCCTTTTCAAATGTGGTGCAGTACCTTTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGAc                                            <  1:4653715/150‑1 (MQ=38)
                                                                                                             tgCCGGCCCTTTCGCGGGTAAAACCGCTCCGACACGTATTACGCCCAATTCAAATGGGGGGCAGGACCGTTAGGAGCGGGGTTAACCGCGAAAGGGCCGGCGCCGGTTGTATAGTAGCAACCTTTTTTAGACGGGCTTTATTAcgcagag                                      >  2:3364191/1‑144 (MQ=255)
                                                                                                             tgCCGGCCCTTTCGCGGGTAAAACCGCTCCAAAAGGTATTACGCCCAGTTCAAAAGTTGTTCATTACCGGTATGAGCGGGTTTAACCGCGAAAGGGCCGGCGCAGTTTTTAGAGATGGATCAATTTTTCGTCGGGCGTTATGACTGTGAt                                      >  2:3364205/1‑150 (MQ=255)
                                                                                                             tgCCGGCCCTTTAGCGGGTAAACCCCCTCCTAAAGGTATTAAGCCCAATTCAAATGTGGTGCACTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGGTTGTAGAGAGGCATCAATTTTTAGAAGGGCGTTATGACTCTGAt                                      >  2:3364236/1‑150 (MQ=25)
                                                                                                              gCCGGCCCTTTCGGGGGTAAACCCGCTCCCCCCGGTATTCCGCCCTTTACAATTGTGGTGCTGTACCTGTAGGCGGGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGAtt                                     <  2:1220796/150‑1 (MQ=37)
                                                                                                                 gccccTTGCGCGTGTAAACCCGCTCCTACCGGAATACCCCCATTTACACATGTGGTGCTTTCCCGTTAGTTGGGGTTTACCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                  <  2:2599794/148‑1 (MQ=255)
                                                                                                                      tttCGCGGGTAAAACCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTTCAGTACCTGTATGATCGGGTTTACCACCGAAAGGGCCGGCGCAGTTTGTAGAGATGAATCATTTTTTCGAAGGGCTTTATTACTCTGATTCTAgtgat                             >  2:3979897/1‑146 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                       >  1:172572/1‑150 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                       >  1:172278/1‑150 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                       >  1:4672718/1‑150 (MQ=17)
                                                                                                                            gggTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                       >  1:4672724/1‑150 (MQ=17)
                                                                                                                            gggTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                       >  1:4672725/1‑150 (MQ=17)
                                                                                                                              gTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCg                     >  1:2475590/1‑150 (MQ=17)
                                                                                                                                aaaCCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGc                   >  1:127731/1‑150 (MQ=255)
                                                                                                                                aaaCCCGCGCCTCCCGTTATTACGCCCATTTCAACGGGGGTGATGTACCTGTTGGAGGGGGTTTACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGc                   <  2:3416212/150‑1 (MQ=255)
                                                                                                                                aaaCCCGCCCCTCCCGTTTTTACGCCATTTTCAATTTGGGTGCAGTACCTGTTGGTGGGGTTTTCCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGc                   <  2:844329/150‑1 (MQ=255)
                                                                                                                                 acccggCTCCTACCGTTTTCCCGCCATTTTCAATTGTGTTGCTGTCCGTTTTGTAGGGGTTTAACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCt                  <  2:254452/145‑1 (MQ=255)
                                                                                                                                  aCCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                 >  1:3828960/1‑150 (MQ=18)
                                                                                                                                   cccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCCGTGTGTAGAGAGGCCTCATTTTTTCGACGGGCTTTCTGACTCTGATTCTAGGTATCGACCCCGGCTcg                >  2:2950655/1‑150 (MQ=38)
                                                                                                                                   cccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                >  2:2950464/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTCCCGGTATTACGCCCTTTTCAATTGTGGTGCAGTACCTGTTGGAGGGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc               <  2:832156/150‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc               >  2:1627435/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc               >  1:674973/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc               >  2:4638969/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc               >  2:827862/1‑150 (MQ=18)
                                                                                                                                       cTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCa            >  1:1550922/1‑150 (MQ=18)
                                                                                                                                        tCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCAt           >  2:374930/1‑150 (MQ=255)
                                                                                                                                         ccTACGGGTTTAACGCCCTTTAAAAATGTGTTGACGTACCTTTTGGTGGGGTTTTCCCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATc          <  2:2576056/150‑1 (MQ=255)
                                                                                                                                          cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  2:5095011/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  2:4894957/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  1:2887395/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  1:4520729/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  2:5094645/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  2:3175540/1‑150 (MQ=18)
                                                                                                                                          cTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa         >  1:2944928/1‑150 (MQ=18)
                                                                                                                                           tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc        >  1:3414545/1‑150 (MQ=21)
                                                                                                                                           tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc        >  1:3414434/1‑150 (MQ=21)
                                                                                                                                           tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc        >  1:576650/1‑150 (MQ=21)
                                                                                                                                           tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc        >  2:1377455/1‑150 (MQ=21)
                                                                                                                                            aCAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc       >  1:3271563/1‑150 (MQ=255)
                                                                                                                                            aCAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc       >  1:4274043/1‑150 (MQ=255)
                                                                                                                                            aCAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc       >  1:1023941/1‑150 (MQ=255)
                                                                                                                                             ccgttATTACGCCCTTTTCAACTGTGGTGCAGTACCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCg      <  2:2104052/146‑1 (MQ=255)
                                                                                                                                             cAGGTATTACGCCCATTTCAATTGTGGTGCAGTACCTTTTGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCg      <  2:2104096/150‑1 (MQ=255)
                                                                                                                                              aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg     >  2:489351/1‑150 (MQ=255)
                                                                                                                                              aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg     >  1:5070940/1‑150 (MQ=255)
                                                                                                                                                 tATTACGCCCATTTCAAATGTGCTGCAGTACCGGGAGGAGCCGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTa  >  2:564124/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTA  >  NC_002947/1392774‑1393068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: