Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 563,005 C→A 7.1% P30T (CCT→ACT) ‡ rpoA → DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947563,0050CA 7.1% 221.9 / 4.2 77P30T (CCT→ACT) ‡rpoADNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base A (4/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

ATCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGAC  >  NC_002947/562857‑563148
                                                                                                                                                    |                                                                                                                                                
aTCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                  >  4:145058/1‑149 (MQ=255)
   aaCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctc                                                                                                                                               >  8:119884/1‑149 (MQ=255)
    aCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctct                                                                                                                                              >  4:160529/1‑149 (MQ=255)
      gAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGATCGAGCCtctc                                                                                                                                             >  4:167042/1‑148 (MQ=255)
         aaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCAAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCACCGACACATTGAAGTGCAGGTAGACAGTCCGACCCGCGCCAAAATTACTCTCGATCATCTCGAGc                                                                                                                                         >  4:201129/1‑149 (MQ=255)
           ccTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGt                                                                                                                                       <  5:2989/149‑1 (MQ=255)
            cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                     <  1:108666/149‑1 (MQ=255)
            cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                     <  3:160529/149‑1 (MQ=255)
            cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                     >  2:5292/1‑149 (MQ=255)
               aTTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                  <  1:164287/149‑1 (MQ=255)
                   tCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCgg                                                                                                                              >  3:146146/1‑149 (MQ=255)
                      aGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                  <  3:78310/127‑1 (MQ=255)
                      aGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                  >  4:78310/1‑127 (MQ=255)
                        ctctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATa                                                                                                                         >  7:212138/1‑149 (MQ=255)
                          ctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAGATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAcc                                                                                                                       >  7:159331/1‑149 (MQ=255)
                           tACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAccc                                                                                                                      >  1:282991/1‑149 (MQ=255)
                             cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCGCATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                     <  7:77787/148‑1 (MQ=255)
                             cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTg                                                                                                                    <  8:212138/149‑1 (MQ=255)
                              tCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                   >  2:278525/1‑149 (MQ=255)
                              tCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                   >  2:250430/1‑149 (MQ=255)
                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                   <  3:145058/147‑1 (MQ=255)
                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGc                                                                                                                 <  2:99929/149‑1 (MQ=255)
                                  aGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCaa                                                                                                               <  6:191193/149‑1 (MQ=255)
                                      aGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAAcgcg                                                                                                           <  4:146146/149‑1 (MQ=255)
                                             gAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGTTGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGcgcc                                                                                                    <  1:250430/149‑1 (MQ=255)
                                             gAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGcgcc                                                                                                    >  3:11862/1‑149 (MQ=255)
                                             caaaaTAGTTGAATCCATGAAGAATTCGGTAAATCAGTACCTGACTACCCCCCACATTGAAGTGCAGGAAGTCAGTCCGACCCGCGCCAAAATAACGATCTAGCATCTCGACCGAGTATACGGCAATACCCTGCTCAACGAGAAGcgcc                                                                                                    >  8:223477/2‑149 (MQ=255)
                                                     gttcATCCATGCAGATTTCGTTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGCACCGCGCCAAAATTACGCTCGAGACTCTCGAGCGTGGTTTCGGCCAGACCCTGGGCAACGCGCTGCGCCGCATccag                                                                                            >  6:270962/4‑147 (MQ=255)
                                                     gTGCATCCATGCAGATTTCGGTACATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtg                                                                                            <  8:159331/149‑1 (MQ=255)
                                                       gCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAAAGCGCTGCGCCGCATCCtgtt                                                                                          >  1:99909/1‑149 (MQ=255)
                                                        catgcatgCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                          <  5:248655/148‑1 (MQ=255)
                                                        cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                          <  7:119884/148‑1 (MQ=255)
                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                        >  5:45161/1‑149 (MQ=255)
                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                        >  7:14690/1‑149 (MQ=255)
                                                          tCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtc                                                                                       <  4:41773/149‑1 (MQ=255)
                                                            cATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcct                                                                                     <  1:278525/149‑1 (MQ=255)
                                                                 aGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTTAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATg                                                                                >  7:235940/1‑149 (MQ=255)
                                                                     ttCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTTAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                             <  8:235940/148‑1 (MQ=255)
                                                                     ttCGGTAAATGAGTTCCTGACTCACCTCCACATTGATTTGCAGGTTGTCAGTCCGACCCGCGCCAAAATTACGCTGGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                             <  3:201129/148‑1 (MQ=255)
                                                                        ggTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGc                                                                          >  4:156532/1‑148 (MQ=255)
                                                                        ggTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGc                                                                          <  3:156532/148‑1 (MQ=255)
                                                                            aaTGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCAAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgtg                                                                      <  6:45161/148‑1 (MQ=255)
                                                                             aTGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCa                                                                    <  4:11862/149‑1 (MQ=255)
                                                                                aGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCagta                                                                 <  3:167042/149‑1 (MQ=255)
                                                                                 gTTCCTTACTCCACGCCACATTTCTGTGCAGGTAGTCAGTCCGACCCTCGCCAAAATTACGCTGGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCGTGGCTGTGCagtag                                                                <  7:223477/149‑1 (MQ=255)
                                                                                    ccTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                              >  6:150994/1‑148 (MQ=255)
                                                                                     cTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                              >  4:145142/1‑147 (MQ=255)
                                                                                            cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTg                                                                            <  4:244969/126‑1 (MQ=255)
                                                                                            cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTg                                                                            >  3:244969/1‑126 (MQ=255)
                                                                                                  acaTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGAt                                               >  1:132874/1‑149 (MQ=255)
                                                                                                   caTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATg                                              <  3:145142/149‑1 (MQ=255)
                                                                                                                  agtcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGAGTACTcatggag                                >  8:112497/1‑142 (MQ=255)
                                                                                                                    tcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTAc                             >  1:172292/1‑149 (MQ=255)
                                                                                                                      agtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTc                           >  4:141459/1‑149 (MQ=255)
                                                                                                                        tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcg                         <  5:150994/149‑1 (MQ=255)
                                                                                                                        tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGAGTACTCATGGAGTACGcca                         <  7:112497/149‑2 (MQ=255)
                                                                                                                          cGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcgc                        >  2:233291/1‑148 (MQ=255)
                                                                                                                                gcgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCg                   >  1:192727/1‑148 (MQ=255)
                                                                                                                                 cgcCAAAATTAGGCTCGAGCCTCTCGTGCGTTGTTTGGGCCATTCCCTGGGCAACGCGCTGCGCCGCATCCTGTTTTCCTCCATGCCTGTCTGTGCAGTAGTCGAGGCCTAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                <  5:270962/149‑1 (MQ=255)
                                                                                                                                 cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                <  5:259532/149‑1 (MQ=255)
                                                                                                                                 cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                <  4:6663/149‑1 (MQ=255)
                                                                                                                                 cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                <  7:279810/149‑1 (MQ=255)
                                                                                                                                 cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                <  2:192727/149‑1 (MQ=255)
                                                                                                                                   ccAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTAGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgt              >  2:113290/1‑149 (MQ=255)
                                                                                                                                    cgaaaccaCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCCTCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtg             >  7:207931/8‑149 (MQ=255)
                                                                                                                                    cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgt                                                                       >  4:49677/1‑91 (MQ=255)
                                                                                                                                    cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgt                                                                       <  3:49677/91‑1 (MQ=255)
                                                                                                                                    cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATc                                                  >  7:123583/1‑112 (MQ=255)
                                                                                                                                    cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATc                                                  <  8:123583/112‑1 (MQ=255)
                                                                                                                                    cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtg             <  3:141459/149‑1 (MQ=255)
                                                                                                                                     aaaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtgt            <  2:282991/149‑1 (MQ=255)
                                                                                                                                      aaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtgt            >  4:139607/1‑148 (MQ=255)
                                                                                                                                         ttACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg        >  6:204193/1‑149 (MQ=255)
                                                                                                                                          caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGc                                            >  2:307949/2‑112 (MQ=255)
                                                                                                                                          caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGc                                            <  1:307949/111‑1 (MQ=255)
                                                                                                                                           aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTGCATGAGTACTCCGCGATCGAAGGTGTTCAg        <  6:197816/147‑1 (MQ=255)
                                                                                                                                           aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCAGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg        >  5:197816/1‑147 (MQ=255)
                                                                                                                                             gCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGa                                 <  6:122846/120‑1 (MQ=255)
                                                                                                                                             gCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGa                                 >  5:122846/1‑120 (MQ=255)
                                                                                                                                                cGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTTGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTACGCGATCGAAGGGGTTCAGGAAGAc  >  8:269999/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGAC  >  NC_002947/562857‑563148

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: