Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,373,850 T→G 8.8% E150D (GAA→GAC)  dapB ← 4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,373,8500TG 8.8% 199.1 / 9.2 68E150D (GAA→GAC) dapB4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base G (3/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAAC  >  NC_002947/5373703‑5373998
                                                                                                                                                   |                                                                                                                                                    
gCGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGAtt                                                                                                                                                      <  7:226444/148‑1 (MQ=255)
 cGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGAttt                                                                                                                                                     <  7:116730/148‑1 (MQ=255)
 cGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTc                                                                                                                                                    <  1:241861/149‑1 (MQ=255)
 cGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGCTGGTGCGCCTCGATGATTTc                                                                                                                                                    >  4:106840/1‑149 (MQ=255)
  gATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTc                                                                                                                                                    >  1:110699/1‑148 (MQ=255)
  gATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTc                                                                                                                                                    <  2:110699/148‑1 (MQ=255)
      aCGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGata                                                                                                                                                <  8:107373/148‑1 (MQ=255)
      aCGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGata                                                                                                                                                >  4:17227/1‑148 (MQ=255)
       cgcgcgCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATc                                                                                                                                              >  2:73113/1‑149 (MQ=255)
              cGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGAcatca                                                                                                                                       >  2:232040/1‑149 (MQ=255)
                    gCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGcc                                                                                                                                  <  5:249277/148‑1 (MQ=255)
                    gCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGcc                                                                                                                                  >  4:154766/1‑148 (MQ=255)
                     cccTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCa                                                                                                                                <  3:17227/149‑1 (MQ=255)
                     cccTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCa                                                                                                                                <  3:154766/149‑1 (MQ=255)
                      ccTTCACGTCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATTAGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAg                                                                                                                               >  2:309918/1‑149 (MQ=255)
                      ccTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAg                                                                                                                               >  6:242282/1‑149 (MQ=255)
                           aCGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCAccc                                                                                                                          <  1:232040/149‑1 (MQ=255)
                              gCCATACACGGCCACTTCCTGCAGATCCCGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGCCCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCcgcg                                                                                                                       >  5:339840/1‑149 (MQ=255)
                                       tgCCAATTCCTTCAGATCACGGCCCAGCGCATTGTCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCTCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATGTCGATATCGAAATCATCgaggc                                                                                                                                <  4:252632/130‑5 (MQ=255)
                                       ggCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATGTCGATATCGAAATCATCgaggc                                                                                                                                >  3:252632/1‑127 (MQ=255)
                                        gCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                             >  2:182348/1‑149 (MQ=255)
                                        gCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                             <  1:73113/149‑1 (MQ=255)
                                        gCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                             <  3:106840/149‑1 (MQ=255)
                                          cACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                             >  1:246769/1‑147 (MQ=255)
                                                cTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCagcagc                                                                                                      <  5:242282/148‑1 (MQ=255)
                                                 tGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATGTCGATATCGAAatcatc                                                                                                                                     <  2:216176/116‑1 (MQ=255)
                                                 tGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATGTCGATATCGAAatcatc                                                                                                                                     >  1:216176/1‑116 (MQ=255)
                                                   cAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCtt                                                                                                    >  3:253169/1‑147 (MQ=255)
                                                   cAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTg                                                                                                   <  2:27785/148‑1 (MQ=255)
                                                   cAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGa                                                                                                  <  2:365215/149‑1 (MQ=255)
                                                    aGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAg                                                                                                 <  4:253169/149‑1 (MQ=255)
                                                      aTCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGc                                                                                               >  6:18283/1‑149 (MQ=255)
                                                           ggCCCCGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAgg                                                                                          <  6:123342/149‑1 (MQ=255)
                                                             cccAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGtt                                                                                        >  2:136238/1‑149 (MQ=255)
                                                               cAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGa                                                                                      >  1:45427/1‑149 (MQ=255)
                                                                aGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGAc                                                                                     >  7:249936/1‑149 (MQ=255)
                                                                     aTTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCa                                                                                <  2:246769/149‑1 (MQ=255)
                                                                         gCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCt                                                                            <  4:87032/149‑1 (MQ=255)
                                                                           acccACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTg                                                                           >  4:233969/3‑148 (MQ=255)
                                                                            aCCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTg                                                                           >  3:178887/1‑147 (MQ=255)
                                                                               aCTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGtt                                                                      >  2:197352/1‑149 (MQ=255)
                                                                                 ttCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGGCCCCCACGCTGAAGTTgg                                                                    >  4:131592/1‑149 (MQ=255)
                                                                                    aCCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTggcgg                                                                 >  6:209408/1‑149 (MQ=255)
                                                                                      ccATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCGAGCACCCGCGCCGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTggcggc                                                                >  6:272851/1‑148 (MQ=255)
                                                                                      ccATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTggcggc                                                                >  5:304434/1‑148 (MQ=255)
                                                                                        aagagcagtgcagtACCCGACGGCGCATCGACCTTGTGACGGTGGTGCGCCTCGATCATTTCGACAACCACATCATCGCCCAGCACCCGCGCAGCCATATCCCGCAGCTAGAGGCGAAGCTTGACCCCCACCCTGCAGTTGGCGGCGaa                                                             >  7:355900/5‑149 (MQ=255)
                                                                                               gtgcagtACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGAtt                                                      <  3:233969/149‑1 (MQ=255)
                                                                                                    gtACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAAt                                                 <  1:197352/149‑1 (MQ=255)
                                                                                                        ccGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAAt                                                 <  6:311577/145‑1 (MQ=255)
                                                                                                        ccGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAAt                                                 >  5:311577/1‑145 (MQ=255)
                                                                                                              gcgcATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCa                                       <  1:182348/149‑1 (MQ=255)
                                                                                                               cgcATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACACGAGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAAGAAGAAGTAGGCGGCGAACACGATTGGAATGTCCTTGCCAg                                      >  8:205431/1‑149 (MQ=255)
                                                                                                                  aTCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGGGCAGCCATATCCAGCAGCTAGAGGCTAAGGTTGACACCCACGCTGAAGTTGGCGGCGACCACGATTGGAATGTGCTTGCCAGCCt                                   >  8:316821/1‑149 (MQ=255)
                                                                                                                    cGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCt                                   >  3:124506/1‑147 (MQ=255)
                                                                                                                     gACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCa                                 <  8:249936/148‑1 (MQ=255)
                                                                                                                     gACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAg                                >  2:329451/1‑149 (MQ=255)
                                                                                                                      aCCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGTCTCAGc                               <  2:76059/149‑1 (MQ=255)
                                                                                                                      aCCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGc                               <  6:304434/149‑1 (MQ=255)
                                                                                                                           gtgCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCcaa                            <  3:349714/147‑1 (MQ=255)
                                                                                                                           gtgCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGCTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCcaac                           <  1:309918/148‑1 (MQ=255)
                                                                                                                             gCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCaacaac                        >  7:64930/1‑149 (MQ=255)
                                                                                                                              ccGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCaacaac                        <  1:329451/148‑1 (MQ=255)
                                                                                                                              ccGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACt                       <  5:209408/149‑1 (MQ=255)
                                                                                                                               cGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTg                      <  5:18283/149‑1 (MQ=255)
                                                                                                                                      gcgcCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTctct               >  7:324353/1‑149 (MQ=255)
                                                                                                                                           gcGATGATGTCGATATCGAAATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCa                                       <  3:156336/119‑1 (MQ=255)
                                                                                                                                           gcGATGATGTCGATATCGAAATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCa                                       >  4:156336/2‑120 (MQ=255)
                                                                                                                                            cGATGATTTCGATATCGACATCATCGCACAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGAGGCGAACACTATTGGAATGTCCTTGACAGCCTCAGCCAACAACTGCTTCTCTTCGaag         >  6:343634/1‑147 (MQ=255)
                                                                                                                                             gatgatTTCGATATCGACATCATCGCCCAGCACCCGCGCAGACATATACAGCAGCTTGAGGCTAAGGTTGACGCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGACAACAACTGCTTCTCTTCGACt         >  7:292340/1‑148 (MQ=255)
                                                                                                                                              atgatTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGt       <  5:343634/149‑1 (MQ=255)
                                                                                                                                                   ttCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGAATGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAAc  >  5:7501/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGATCACGCGCGCCGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACGCTGAAGTTGGCGGCGAACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAAC  >  NC_002947/5373703‑5373998

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: ‑