Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,442 Δ1 bp 8.2% intergenic (+14/+65) ftnB → / ← yecJ ferritin B, putative ferrous iron reservoir/DUF2766 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4420T. 8.2% 256.5 / 16.3 85intergenic (+14/+65)ftnB/yecJferritin B, putative ferrous iron reservoir/DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base T (38/40);  new base . (3/4);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTG  >  NC_000913/1987296‑1987579
                                                                                                                                                  |                                                                                                                                          
aaaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                          <  3:280461/149‑1 (MQ=255)
 aaaGAACAGGAGAATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAAt                                                                                                                                         <  2:84908/149‑1 (MQ=255)
      aCAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAtt                                                                                                                                    >  6:262366/1‑149 (MQ=255)
       cagcagCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑g                                                                                                                                  >  5:381778/1‑149 (MQ=255)
           agcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGcc                                                                                                                              <  8:119773/149‑1 (MQ=255)
              atgatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt                                                                                                                           >  1:41352/1‑149 (MQ=255)
               tgatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt                                                                                                                           >  1:253089/1‑148 (MQ=255)
                       tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                  >  1:413869/1‑149 (MQ=255)
                         ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTga                                                                                                                <  2:253089/149‑1 (MQ=255)
                         ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAACGtttt                                                                                                                  >  6:274687/1‑147 (MQ=255)
                         ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGAGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                  <  8:361668/147‑1 (MQ=255)
                             tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                            >  6:251640/1‑149 (MQ=255)
                             tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                            >  5:309716/1‑149 (MQ=255)
                               cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                            <  5:262366/147‑1 (MQ=255)
                               cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAAcc                                                                                                          >  6:36103/1‑149 (MQ=255)
                                 aaCCATTCTTGATGAAGTGCGCAGTGCGAAAATAGCGCGTATGTGCCCTGTGCACACCGACCAACATGTTCTGAAAGTCGTGTCACCCCAGCTGCATTGATCATCATCGGCGCTAAAGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGc                                                                                                        >  5:109152/1‑149 (MQ=255)
                                    cATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGc                                                                                                     >  4:252262/1‑149 (MQ=255)
                                         ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                >  4:373968/1‑149 (MQ=255)
                                           gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                <  2:113264/147‑1 (MQ=255)
                                                        gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGCCCGACCAACAGGTTCTGAATGTCGTGTCAGCGCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGCTTTTGAGAAACCGCTCCCTCCTCGGATAGAagcaggg                                                                                 >  5:92555/1‑147 (MQ=255)
                                                         tGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCgg                                                                                >  6:49028/1‑149 (MQ=255)
                                                              aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGaa                                                                                       <  8:184953/137‑1 (MQ=255)
                                                              aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGaa                                                                                       >  7:184953/1‑137 (MQ=255)
                                                              aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttt                                                                           <  2:413869/149‑1 (MQ=255)
                                                              aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttt                                                                           <  4:328282/149‑1 (MQ=255)
                                                               aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                          <  6:381778/149‑1 (MQ=255)
                                                               aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                          <  3:252262/149‑1 (MQ=255)
                                                                cTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt                                                                         <  2:12667/149‑1 (MQ=255)
                                                                 ttGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt                                                                         >  3:6386/1‑148 (MQ=255)
                                                                    cGGGTATGTGGCCTGTGCAGACCGACCAACTTGTTCTGAATGTCGTGTAACACCAGCTGCATTGATCATCTTGGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGTAGCAGCTGTTTTTTTAAt                                                                      <  6:109152/148‑1 (MQ=255)
                                                                    cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  8:240597/148‑1 (MQ=255)
                                                                    cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  3:326946/148‑1 (MQ=255)
                                                                    cGGGGATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGACGCAGCGGTTTTTTTAATg                                                                     <  2:418225/149‑1 (MQ=255)
                                                                     gggTATGTGCCCTGTGCAGACCGAGCAACATGTTCTGAATGTCGTGTCACACCAGCTGCTTTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg                                                                    <  6:92555/149‑1 (MQ=255)
                                                                      ggTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg                                                                   <  3:373968/149‑1 (MQ=255)
                                                                       gTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGa                                                                  <  7:302218/149‑1 (MQ=255)
                                                                          tgtcCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTc                                                               <  5:49028/149‑1 (MQ=255)
                                                                           gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                              >  5:204291/1‑149 (MQ=255)
                                                                            tgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAc                                                             >  7:332384/1‑149 (MQ=255)
                                                                             gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                                                            >  6:341542/1‑149 (MQ=255)
                                                                                   tgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                       <  3:223043/148‑1 (MQ=255)
                                                                                   tgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgg                                                      >  2:329459/1‑149 (MQ=255)
                                                                                        accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttt                                                                           <  3:354615/123‑1 (MQ=255)
                                                                                        accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttt                                                                           >  4:354615/1‑123 (MQ=255)
                                                                                        accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTa                                                  >  6:242074/1‑148 (MQ=255)
                                                                                          cgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAAc                                               <  5:36103/149‑1 (MQ=255)
                                                                                              cAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGGAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                           <  2:41352/149‑1 (MQ=255)
                                                                                              cAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                           >  3:81075/1‑149 (MQ=255)
                                                                                               aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg                                          <  5:242074/149‑1 (MQ=255)
                                                                                                aCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAt                                                                                                                                     <  2:288001/58‑1 (MQ=255)
                                                                                                aCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAt                                                                                                                                     >  1:288001/1‑58 (MQ=255)
                                                                                                 cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt                                        <  5:341542/149‑1 (MQ=255)
                                                                                                  aTGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt                                        >  4:208729/1‑148 (MQ=255)
                                                                                                     ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                    <  6:204291/149‑1 (MQ=255)
                                                                                                       cTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATa                                   <  4:6386/148‑1 (MQ=255)
                                                                                                        tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGttt                                                                                                                   >  4:213687/1‑67 (MQ=255)
                                                                                                        tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGttt                                                                                                                   <  3:213687/67‑1 (MQ=255)
                                                                                                          aaTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTATACTTGAGTTCAATAAgc                                <  6:309716/148‑1 (MQ=255)
                                                                                                          aaTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgcg                               >  7:108129/1‑149 (MQ=255)
                                                                                                            tGTCGTGTCACACCAGCTGCATTGATCATCATCTGCGCTAATGGATT‑GCGACGATGAAGGTTTTGAGAAACCGCTGCCTCATGTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCTCCCTGTGGGGTTTACTTGTGTTCAAt                                    >  5:291430/1‑142 (MQ=255)
                                                                                                            tGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATA‑GCGCCTATGAAGGTTTTGAGAAACCGCTGCCACATCTGTTTTAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                    <  6:291430/142‑1 (MQ=255)
                                                                                                                 tgtCACACCAGCTGCTTTGATCATCATCGGCGCTTATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGTAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGc                        <  5:274687/149‑1 (MQ=255)
                                                                                                                 tgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGc                        <  5:251640/149‑1 (MQ=255)
                                                                                                                     acacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                    >  1:228723/1‑149 (MQ=255)
                                                                                                                       acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                  >  5:233735/1‑149 (MQ=255)
                                                                                                                          aGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt               <  8:108129/149‑1 (MQ=255)
                                                                                                                           gCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                            <  8:406473/135‑1 (MQ=255)
                                                                                                                           gCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                            >  7:406473/1‑135 (MQ=255)
                                                                                                                              gCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGcc           >  4:81108/1‑149 (MQ=255)
                                                                                                                               cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGccc          >  1:123902/1‑149 (MQ=255)
                                                                                                                               cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGc            <  3:208729/147‑1 (MQ=255)
                                                                                                                                 ttGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg         >  6:222693/1‑148 (MQ=255)
                                                                                                                                    atcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCAt      <  6:233735/148‑1 (MQ=255)
                                                                                                                                    atcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATg     >  4:114633/1‑149 (MQ=255)
                                                                                                                                     tcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  5:297689/83‑1 (MQ=255)
                                                                                                                                     tcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      >  6:297689/1‑83 (MQ=255)
                                                                                                                                     tcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATgg    <  2:228723/149‑1 (MQ=255)
                                                                                                                                      cttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  4:176761/80‑1 (MQ=37)
                                                                                                                                      cttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  4:170658/80‑1 (MQ=37)
                                                                                                                                      cttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      >  3:176761/3‑82 (MQ=37)
                                                                                                                                      cttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      >  3:170658/3‑82 (MQ=37)
                                                                                                                                       atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGt   >  5:364399/1‑148 (MQ=255)
                                                                                                                                       atcatcGGCGCTAATGAATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCCTCTGTTTGAAGCAGCGGTTTTTTTAATTGGATTCACCCTGTGGGGTAAAATTGAGTTCAATAAGCTCGATGGCTTTTTGTATTGCCCGCATGGTg  >  6:168259/1‑149 (MQ=255)
                                                                                                                                          atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                           <  2:20159/105‑1 (MQ=39)
                                                                                                                                          atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                           >  1:20159/1‑105 (MQ=39)
                                                                                                                                                  |                                                                                                                                          
AAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTG  >  NC_000913/1987296‑1987579

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑