Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,071,883 C→G 6.0% L214F (TTG→TTC pstB‑II ← phosphate ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,071,8830CG 6.0% 274.5 / 6.1 83L214F (TTG→TTCpstB‑IIphosphate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (42/36);  new base G (3/2);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTC  >  NC_002947/6071741‑6072029
                                                                                                                                              |                                                                                                                                                  
tatcggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTc                                                                                                                                              <  4:116429/149‑1 (MQ=255)
tatcggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTc                                                                                                                                              <  4:82798/149‑1 (MQ=255)
      tatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCgaagaa                                                                                                                                        <  3:66233/149‑5 (MQ=255)
        tcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg                                                                                                                                      <  6:254618/149‑1 (MQ=255)
        atcccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg                                                                                                                                      <  5:202883/147‑1 (MQ=255)
         ccccGAATTCGACCAGTTTGCCCATGTACATGTACGGGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCggg                                                                                                                                     <  3:267162/149‑1 (MQ=255)
          cccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGt                                                                                                                                    >  7:313076/1‑149 (MQ=255)
          cccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGt                                                                                                                                    >  6:338123/1‑149 (MQ=255)
           ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTc                                                                                                                                   >  8:49294/1‑149 (MQ=255)
             gAATTCGACCAGTTTGCCCATGTACATGAACGCGTTGTAATCCGAAACACGCGCCTCTTGCTTCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCa                                                                                                                                 <  5:83582/149‑1 (MQ=255)
             gAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg                                                                                                                                      >  6:202883/1‑144 (MQ=255)
                     ccAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAg                                                                                                                         <  8:59229/149‑1 (MQ=255)
                              ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGa                                                                                                                          <  4:319668/139‑1 (MQ=255)
                              ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGa                                                                                                                          >  3:319668/1‑139 (MQ=255)
                               cATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                               <  6:365088/149‑1 (MQ=255)
                                aTGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                               >  7:159158/1‑148 (MQ=255)
                                aTGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTc                                                                                                              <  8:311122/149‑1 (MQ=255)
                                 tGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTc                                                                                                              >  2:210999/1‑148 (MQ=255)
                                   tACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAag                                                                                                           >  1:104899/1‑149 (MQ=255)
                                    aCATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAag                                                                                                           >  3:267174/1‑148 (MQ=255)
                                      aTGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                               <  5:76058/142‑1 (MQ=255)
                                      aTGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                               >  6:76058/1‑142 (MQ=255)
                                      aTGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcag                                                                                                        >  3:6508/1‑149 (MQ=255)
                                         aaCGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATc                                                                                                                                        <  8:91804/114‑1 (MQ=255)
                                         aaCGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATc                                                                                                                                        >  7:91804/1‑114 (MQ=255)
                                                 gtAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGtt                                                                                                                  >  2:218239/1‑128 (MQ=255)
                                                 gtAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGtt                                                                                                                  <  1:218239/128‑1 (MQ=255)
                                                     tCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTc                                                                                         >  2:297177/1‑149 (MQ=255)
                                                        gAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCATGGTCGAGATCGGGTGCCGCTGCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCAc                                                                                      >  4:17738/1‑149 (MQ=255)
                                                             acGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCaa                                                                                 >  8:127426/1‑149 (MQ=255)
                                                               gcgcCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATg                                                                               <  5:338123/149‑1 (MQ=255)
                                                               gcgcCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGAGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATg                                                                               >  1:287597/1‑149 (MQ=255)
                                                                     cTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCgg                                                                         <  8:313076/149‑1 (MQ=255)
                                                                      ttGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACgg                                                                                    <  2:11008/137‑1 (MQ=255)
                                                                      ttGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACgg                                                                                    >  1:11008/1‑137 (MQ=255)
                                                                        gctgctTGTTGTGGGTCACTATGACTATGGTGTACTTTGATTTGAATTCATAGATCAGTTCTTCGACCTTCAATGTCGTGATCGGGTCCATCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCTCGGCAATGGTGCGGGc                                                                       <  2:287597/148‑1 (MQ=255)
                                                                           gcATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGAt                                                                    >  7:290732/1‑148 (MQ=255)
                                                                                  gtgGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCTAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGa                                                             <  8:290732/148‑1 (MQ=255)
                                                                                     ggTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAcc                                                                <  8:242982/142‑1 (MQ=255)
                                                                                     ggTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAcc                                                                >  7:242982/1‑142 (MQ=255)
                                                                                     ggTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgc                                                          <  4:6508/148‑1 (MQ=255)
                                                                                           gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgctgctgct                                                   >  1:65065/1‑149 (MQ=255)
                                                                                           gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgctgctgct                                                   >  2:85085/1‑149 (MQ=255)
                                                                                                 gatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCAc                                             <  1:55281/149‑1 (MQ=255)
                                                                                                     gtgtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGa                                         <  7:127426/149‑1 (MQ=255)
                                                                                                      tgtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGGGATGACCAGACGCTGCTGCTGGCCACCGGAc                                        <  4:267174/149‑1 (MQ=255)
                                                                                                       gtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGa                                         >  7:126052/1‑147 (MQ=255)
                                                                                                       gtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGACGACCAGACGCTGCTGCTGGCCACCGGa                                         >  2:240944/1‑147 (MQ=255)
                                                                                                         aCTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGGGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcagcag                                                                                                     <  8:48557/85‑1 (MQ=255)
                                                                                                         aCTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGGGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcagcag                                                                                                     >  7:48557/1‑85 (MQ=255)
                                                                                                          cTTGGATTTCAATTCATAGATCAGTTCTTCGCCCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGc                                    <  6:371689/149‑1 (MQ=255)
                                                                                                          cTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAgg                                     <  2:104899/148‑1 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCTCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCAGACAGGcc                                   >  3:268724/1‑149 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGc                                                                           >  5:244676/1‑109 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGc                                                                           <  6:244676/109‑1 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc                                   >  7:292950/1‑149 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc                                   >  4:211225/1‑149 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc                                   >  3:102510/1‑149 (MQ=255)
                                                                                                           ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGACCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc                                   >  4:75947/1‑149 (MQ=255)
                                                                                                             ggATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGccc                                  <  1:85085/148‑1 (MQ=255)
                                                                                                              gATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGACGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCa                                 <  1:240944/148‑1 (MQ=255)
                                                                                                               aTTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAgc                               >  1:184223/1‑149 (MQ=255)
                                                                                                                tttCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAgcg                              >  7:314124/1‑149 (MQ=255)
                                                                                                                  tCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCTATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAgcgc                             <  2:65065/148‑1 (MQ=255)
                                                                                                                  tCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCCCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCAGACAGGCCCAgcgc                             <  4:268724/148‑1 (MQ=255)
                                                                                                                        cATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTc                       <  8:126052/148‑1 (MQ=255)
                                                                                                                        cATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCAGGCCACCGGACAGGCCCAGCGCCGACTc                       <  1:210999/148‑1 (MQ=255)
                                                                                                                               cAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGc               >  4:33946/1‑149 (MQ=255)
                                                                                                                                ccgatctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCACCGACTCGTGCAGGCg              >  6:138872/7‑149 (MQ=255)
                                                                                                                                 gTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCgg             >  1:104273/1‑149 (MQ=255)
                                                                                                                                  ttcttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGt            >  1:251804/1‑149 (MQ=255)
                                                                                                                                  gatctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGTAGGCGGt            >  6:68845/5‑149 (MQ=255)
                                                                                                                                   atctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTc           >  1:302833/4‑149 (MQ=255)
                                                                                                                                    cttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTcc          >  4:153704/1‑149 (MQ=255)
                                                                                                                                       cGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt        <  4:102510/148‑1 (MQ=255)
                                                                                                                                       cGACCTTCAAGGTCGAGATCGGATCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt        >  2:179323/1‑148 (MQ=255)
                                                                                                                                       cGACCTTCAAGGTCGAGACCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt        <  3:75947/148‑1 (MQ=255)
                                                                                                                                        gACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGa      >  5:58468/1‑149 (MQ=255)
                                                                                                                                         aCCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGTAGGCGGTCCTTGAc     <  5:68845/149‑1 (MQ=255)
                                                                                                                                         aCCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGAc     >  1:309711/1‑149 (MQ=255)
                                                                                                                                            ttCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGAt                                                                    <  2:336703/83‑1 (MQ=255)
                                                                                                                                            ttCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGAt                                                                    >  1:336703/1‑83 (MQ=255)
                                                                                                                                            ttCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTc  >  6:182509/1‑149 (MQ=255)
                                                                                                                                              cAAGGTCGAGATCGGGTCCAGCGCCGAGCACTGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTc  <  3:211225/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTC  >  NC_002947/6071741‑6072029

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: