Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,162 T→A 6.2% intergenic (+5/+66) smrB → / ← yfcO putative DNA endonuclease/DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,1620TA 6.2% 229.1 / 4.8 64intergenic (+5/+66)smrB/yfcOputative DNA endonuclease/DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (24/36);  new base A (2/2);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGG  >  NC_000913/2449015‑2449226
                                                                                                                                                   |                                                                
ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCATCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGtt                                                                 <  3:1248967/149‑1 (MQ=255)
ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGtt                                                                 >  2:1870126/1‑149 (MQ=255)
ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGtt                                                                 <  2:607384/149‑1 (MQ=255)
ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGtt                                                                 <  1:1688927/149‑1 (MQ=255)
 tttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTg                                                                <  2:1794282/149‑1 (MQ=255)
    gAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGTCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAgc                                                             >  1:1121958/1‑149 (MQ=255)
     aaGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAgcg                                                            <  1:1870126/149‑1 (MQ=255)
       gCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGcgcg                                                          <  4:1454960/149‑1 (MQ=255)
        caacaaACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGt                                                         <  1:354468/149‑1 (MQ=255)
         aacaaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGtt                                                        >  3:882377/1‑149 (MQ=255)
         aacaaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGATGCCGCCGGAGTTGCCCTGAGGAGTTGa                                                               <  3:1463976/142‑1 (MQ=255)
         aacaaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGATGCCGCCGGAGTTGCCCTGAGGAGTTGa                                                               >  4:1463976/1‑142 (MQ=255)
           caaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTgc                                                      >  1:1321950/1‑149 (MQ=255)
                acCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAg                                                 >  2:600673/1‑149 (MQ=255)
                    cTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTAc                                             <  4:882377/149‑1 (MQ=255)
                     tgtgGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACg                                            <  2:112302/149‑1 (MQ=255)
                          ctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTaa                                       <  2:1321950/149‑1 (MQ=255)
                           tggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAAt                                      <  1:477579/149‑1 (MQ=255)
                              cacaACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATg                                     <  2:158383/147‑1 (MQ=255)
                              cacaACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATg                                     >  1:158383/1‑147 (MQ=255)
                               acaacaTCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCa                                  <  2:1189005/149‑1 (MQ=255)
                                 aacaTCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgg                                <  2:101517/149‑1 (MQ=255)
                                 aacaTCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGATGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgg                                >  1:486124/1‑149 (MQ=255)
                                 aacaTCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGATGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgg                                >  1:224589/1‑149 (MQ=255)
                                      cggCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGa                           <  2:1121958/147‑1 (MQ=255)
                                       cGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGATGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          <  2:486124/149‑1 (MQ=255)
                                       cGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGATGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          <  2:224589/149‑1 (MQ=255)
                                        gCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          <  3:1792674/148‑1 (MQ=255)
                                        gCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          >  4:1792674/1‑148 (MQ=255)
                                        gCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                         <  4:282470/149‑1 (MQ=255)
                                          aTGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGc                       <  2:611086/149‑1 (MQ=255)
                                           tgtgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcc                      <  3:1518388/149‑1 (MQ=255)
                                            gtgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGccc                     <  1:1028685/149‑1 (MQ=255)
                                               aTGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGCTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTgcgc                                                    >  2:44335/1‑115 (MQ=255)
                                               aTGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGCTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTgcgc                                                    <  1:44335/115‑1 (MQ=255)
                                                 ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAg                                                 >  4:720055/1‑116 (MQ=255)
                                                 ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAg                                                 >  4:683263/1‑116 (MQ=255)
                                                 ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAg                                                 <  3:720055/116‑1 (MQ=255)
                                                 ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAg                                                 <  3:683263/116‑1 (MQ=255)
                                                 ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcc                <  1:600673/149‑1 (MQ=255)
                                                     tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgg                                <  4:612307/129‑1 (MQ=255)
                                                     tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgg                                >  3:612307/1‑129 (MQ=255)
                                                         atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAAt                                      >  4:245627/1‑119 (MQ=255)
                                                         atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAAt                                      <  3:245627/119‑1 (MQ=255)
                                                          tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTg                                                                <  2:68792/92‑1 (MQ=255)
                                                          tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTg                                                                >  1:68792/1‑92 (MQ=255)
                                                             ggCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          >  2:899581/1‑127 (MQ=255)
                                                             ggCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                          <  1:899581/127‑1 (MQ=255)
                                                              gCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCa                                                  <  2:48161/102‑1 (MQ=255)
                                                              gCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCa                                                  >  1:48161/1‑102 (MQ=255)
                                                                   aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgggg                              <  3:168871/117‑1 (MQ=255)
                                                                   aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAgggg                              >  4:168871/1‑117 (MQ=255)
                                                                    aaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                         >  2:1268225/1‑121 (MQ=255)
                                                                    aaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                         <  1:1268225/121‑1 (MQ=255)
                                                                            tGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg  >  1:1191408/1‑136 (MQ=255)
                                                                            tGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg  <  2:1191408/136‑1 (MQ=255)
                                                                                      gcggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTaa                                       <  3:1663169/89‑1 (MQ=255)
                                                                                      gcggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTaa                                       >  4:1663169/1‑89 (MQ=255)
                                                                                               ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAgc                                                             <  4:1708539/58‑1 (MQ=255)
                                                                                               ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAgc                                                             >  3:1708539/1‑58 (MQ=255)
                                                                                               ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTaa                                       >  1:1011360/1‑80 (MQ=255)
                                                                                               ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTaa                                       <  2:1011360/80‑1 (MQ=255)
                                                                                                         ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg  >  2:1248269/1‑107 (MQ=255)
                                                                                                         ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg  <  1:1248269/107‑1 (MQ=255)
                                                                                                               tCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGc                             <  3:942229/74‑1 (MQ=255)
                                                                                                               tCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGc                             >  4:942229/1‑74 (MQ=255)
                                                                                                                                           cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGc                             <  1:1688020/46‑1 (MQ=255)
                                                                                                                                           cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGc                             >  2:1688020/1‑46 (MQ=255)
                                                                                                                                                   |                                                                
TTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGG  >  NC_000913/2449015‑2449226

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