Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,373,861 C→A 9.1% V147F (GTC→TTC)  dapB ← 4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,373,8610CA 9.1% 198.3 / 9.6 66V147F (GTC→TTC) dapB4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base A (3/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAACCGGTCGTGCC  >  NC_002947/5373717‑5374008
                                                                                                                                                |                                                                                                                                                        
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                                                                                                                            <  1:73113/149‑1 (MQ=255)
                          gCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                                            <  3:106840/149‑1 (MQ=255)
                          gCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATc                                                                                                                            >  2:182348/1‑149 (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)
                                                           gCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑Ct                                                                                        <  4:87032/149‑1 (MQ=255)
                                                             acccACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTg                                                                                       >  4:233969/3‑148 (MQ=255)
                                                              aCCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTg                                                                                       >  3:178887/1‑147 (MQ=255)
                                                                 aCTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGtt                                                                                  >  2:197352/1‑149 (MQ=255)
                                                                   ttCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGGCCCCCAC‑‑G‑CTGAAGTTgg                                                                                >  4:131592/1‑149 (MQ=255)
                                                                      aCCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTggcgg                                                                             >  6:209408/1‑149 (MQ=255)
                                                                        ccATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCGAGCACCCGCGCCGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTggcggc                                                                            >  6:272851/1‑148 (MQ=255)
                                                                        ccATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTggcggc                                                                            >  5:304434/1‑148 (MQ=255)
                                                                          aagagcagtgcagtACCCGACGGCGCATCGACCTTGTGACGGTGGTGCGCCTCGATCATTTCGACAACCACATCATCGCCCAGCACCCGCGCAGCCATATCCCGCAGCTAGAGGCGAAGCTTGACCCCCAC‑‑C‑CTGCAGTTGGCGGC‑G‑aa                                                                       >  7:355900/5‑149 (MQ=255)
                                                                                 gtgcagtACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGAtt                                                                <  3:233969/149‑1 (MQ=255)
                                                                                      gtACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAAt                                                           <  1:197352/149‑1 (MQ=255)
                                                                                          ccGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAAt                                                           >  5:311577/1‑145 (MQ=255)
                                                                                          ccGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAAt                                                           <  6:311577/145‑1 (MQ=255)
                                                                                                gcgcATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCa                                                 <  1:182348/149‑1 (MQ=255)
                                                                                                 cgcATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACACGAGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAA‑‑G‑AAGAAGTAGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAg                                                >  8:205431/1‑149 (MQ=255)
                                                                                                    aTCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGGGCAGCCATATCCAGCAGCTAGAGGCTAAGGTTGACACCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑ACCACGATTGGAATGTGCTTGCCAGCCt                                             >  8:316821/1‑149 (MQ=255)
                                                                                                      cGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCt                                             >  3:124506/1‑147 (MQ=255)
                                                                                                       gACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCa                                           <  8:249936/148‑1 (MQ=255)
                                                                                                       gACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAg                                          >  2:329451/1‑149 (MQ=255)
                                                                                                        aCCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGTCTCAGc                                         <  2:76059/149‑1 (MQ=255)
                                                                                                        aCCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGc                                         <  6:304434/149‑1 (MQ=255)
                                                                                                             gtgCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCcaa                                      <  3:349714/147‑1 (MQ=255)
                                                                                                             gtgCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGCTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCcaac                                     <  1:309918/148‑1 (MQ=255)
                                                                                                               gCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCaacaac                                  >  7:64930/1‑149 (MQ=255)
                                                                                                                ccGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCaacaac                                  <  1:329451/148‑1 (MQ=255)
                                                                                                                ccGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACt                                 <  5:209408/149‑1 (MQ=255)
                                                                                                                 cGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTg                                <  5:18283/149‑1 (MQ=255)
                                                                                                                        gcgcCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTctct                         >  7:324353/1‑149 (MQ=255)
                                                                                                                             gcGATGATGTCGATATCGAAATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCa                                                 <  3:156336/119‑1 (MQ=255)
                                                                                                                             gcGATGATGTCGATATCGAAATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCa                                                 >  4:156336/2‑120 (MQ=255)
                                                                                                                              cGATGATTTCGATATCGACATCATCGCACAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGAGGC‑G‑AACACTATTGGAATGTCCTTGACAGCCTCAGCCAACAACTGCTTCTCTTCGaag                   >  6:343634/1‑147 (MQ=255)
                                                                                                                               gatgatTTCGATATCGACATCATCGCCCAGCACCCGCGCAGACATATACAGCAGCTTGAGGCTAAGGTTGACGCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGACAACAACTGCTTCTCTTCGACt                   >  7:292340/1‑148 (MQ=255)
                                                                                                                                atgatTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGt                 <  5:343634/149‑1 (MQ=255)
                                                                                                                                     ttCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGAATGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAAc            >  5:7501/1‑149 (MQ=255)
                                                                                                                                      tCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCACTGGAGTGAAGTTGGCGGCAGAAACACGATTGGTATGTACTTGCCAGCCTCAGCCAACAACTGCTTCTAATCGACTGt                 >  1:303016/1‑148 (MQ=255)
                                                                                                                                         atatCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAACCGGt        >  7:7829/1‑149 (MQ=255)
                                                                                                                                            tCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAACCGGTCGt     <  1:136238/149‑1 (MQ=255)
                                                                                                                                              gACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAACCGGTCGTGc   >  3:192044/1‑149 (MQ=255)
                                                                                                                                               aCATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAACTTGGCGGC‑G‑AACACCATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCATCGACTGTAAAAACGCTCGTGcc  >  7:292262/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                        
CGGTTTGCCCTTCACGGCCATACACGGCCACTTCCTGCAGATCACGGCCCAGCGCATTGGCAACCACTTCACCCATGCGCAGTGCAGTACCCGACGGCGCATCGACCTTGTGCCGGTGGTGCGCCTCGATGATTTCGATATCGACATCATCGCCCAGCACCCGCGCAGCCATATCCAGCAGCTTGAGGCTAAGGTTGACCCCCAC‑‑G‑CTGAAGTTGGCGGC‑G‑AACACGATTGGAATGTCCTTGCCAGCCTCAGCCAACAACTGCTTCTCTTCGACTGTAAAACCGGTCGTGCC  >  NC_002947/5373717‑5374008

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