Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 10.9% 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,5120GA10.9% 383.3 / 29.5 128M1M (GTG→ATG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (60/54);  new base A (8/6);  total (68/60)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 137 total aligned reads displayed.

GCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291364‑4291660
                                                                                                                                                    |                                                                                                                                                    
     tCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                                                                   <  3:1594798/147‑1 (MQ=255)
     tCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                                                                   >  4:1594798/1‑147 (MQ=255)
             gCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                                                                         <  2:547879/149‑1 (MQ=255)
              cGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGt                                                                                                                                          <  4:1747975/147‑1 (MQ=255)
                tgtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGAc                                                                                                                                      <  2:1776607/149‑1 (MQ=255)
                    gATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                                                                   <  4:563685/132‑1 (MQ=255)
                    gATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                                                                   >  3:563685/1‑132 (MQ=255)
                         tgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  >  4:644571/1‑128 (MQ=255)
                             cTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           <  3:134123/147‑1 (MQ=255)
                              tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                                                                        >  3:652974/1‑149 (MQ=255)
                                    aTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                  >  1:1853593/1‑149 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  <  1:600368/116‑1 (MQ=255)
                                     ttCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  >  2:600368/1‑116 (MQ=255)
                                      tctcTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGcgcgaact                                                                                                                >  2:1765935/1‑145 (MQ=255)
                                          tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                            >  4:537981/1‑149 (MQ=255)
                                           aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                                                                            <  2:946516/116‑1 (MQ=255)
                                           aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                                                                            >  1:946516/1‑116 (MQ=255)
                                               cGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                                                                  <  4:1559059/122‑1 (MQ=255)
                                                  aCAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACTCCGTTGACGGCCTTCGCgg                                                                                                                           <  2:671130/126‑1 (MQ=255)
                                                  aCAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACTCCGTTGACGGCCTTCGCgg                                                                                                                           >  1:671130/1‑126 (MQ=255)
                                                      ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           <  1:1145246/122‑1 (MQ=255)
                                                           ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGa                                                                                           <  4:1210852/149‑1 (MQ=255)
                                                            gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  >  4:134620/1‑93 (MQ=255)
                                                            gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  <  3:1605839/93‑1 (MQ=255)
                                                            gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  <  3:134620/93‑1 (MQ=255)
                                                             tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGa                                                                                                                         <  2:377044/117‑1 (MQ=255)
                                                             tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGa                                                                                                                         >  1:377044/1‑117 (MQ=255)
                                                                  tAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                    >  1:314227/1‑149 (MQ=255)
                                                                    aGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           <  1:713102/108‑1 (MQ=255)
                                                                    aGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGa                                                                                           >  4:262227/1‑140 (MQ=255)
                                                                    aGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGa                                                                                           <  3:262227/140‑1 (MQ=255)
                                                                     gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCa                                                                                 <  2:314227/149‑1 (MQ=255)
                                                                        aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                                                                              <  3:424489/149‑1 (MQ=255)
                                                                        aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                                                                              >  3:1037752/1‑149 (MQ=255)
                                                                         gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                                                                    >  3:1414490/1‑126 (MQ=255)
                                                                         gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                                                                    <  4:1414490/126‑1 (MQ=255)
                                                                         gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             <  4:1037752/149‑1 (MQ=255)
                                                                             agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCAtt                                                                         >  1:1238689/1‑149 (MQ=255)
                                                                               agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                                                                       >  4:1272937/1‑149 (MQ=255)
                                                                                 gtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                                                                         >  1:413002/1‑113 (MQ=255)
                                                                                 gtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                                                                         <  2:413002/113‑1 (MQ=255)
                                                                                     tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                    <  1:489910/130‑1 (MQ=255)
                                                                                     tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                    >  2:489910/1‑130 (MQ=255)
                                                                                     tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa                                                                 >  3:688893/1‑149 (MQ=255)
                                                                                         tGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                                                                   <  2:769840/79‑1 (MQ=255)
                                                                                         tGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGCGTTGACCCCGTTGACGGc                                                                                                                                   >  1:769840/1‑79 (MQ=255)
                                                                                          ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg                                                            >  1:1135393/1‑149 (MQ=255)
                                                                                           gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                           <  2:1853593/149‑1 (MQ=255)
                                                                                              gTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                          >  4:1090217/1‑115 (MQ=255)
                                                                                              gTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                          <  3:1090217/115‑1 (MQ=255)
                                                                                                 ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTGGGGGTCAACgtgt                                                                                                                                                   >  2:603086/1‑55 (MQ=255)
                                                                                                 ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  >  2:1832823/1‑56 (MQ=255)
                                                                                                 ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  <  1:1832823/56‑1 (MQ=255)
                                                                                                   ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                                                                       >  2:89805/1‑113 (MQ=255)
                                                                                                   ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                                                                       <  1:89805/113‑1 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAc                                                                                                                                               >  1:550475/1‑56 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgat                                                                                               <  1:601539/104‑1 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc                                                         >  3:796760/1‑142 (MQ=255)
                                                                                                    cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgct                                                  <  3:537981/149‑1 (MQ=255)
                                                                                                     ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGATTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                    <  2:584825/114‑1 (MQ=255)
                                                                                                        ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                              <  2:1135393/149‑1 (MQ=255)
                                                                                                        ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTACTGCTGCTc                                              <  2:1238689/149‑1 (MQ=255)
                                                                                                         tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGCGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                >  4:1792041/1‑98 (MQ=255)
                                                                                                         tGAAGTCAGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                <  3:1792041/98‑1 (MQ=255)
                                                                                                                ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                     <  4:722264/118‑1 (MQ=255)
                                                                                                                 gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                     >  1:1721863/1‑149 (MQ=255)
                                                                                                                 gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                     <  4:688893/149‑1 (MQ=255)
                                                                                                                  ttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCg                                    <  3:1272937/149‑1 (MQ=255)
                                                                                                                  ttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCg                                    <  3:119842/149‑1 (MQ=255)
                                                                                                                           gttgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGt                           >  1:1654453/1‑149 (MQ=255)
                                                                                                                            ttgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTa                          >  3:733616/1‑149 (MQ=255)
                                                                                                                              gttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTg                                       >  1:1281762/1‑134 (MQ=255)
                                                                                                                                  ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                          >  3:986743/1‑79 (MQ=255)
                                                                                                                                  ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                    >  4:1558743/1‑149 (MQ=255)
                                                                                                                                  ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                    <  2:1241353/149‑1 (MQ=255)
                                                                                                                                    aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                  <  4:915365/149‑1 (MQ=255)
                                                                                                                                    aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                  <  2:1371652/149‑1 (MQ=255)
                                                                                                                                       tctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca               >  1:804884/1‑149 (MQ=255)
                                                                                                                                       tctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca               >  4:402076/1‑149 (MQ=255)
                                                                                                                                         tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAtt                        <  4:1306236/138‑1 (MQ=255)
                                                                                                                                         tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc             <  2:1721863/149‑1 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                                                                 <  2:762293/47‑1 (MQ=37)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                                                                    <  1:700584/92‑1 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                                                                    >  2:700584/1‑92 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTg                                       >  1:941497/1‑121 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTa                          >  3:1628558/1‑134 (MQ=255)
                                                                                                                                           ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGAACTGACGTa                          <  4:1628558/134‑1 (MQ=255)
                                                                                                                                            gggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGc          <  4:733616/149‑1 (MQ=255)
                                                                                                                                            gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                     <  2:304233/122‑1 (MQ=255)
                                                                                                                                            gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGc                                     >  1:304233/1‑122 (MQ=255)
                                                                                                                                             ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCg         >  3:282257/1‑149 (MQ=255)
                                                                                                                                             ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCg         >  1:129325/1‑149 (MQ=255)
                                                                                                                                             ggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCg         >  3:1777082/1‑149 (MQ=255)
                                                                                                                                                 aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                      <  3:874893/100‑1 (MQ=255)
                                                                                                                                                 aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                      >  4:874893/1‑100 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGa                                                                                           <  1:1288378/62‑1 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGa                                                                                           >  2:1288378/1‑62 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc    >  2:1414658/1‑149 (MQ=255)
                                                                                                                                                  aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc    >  3:229293/1‑149 (MQ=255)
                                                                                                                                                    gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc  <  1:1414658/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291364‑4291660

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