Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 501,674 A→T 8.9% intergenic (+23/‑38) 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,6740AT 8.9% 153.4 / 5.7 56intergenic (+23/‑38)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 = 9.16e-01

GAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGCGCCGCTG  >  NC_002947/501531‑501822
                                                                                                                                                |                                                                                                                                                     
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                                                                                                                          <  4:235867/149‑1 (MQ=255)
                       agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCg                                                                                                                          <  6:236038/149‑1 (MQ=255)
                       agaACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCg                                                                                                                          <  7:161356/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)
                                                                                                                                                aaGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGCGCCGCTg  >  8:287549/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGTGCAGATCGACGTGGCGTTCTGGGCCGACAACAGCGAGCAGCTGGAGCAACGCTTCAACGCCTGGGCGGCGAAGAAGTAAACGTCCAGTGTGTGGGGCCGCAAAGCGGCCCCCAATTCAGCCTGTATTGCGGAGTTCGCCATGGCCATTGCAGTGCCCCTCAACGAAGGCGCAGGTCCAAGTCTCAAGCAGCGCCTCAAACACGCCGAGCGGGTCAACCGCTGGAAGGCGCAGGCGTTGATCGCGCCGCTG  >  NC_002947/501531‑501822

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: