Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,290,993 A→C Y173S (TAT→TCT)  nrfD → formate‑dependent nitrite reductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,290,9930AC95.2% 192.9 / ‑3.7 64Y173S (TAT→TCT) nrfDformate‑dependent nitrite reductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base C (27/33);  minor base T (2/1);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CAACAACGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTATCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCGTATGGAAA  >  NC_000913/4290845‑4291139
                                                                                                                                                    |                                                                                                                                                   
caacaaCGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGtctc                                                                                                                                                   >  1:856246/1‑151 (MQ=255)
caacaaCGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGtctc                                                                                                                                                   <  1:1087943/151‑1 (MQ=255)
   caacGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGt                                                                                                                                                <  1:773659/151‑1 (MQ=255)
    aacGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGGAATCGTCTCCGtt                                                                                                                                               >  2:588631/1‑151 (MQ=255)
    aacGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGtt                                                                                                                                               <  2:190958/151‑1 (MQ=255)
     acGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTc                                                                                                                                              <  2:1001097/151‑1 (MQ=255)
      cGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTcc                                                                                                                                             <  1:675651/151‑1 (MQ=255)
       gCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGGCTATACCGGCCTTCTGCTGTCGGCGCTGAAATTGTCTCCGTTcc                                                                                                                                             >  1:593180/1‑150 (MQ=255)
               ccGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAAcc                                                                                                                                    <  2:883663/151‑1 (MQ=255)
                       ggggATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCt                                                                                                                             <  1:67324/150‑1 (MQ=255)
                       ggggATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCt                                                                                                                             <  2:555994/150‑1 (MQ=255)
                        gggATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTg                                                                                                                            >  1:580906/1‑150 (MQ=255)
                           aTCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTTGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGGAATCGTCTCCGTTCCTCAATAACCCGGTCCTGCCgg                                                                                                                        >  2:853946/1‑151 (MQ=255)
                              gTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTTGGGGCTTCTACCGGCTTTTTGCTGTCGGCGCTGCAATCGGCTCCGTTCCTCAATAACCCGATCCTGCCCGTg                                                                                                                      >  1:420144/1‑150 (MQ=255)
                              gTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGc                                                                                                                     <  1:794000/151‑1 (MQ=255)
                                     aGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCct                                                                                                              <  1:352147/151‑1 (MQ=255)
                                         tCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTCCTGGCGGCGTTCTTGGGGGCCTTTACCGGCTTTCTGCTGAAGGGGCTTAACCCATTGCAGTTTCTCAATAACACCATCCTTCCGGTGCTATTCctcctg                                                                                                          >  2:854726/1‑148 (MQ=255)
                                          cTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTTTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCtctt                                                                                                           >  1:123454/1‑149 (MQ=255)
                                            gAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCtcttct                                                                                                         <  1:451278/149‑1 (MQ=255)
                                               ccTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGGGCTGTTCCTCTTCTCCGGc                                                                                                    >  1:1023948/1‑151 (MQ=255)
                                                  gctgACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCCGCAt                                                                                                  >  2:82066/1‑150 (MQ=255)
                                                    tgACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTTTTCCCGCTTTTTGCTGTCGGCCCTGAAATCGTCTCCGTTCCTCAAATACCCGCTCCTTCCGGGGCTGTTTCTCTTCTCCGGcct                                                                                                  >  1:575689/1‑146 (MQ=255)
                                                        gCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAGCCCGATCCTGCCAGTGCTGTTCCTCTTCTCCGGCATCTCGTc                                                                                            <  2:160787/150‑1 (MQ=255)
                                                        gCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTcc                                                                                           <  2:483148/151‑1 (MQ=255)
                                                          ccGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCgg                                                                                         <  1:217709/151‑1 (MQ=255)
                                                          ccGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCgg                                                                                         <  1:719668/151‑1 (MQ=255)
                                                          ccGTTCATCGCGGACTGGAAACATTGATGCAGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCgg                                                                                         <  1:678403/151‑1 (MQ=255)
                                                           cGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGGGGTGTTGGGGGGTTTTAACGGGTTTTTTGTGTCGGCGGTGGAATCGGCTCCGTTCCTCAAAAACCCCAACCTGCCGGGGGCGGTCCTCCTCTTCGGGATTTCTTCCg                                                                                          >  2:1196267/1‑149 (MQ=255)
                                                            gTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGt                                                                                        <  2:794404/150‑1 (MQ=255)
                                                               cATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTgcgg                                                                                    <  1:45171/151‑1 (MQ=255)
                                                               cATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTgcgg                                                                                    <  2:358383/151‑1 (MQ=255)
                                                                aTCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTgcg                                                                                     <  1:836046/149‑1 (MQ=255)
                                                                  cgcgGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCGTATACCGGCTTTCTGCTGGCGGGGGTGAAGTCGTCTCCGGTCCTCAATAACCCGATCCTGCCCGTGCTGTTCCTCTTCTCCGGCACCTCGTCCGGTGcggcgg                                                                                 >  2:909617/1‑151 (MQ=255)
                                                                   gcgGTCTGGAAAACTTGATGCTGGTGCTAGGGGTGGTGTAGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAAGCGCTTCCGTTAGTCACGAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCTTCCGGTgcggag                                                                                  <  2:319824/149‑3 (MQ=255)
                                                                   gcgGACTGGAAACATTGATGCTCGTGCTGGCGGAGTTGTTGGGGGCTTCCACCGGCTTTCTGCTGTCGGCGCCGGAATCCTCTCCGGTCCCTAAAAACCCGCTCCTGCCCGCGCTGCTCCTCTTCTCCGGCATAGCGTCCCGTTtggcgg                                                                                 >  2:105196/1‑150 (MQ=255)
                                                                       aCTGGAAACATTGCTGCTGGTGCTGGCGGTGGTGTTGGGGGCTTCTCCCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCCATCCTGCCGGGGCTGTTCCTCTTCTCCGGCATCTCGTCCGGGGGGGCGgggg                                                                             >  2:303303/1‑149 (MQ=255)
                                                                                      gctggtgctggCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCg                                                            <  1:929629/151‑1 (MQ=255)
                                                                                        tggtgctggCGGCGGTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTTTCCGTTCCTCAATAACCCGATCCCGCCGGTGGTGGTCCTCTTCTCCGGCATCTCCTCCGGGGCGGCGGGGGGGCTggt                                                                      >  2:474890/1‑137 (MQ=255)
                                                                                                gCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAAc                                                  <  1:964325/151‑1 (MQ=255)
                                                                                                  ggTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGGGGAGCTGAACGGCAAGGCGAttagccca                                                   >  2:1085940/1‑141 (MQ=255)
                                                                                                        gttgGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTaa                                           <  2:231401/150‑1 (MQ=255)
                                                                                                        gttgGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAAc                                          <  1:449849/151‑1 (MQ=255)
                                                                                                            ggggCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCCCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGc                                      <  1:552134/151‑1 (MQ=255)
                                                                                                             gggCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATACCCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTATCCCGCa                                     <  2:307857/151‑1 (MQ=255)
                                                                                                                   ataCCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGATGCCCCGGTCGTCCTGGCGGTACGCCAACGGCGGAACACGCATTcc                                 >  2:563057/1‑149 (MQ=255)
                                                                                                                     aCCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGGGCTGATCGCCCTGGGGGTACGCCAACGCAGTAACCCGCCTTCCACgg                             >  2:450919/1‑151 (MQ=255)
                                                                                                                     aCCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACg                              <  1:1068803/150‑1 (MQ=255)
                                                                                                                     aCCGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACg                              <  1:1086779/150‑1 (MQ=255)
                                                                                                                      ccGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGGGGCGCTGGGCGCCCCGGGGACACGCCCACGCAGTAAACCGCATTCCCCgg                             >  2:976276/1‑150 (MQ=255)
                                                                                                                       cGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGaa                           <  1:542634/151‑1 (MQ=255)
                                                                                                                       cGGCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCACTCCACgg                             >  2:353146/1‑149 (MQ=255)
                                                                                                                        ggCTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCCTGGCGGTACGCCAACGCAGTAACCCCCATTCCACg                              >  2:470518/1‑147 (MQ=255)
                                                                                                                          cTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAgc                         >  1:1058877/1‑150 (MQ=255)
                                                                                                                          cTTTCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGAGGGGCTGGTCGCCCTTGCGATACGCCAACGCAGTAAACCCCATTCCACGGaa                           >  1:168789/1‑148 (MQ=255)
                                                                                                                            ttCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGCATGCCCTTGCGATACGCCAACGGAGTAACCCGCACTCCACGGAAgcgc                       >  1:259358/1‑150 (MQ=255)
                                                                                                                             tCTGCTGTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGGTCGCCATGGCGATACGCCAACGCAGTAACCCGCCTTCCACGGGAgcgc                       >  1:563092/1‑149 (MQ=255)
                                                                                                                                  tgTCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGCGGGCGCGAGCGCCGTCGAAATACTGCAACGCAGTAACCCGCACTCCACgacggcccgcactg                 >  2:501063/1‑137 (MQ=255)
                                                                                                                                    tCGGCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGGAGCGCAGTTTGt                >  1:247933/1‑149 (MQ=255)
                                                                                                                                      ggCGCTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTacac            <  1:124143/151‑1 (MQ=255)
                                                                                                                                       gcgcTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACAcc           <  1:156993/151‑1 (MQ=255)
                                                                                                                                        cgcTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCg          <  1:884235/151‑1 (MQ=255)
                                                                                                                                         gcTGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCGt         <  1:450940/151‑1 (MQ=255)
                                                                                                                                           tGAAATCGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAAACCGCATTCCACGGAAGCGCAGTTTGTACACCGTa        >  1:593127/1‑150 (MQ=255)
                                                                                                                                            gAAATCGTCTCCGTTTCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCGTAt       <  1:893274/150‑1 (MQ=255)
                                                                                                                                            gAAATCGTATCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCGTa        <  1:1052202/149‑1 (MQ=255)
                                                                                                                                                 cGTCTCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGGGATACGCCAACGCAGTAAACCGCATTCCACGGGAGCGCAGTTTGTACACCGTATGGaaa  >  1:342864/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAACAACGCTGGTTGCCGAAGCTGGGGATCGTGCAAAAGGTTCTGAGCCTGCTGACGCCCGTTCATCGCGGACTGGAAACATTGATGCTGGTGCTGGCGGTGTTGTTGGGGGCTTATACCGGCTTTCTGCTGTCGGCGCTGAAATCGTATCCGTTCCTCAATAACCCGATCCTGCCGGTGCTGTTCCTCTTCTCCGGCATCTCGTCCGGTGCGGCGGTGGCGCTGATCGCCATGGCGATACGCCAACGCAGTAACCCGCATTCCACGGAAGCGCAGTTTGTACACCGTATGGAAA  >  NC_000913/4290845‑4291139

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