Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 501,673 A→T 8.8% intergenic (+22/‑39) PP_0412 → / → PP_0413 polyamine ABC transportersubstrate‑binding protein/polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947501,6730AT 8.8% 144.9 / 3.3 56intergenic (+22/‑39)PP_0412/PP_0413polyamine ABC transportersubstrate‑binding protein/polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base T (3/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGCGCCGCT  >  NC_002947/501525‑501821
                                                                                                                                                     |                                                                                                                                                     
aGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCa                                                                                                                                                       <  4:176824/149‑1 (MQ=255)
      gaaagaCATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCgg                                                                                                                                                 >  1:87719/1‑149 (MQ=255)
       aaagaCATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGc                                                                                                                                                >  3:163694/1‑149 (MQ=255)
       aaagaCATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGc                                                                                                                                                <  2:95254/149‑1 (MQ=255)
        aagaCATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTTGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGcc                                                                                                                                               <  1:33795/149‑1 (MQ=255)
         agaCATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCTTTGCGGccc                                                                                                                                              <  5:97619/149‑1 (MQ=255)
            cATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCTTTGCGGccccac                                                                                                                                           <  7:171949/149‑3 (MQ=255)
             aTGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCaa                                                                                                                                         <  1:240285/149‑1 (MQ=255)
                     cACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCt                                                                                                                                 <  8:266471/149‑1 (MQ=255)
                      aCGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTg                                                                                                                                >  1:75784/1‑149 (MQ=255)
                          cTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTAtt                                                                                                                            >  8:275635/1‑149 (MQ=255)
                             agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGc                                                                                                                          <  4:145534/148‑1 (MQ=255)
                             agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCg                                                                                                                         <  6:236038/149‑1 (MQ=255)
                             agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCg                                                                                                                         <  7:161356/149‑1 (MQ=255)
                             agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCg                                                                                                                         <  4:235867/149‑1 (MQ=255)
                              gaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCgg                                                                                                                        <  2:290934/149‑1 (MQ=255)
                                aCATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGa                                                                                                                       <  4:163694/148‑1 (MQ=255)
                                  aTTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGACGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAACCGGCCACCAATTAAGCCAGTATTGCGGAGtt                                                                                                                    >  7:191531/1‑149 (MQ=255)
                                      ccAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGGCCGCTTTGCGGccccaca                                                                                                                                          >  6:235775/1‑122 (MQ=255)
                                        aaCCAGGTGCAGAGCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCAt                                                                                                              <  6:228241/149‑1 (MQ=255)
                                          ccAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATgg                                                                                                            <  2:87719/149‑1 (MQ=255)
                                                 cAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTg                                                                                                      <  7:115202/148‑1 (MQ=255)
                                                    aTCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTAAGCCTGTTTTGGGGAGTTAGCCATGGCCATTGCAGt                                                                                                  >  3:317569/1‑149 (MQ=255)
                                                          gTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCt                                                                                            <  6:277714/149‑1 (MQ=255)
                                                              cGTTCTGGGCCGACAACAGCGAGCAGGTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAAc                                                                                        <  1:196780/149‑1 (MQ=255)
                                                               gTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAAc                                                                                        >  8:279033/1‑148 (MQ=255)
                                                                   tGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAgg                                                                                   >  3:88994/1‑149 (MQ=255)
                                                                     ggCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAg                                                                                    <  4:143938/146‑1 (MQ=255)
                                                                        cGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAgg                                                                                   <  7:94193/144‑1 (MQ=255)
                                                                        cGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAgg                                                                                   >  8:94193/1‑144 (MQ=255)
                                                                         gACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCACAAAGCGGCCCCCGAAATCACCCAGTATTGCGGAGTTAGCCAATGCCAATGCAGTGCCCCTCAACAAAGGCGCAg                                                                              >  1:102265/1‑149 (MQ=255)
                                                                           caacaGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTc                                                                           >  4:225157/1‑149 (MQ=255)
                                                                                  gAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGtc                                                                     >  3:117376/1‑148 (MQ=255)
                                                                                    gcagcTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGtctc                                                                   >  7:224909/1‑148 (MQ=255)
                                                                                     cagcTGGAGCAACGCTTCAACGCCTGGGTGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCaa                                                                 <  5:54189/149‑1 (MQ=255)
                                                                                       gcTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCCGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAACGCGCAGGTCCAAGTCTCAag                                                                >  1:289049/1‑148 (MQ=255)
                                                                                          ggAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTAAGCCTGTATTGCGGAGTTAGCCATGGCCATTTGACTGCCCCTAAAAGAAGGCCCAGGTGCAAGTCTCAagcagc                                                            >  4:304593/1‑149 (MQ=255)
                                                                                             gCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAAATCAGCCTGTATTGCGGAGTTCGCCATGGCAATAGCACTGCCCCTCAACGAAGGAGCAGGTCCAAGTCTCAATCAGCGcc                                                         >  1:5826/1‑149 (MQ=255)
                                                                                               aaCGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGCGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCt                                                        >  3:276885/1‑148 (MQ=255)
                                                                                                  gCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCaaa                                                    <  3:304593/149‑1 (MQ=255)
                                                                                                       aaCGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGCGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGc                                               <  4:276885/149‑1 (MQ=255)
                                                                                                       aaCGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTAAGCCTGTATTGCGGAGTTGGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTGTCAAGCAGCGCCTCAAACACGc                                               >  2:224488/1‑149 (MQ=255)
                                                                                                            cTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCCCAAATCGGCCCCCAAATGAGCCAGTATTGCGGAGTTCGCCATGGCCATTTCAGTGCCCCAAAAAGAATGCGCAGGGGCAAGACTCAAGCACCGCCGCCCACACGCCGa                                            >  2:285068/1‑147 (MQ=255)
                                                                                                                        gaagTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCt                              <  8:224909/149‑1 (MQ=255)
                                                                                                                        gaagTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCt                              <  2:102265/149‑1 (MQ=255)
                                                                                                                        gaagTAAACGTCCAGGGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCt                              <  3:12417/149‑1 (MQ=255)
                                                                                                                             aaaCGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGaa                          <  2:75784/148‑1 (MQ=255)
                                                                                                                              aaCGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAg                         <  8:160292/148‑1 (MQ=255)
                                                                                                                                    cAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAg                   >  6:326977/1‑148 (MQ=255)
                                                                                                                                    cAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAg                   <  5:326977/148‑1 (MQ=255)
                                                                                                                                    cAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAg                   <  4:317569/148‑1 (MQ=255)
                                                                                                                                      tgaattgggggCCGCTTTGCGGCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCg                >  3:44806/8‑149 (MQ=255)
                                                                                                                                           tgGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGa            <  4:117376/148‑1 (MQ=255)
                                                                                                                                             ggggCCGCTTTGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATc          >  6:121768/1‑149 (MQ=255)
                                                                                                                                              gggCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATc          <  5:252882/148‑1 (MQ=255)
                                                                                                                                              gggCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCAGAGCGGGTCAACCGCTGGAAGGCGCCGGCGTTGAtag         >  6:252882/1‑147 (MQ=255)
                                                                                                                                                   gCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGcgccg    >  6:62969/1‑149 (MQ=255)
                                                                                                                                                     aaaGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCTCAGGCGTTGATCGCGCCGCt  >  8:70571/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGCGCCGCT  >  NC_002947/501525‑501821

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: