Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 11.4% M1M (GTG→ATG) †‡ nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,5120GA11.4% 203.6 / 15.0 70M1M (GTG→ATG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (5/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291364‑4291660
                                                                                                                                                    |                                                                                                                                                    
gCCAGTCTGGTCTGTGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                                                                      >  3:66803/1‑149 (MQ=255)
 ccAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTATGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                                                                      >  1:592290/1‑148 (MQ=255)
  cAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgt                                                                                                                                                     >  1:672125/1‑148 (MQ=255)
      cTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                                                                <  3:864167/149‑1 (MQ=255)
              cGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                                                                         >  1:620906/1‑148 (MQ=255)
                     aTGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGAc                                                                                                                                      >  4:869396/1‑144 (MQ=255)
                     aTGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGAc                                                                                                                                      <  3:869396/144‑1 (MQ=255)
                       gctgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTc                                                                                                                               >  3:772578/1‑149 (MQ=255)
                          gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcgcg                                                                                                                            >  4:450299/1‑149 (MQ=255)
                            gcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                          >  3:464416/1‑149 (MQ=255)
                            gcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           >  1:502170/1‑148 (MQ=255)
                             cTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                          <  1:735978/148‑1 (MQ=255)
                             cTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                          <  2:672125/148‑1 (MQ=255)
                              tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                                                                        >  4:284800/1‑149 (MQ=255)
                                tttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTAACCCCGTTGACGGCCTTCGCGGGAGTg                                                                                                                      <  3:603288/149‑1 (MQ=255)
                                    aTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                                                                   <  4:52705/132‑1 (MQ=255)
                                    aTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                                                                   >  3:52705/1‑132 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                  >  2:579287/1‑148 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                  <  2:502170/148‑1 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                          <  4:347808/124‑2 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                          >  3:166591/1‑123 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                          >  3:347808/1‑123 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                          <  4:166591/124‑2 (MQ=255)
                                       ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                                                                                >  2:1825/1‑148 (MQ=255)
                                         ctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                             <  4:451/149‑1 (MQ=255)
                                         ctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                              <  1:419049/148‑1 (MQ=255)
                                           aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                            <  2:592290/148‑1 (MQ=255)
                                           aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                                                                           <  3:409990/149‑1 (MQ=255)
                                               cGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCt                                                                                                       <  4:772578/149‑1 (MQ=255)
                                               cGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCt                                                                                                       <  3:251514/149‑1 (MQ=255)
                                                ggACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                                                                 <  4:239807/138‑1 (MQ=255)
                                                ggACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                                                                 >  3:239807/1‑138 (MQ=255)
                                                                  tAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                    >  4:485096/1‑149 (MQ=255)
                                                                      ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  3:485096/149‑1 (MQ=255)
                                                                      ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  3:450299/149‑1 (MQ=255)
                                                                       cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                                                                     >  3:484900/1‑95 (MQ=255)
                                                                       cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                                                                     <  4:484900/95‑1 (MQ=255)
                                                                       cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATc                                                                               <  3:235320/149‑1 (MQ=255)
                                                                       cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                                                                           <  1:153661/89‑1 (MQ=255)
                                                                       cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                                                                           >  2:153661/1‑89 (MQ=255)
                                                                          aaaaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             >  2:589394/1‑148 (MQ=255)
                                                                                    gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                                                                   >  2:654206/1‑148 (MQ=255)
                                                                                           gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                           >  3:244536/1‑149 (MQ=255)
                                                                                            aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                           >  1:696370/1‑148 (MQ=255)
                                                                                             cGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc                                                         >  4:285873/1‑149 (MQ=255)
                                                                                                aTTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                      <  4:464416/149‑1 (MQ=255)
                                                                                                 ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                      >  1:333753/1‑148 (MQ=255)
                                                                                                  tccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                     >  2:741579/1‑148 (MQ=255)
                                                                                                   ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc                                                   >  4:680149/1‑149 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             >  3:883014/1‑122 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             <  4:883014/122‑1 (MQ=255)
                                                                                                     ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt                                                              <  1:712425/136‑1 (MQ=255)
                                                                                                     ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt                                                              >  2:712425/1‑136 (MQ=255)
                                                                                                      tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgctg                                                 <  2:314371/148‑1 (MQ=255)
                                                                                                           aaGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGcc                                           <  4:234577/149‑1 (MQ=255)
                                                                                                           aaGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGc                                            <  1:62138/148‑1 (MQ=255)
                                                                                                              tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCg                                                                                                                    >  3:697847/1‑73 (MQ=255)
                                                                                                              tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCg                                                                                                                    <  4:697847/73‑1 (MQ=255)
                                                                                                                  ttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  2:344864/105‑1 (MQ=255)
                                                                                                                  ttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                >  1:344864/1‑105 (MQ=255)
                                                                                                                    ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGt                                   <  2:696370/148‑1 (MQ=255)
                                                                                                                    ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTa                                  <  4:508030/149‑1 (MQ=255)
                                                                                                                                   tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt                   <  3:285873/149‑1 (MQ=255)
                                                                                                                                   tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                    >  1:254425/1‑148 (MQ=255)
                                                                                                                                      gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAtt                        <  2:457125/141‑1 (MQ=255)
                                                                                                                                      gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAtt                        >  1:457125/1‑141 (MQ=255)
                                                                                                                                       caacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtctc                >  1:257352/4‑148 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAg           >  3:258028/1‑149 (MQ=255)
                                                                                                                                              gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGACTGCGGTTGCGCTGGCATGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGa        <  4:293071/149‑1 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAtt      <  1:644386/147‑1 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAtt      >  2:644386/1‑147 (MQ=255)
                                                                                                                                                   cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc    >  1:249078/1‑148 (MQ=255)
                                                                                                                                                   cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc    >  1:601748/1‑148 (MQ=255)
                                                                                                                                                    gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTTACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct   <  1:579287/148‑1 (MQ=255)
                                                                                                                                                    gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc  <  3:284800/149‑1 (MQ=255)
                                                                                                                                                    gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct   <  1:1825/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291364‑4291660

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