Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,623,202 C→T 8.4% intergenic (+38/‑208) PP_5517 → / → PP_2295 hypothetical protein/antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,623,2020CT 8.4% 236.5 / 8.9 83intergenic (+38/‑208)PP_5517/PP_2295hypothetical protein/antirestriction protein
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (4/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAC  >  NC_002947/2623058‑2623345
                                                                                                                                                |                                                                                                                                                
cGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCggg                                                                                                                                               <  8:211085/148‑1 (MQ=255)
cGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCggg                                                                                                                                               <  7:282796/148‑1 (MQ=255)
    gATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGcccc                                                                                                                                           <  1:315734/148‑2 (MQ=255)
     aTGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCttt                                                                                                                                         >  8:18569/1‑149 (MQ=255)
      tGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGCCCGATTGGGCCCcttga                                                                                                                                       <  2:17258/149‑3 (MQ=255)
            tgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGt                                                                                                                                  <  1:17060/149‑1 (MQ=255)
             ggAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGt                                                                                                                                  <  4:124097/148‑1 (MQ=255)
             ggAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGt                                                                                                                                  >  3:124097/1‑148 (MQ=255)
              gAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCt                                                                                                                                >  3:205394/1‑149 (MQ=255)
                ggTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGt                                                                                                                              <  4:205394/149‑1 (MQ=255)
                    ttCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGt                                                                                                                              <  7:262135/145‑1 (MQ=255)
                       cAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGa                                                                                                                       >  3:361312/1‑149 (MQ=255)
                         aCAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAtt                                                                                                                     <  5:42088/149‑1 (MQ=255)
                           aGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCt                                                                                                                   >  7:5586/1‑149 (MQ=255)
                            gAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCtttt                                                                                                                                        <  8:179890/127‑1 (MQ=255)
                            gAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCtttt                                                                                                                                        >  7:179890/1‑127 (MQ=255)
                            gAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTg                                                                                                                  <  6:306564/149‑1 (MQ=255)
                                      gcgGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCg                                                                                                         >  6:324379/1‑148 (MQ=255)
                                        ggTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGGCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAg                                                                                                      >  2:202307/1‑149 (MQ=255)
                                            aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTc                                                                                                                                      <  1:333909/113‑1 (MQ=255)
                                            aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTc                                                                                                                                      >  2:333909/1‑113 (MQ=255)
                                            aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAgggg                                                                                                   <  2:3362/148‑1 (MQ=255)
                                            aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGc                                                                                                  <  7:215454/149‑1 (MQ=255)
                                            aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGc                                                                                                  >  3:109446/1‑149 (MQ=255)
                                             ccttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCt                                                                                                 <  3:355968/149‑1 (MQ=255)
                                              cttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCt                                                                                                 >  5:117686/1‑148 (MQ=255)
                                                     tctcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGccc                                                                                           <  7:18569/147‑1 (MQ=255)
                                                       tcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                       <  1:42085/149‑1 (MQ=255)
                                                         cGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                       <  2:133530/147‑1 (MQ=255)
                                                               gTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                               <  4:131570/149‑1 (MQ=255)
                                                                tCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgcg                                                                              >  1:41814/1‑149 (MQ=255)
                                                                   aGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAg                                                                           <  6:63534/149‑1 (MQ=255)
                                                                    gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGga                                                                                                                >  6:195829/1‑111 (MQ=255)
                                                                    gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGga                                                                                                                <  5:195829/111‑1 (MQ=255)
                                                                    gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAg                                                                           >  2:80819/1‑148 (MQ=255)
                                                                     ccTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGaa                                                                         <  5:246229/149‑1 (MQ=255)
                                                                        tCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgcg                                                                      >  2:276551/1‑149 (MQ=255)
                                                                           acacctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgcgcgc                                                                   >  7:39578/1‑149 (MQ=255)
                                                                             acctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGcc                                                                  >  5:134707/1‑148 (MQ=255)
                                                                               ctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCaa                                                                >  8:65226/1‑148 (MQ=255)
                                                                                   accAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAAtt                                                            >  8:313346/1‑148 (MQ=255)
                                                                                   accAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                           >  1:327104/1‑149 (MQ=255)
                                                                                    ctaaGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAc                                                          <  5:324379/147‑1 (MQ=255)
                                                                                    ccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAc                                                          <  7:65226/149‑1 (MQ=255)
                                                                                     cAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAATCGCGCCAAAATTAc                                                          >  6:323486/1‑148 (MQ=255)
                                                                                     cAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAc                                                          >  4:209450/1‑148 (MQ=255)
                                                                                               aGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGt                                               <  4:109446/148‑1 (MQ=255)
                                                                                                 gAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCAt                                            >  5:161525/1‑149 (MQ=255)
                                                                                                    gcgcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATcc                                                                                 >  5:281749/1‑110 (MQ=255)
                                                                                                    gcgcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATcc                                                                                 <  6:281749/110‑1 (MQ=255)
                                                                                                    gcgcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTc                                         <  4:361312/149‑1 (MQ=255)
                                                                                                    gcgcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTc                                         >  2:316284/1‑149 (MQ=255)
                                                                                                         gACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCcaca                                     >  3:54848/1‑148 (MQ=255)
                                                                                                          aCGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACaa                                    <  7:313346/148‑1 (MQ=255)
                                                                                                              aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                       <  4:29920/94‑1 (MQ=255)
                                                                                                              aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                       >  3:29920/1‑94 (MQ=255)
                                                                                                              aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGa                               <  6:134707/149‑1 (MQ=255)
                                                                                                                 tCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCcac                                      <  8:132393/139‑1 (MQ=255)
                                                                                                                 tCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCcac                                      >  7:132393/1‑139 (MQ=255)
                                                                                                                         aaTCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATgg                                >  1:74228/1‑137 (MQ=255)
                                                                                                                         aaTCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATgg                                <  2:74228/137‑1 (MQ=255)
                                                                                                                          aTCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACtt                   <  2:327104/149‑1 (MQ=255)
                                                                                                                           tCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTg                  >  1:195553/1‑149 (MQ=255)
                                                                                                                               acCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATg              >  6:293372/1‑149 (MQ=255)
                                                                                                                               acCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTACACAAATGGAAATAAATAAATTGCAtt              >  3:229948/1‑148 (MQ=255)
                                                                                                                                ccTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATgg             >  5:323456/1‑149 (MQ=255)
                                                                                                                                 gaaaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  6:364373/7‑149 (MQ=255)
                                                                                                                                 gaaaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  3:36731/7‑149 (MQ=255)
                                                                                                                                 gaaaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  7:4868/7‑149 (MQ=255)
                                                                                                                                 cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGTGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            <  3:209450/149‑1 (MQ=255)
                                                                                                                                 cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  1:79807/1‑149 (MQ=255)
                                                                                                                                 cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg            >  1:342750/1‑149 (MQ=255)
                                                                                                                                 cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCAATTCCACAAATGGAAATAAATCCCTTCCAAggg            >  8:192947/1‑149 (MQ=255)
                                                                                                                                  tCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATgggg           >  1:238726/1‑149 (MQ=255)
                                                                                                                                    aaGGGTCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGTAATAAATCACTTGCATGGGgtg         <  5:323486/149‑1 (MQ=255)
                                                                                                                                    aaGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACTCGTCATTTCCACAAATGGAAATATATCACTGGCATGGtgtg         >  6:122311/1‑149 (MQ=255)
                                                                                                                                    aaGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATgggg           <  8:39578/147‑1 (MQ=255)
                                                                                                                                       gggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCg      >  6:47090/1‑149 (MQ=255)
                                                                                                                                         gCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCATGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGACAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGtt    <  6:323456/149‑1 (MQ=255)
                                                                                                                                         gCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGTGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAAGAAAACTACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGt     <  1:80819/149‑1 (MQ=255)
                                                                                                                                          cccGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTa   <  5:364373/149‑1 (MQ=255)
                                                                                                                                          cccAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTa   >  5:192184/1‑149 (MQ=255)
                                                                                                                                           ccAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAc  <  4:54848/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAC  >  NC_002947/2623058‑2623345

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: