Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,008,318 C→T 7.9% S259S (TCG→TCA)  PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,008,3180CT 7.9% 181.1 / 4.0 63S259S (TCG→TCA) PP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATG  >  NC_002947/3008172‑3008462
                                                                                                                                                  |                                                                                                                                                
tGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATcgacga                                                                                                                                                >  1:346268/1‑149 (MQ=255)
       ccAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCa                                                                                                                                         >  2:234361/1‑149 (MQ=255)
        cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAt                                                                                                                                        >  8:155303/1‑149 (MQ=255)
        cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAt                                                                                                                                        >  7:55764/1‑149 (MQ=255)
         aGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                       <  5:136582/149‑1 (MQ=255)
          ggCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTg                                                                                                                                      <  8:45537/149‑1 (MQ=255)
            cgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                    <  2:92775/149‑1 (MQ=255)
            cgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                    <  3:215648/149‑1 (MQ=255)
              cATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCt                                                                                                                                  <  7:204209/149‑1 (MQ=255)
                    gCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTc                                                                                                                              <  4:342239/147‑1 (MQ=255)
                      gAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                       <  5:209324/136‑1 (MQ=255)
                      gAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                       >  6:209324/1‑136 (MQ=255)
                      gAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGc                                                                                                                          <  7:155303/149‑1 (MQ=255)
                           ccGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGGGCGCTTCGCGGGTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAc                                                                                                                     <  2:153436/149‑1 (MQ=255)
                                     tGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTg                                                                                                                                 >  3:56737/1‑127 (MQ=255)
                                     tGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTg                                                                                                                                 <  4:56737/127‑1 (MQ=255)
                                       aaaTCTTCCATGGCTGCTCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGGGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGcc                                                                                                          <  5:307178/148‑1 (MQ=255)
                                       aaaTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCg                                                                                                         >  1:154737/1‑149 (MQ=255)
                                         aTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCggg                                                                                                       <  4:71671/149‑1 (MQ=255)
                                          tCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACAGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgggg                                                                                                      >  5:319814/1‑149 (MQ=255)
                                           cTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCggggg                                                                                                     >  8:126665/1‑149 (MQ=255)
                                               cATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACg                                                                                                  >  5:211949/1‑148 (MQ=255)
                                                aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGa                                                                                                 >  4:311788/1‑148 (MQ=255)
                                                  ggCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                               <  3:78065/148‑1 (MQ=255)
                                                  ggCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                               >  4:78065/1‑148 (MQ=255)
                                                  ggCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                               >  7:70275/1‑148 (MQ=255)
                                                         cgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCAt                                                                                       <  1:234361/149‑1 (MQ=255)
                                                          gcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATg                                                                                      <  8:70275/149‑1 (MQ=255)
                                                            gcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGc                                                                                    <  3:311788/149‑1 (MQ=255)
                                                               aaGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCaa                                                                                 <  8:55764/149‑1 (MQ=255)
                                                                       ttCGCGGTTGGAGTGCGAAAACAGTTTTATCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGc                                                                                    <  6:241625/138‑1 (MQ=255)
                                                                       ttCGCGGTTGGAGTGCGAAAACAGTTTTATCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGc                                                                                    >  5:241625/1‑138 (MQ=255)
                                                                                   gTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGc                                                                       <  3:157145/139‑1 (MQ=255)
                                                                                   gTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGc                                                                       >  4:157145/1‑139 (MQ=255)
                                                                                      cGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgg                                                                                                        >  6:350876/1‑103 (MQ=255)
                                                                                      cGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCaa                                                          <  6:319814/149‑1 (MQ=255)
                                                                                       gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAAt                                                         >  7:93047/1‑149 (MQ=255)
                                                                                       gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAAt                                                         >  3:39828/1‑149 (MQ=255)
                                                                                         aaaCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATg                                                        >  5:247077/1‑148 (MQ=255)
                                                                                          aaCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCa                                                               >  3:200214/1‑140 (MQ=255)
                                                                                          aaCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCa                                                               <  4:200214/140‑1 (MQ=255)
                                                                                               ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGgca                                                 <  6:211949/149‑1 (MQ=255)
                                                                                                     cGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCAt                                             >  2:36386/1‑147 (MQ=255)
                                                                                                           gATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGCCCCAGGCAATGAAGGGCAGCATCAGttt                                       >  4:80671/1‑147 (MQ=255)
                                                                                                              gTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAgg                                                             >  5:85535/1‑122 (MQ=255)
                                                                                                              gTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAgg                                                             <  6:85535/122‑1 (MQ=255)
                                                                                                                 ggTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACa                                                                         >  6:94436/1‑107 (MQ=255)
                                                                                                                 ggTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACa                                                                         <  5:94436/107‑1 (MQ=255)
                                                                                                                  gtgtAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCCGGGCaaa                              >  6:43648/1‑149 (MQ=255)
                                                                                                                     tAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgcc                           <  2:154737/149‑1 (MQ=255)
                                                                                                                     tAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgcc                           >  8:330697/1‑149 (MQ=255)
                                                                                                                      aaGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgccg                          >  2:28121/1‑149 (MQ=255)
                                                                                                                            cAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGc                     >  8:277594/1‑148 (MQ=255)
                                                                                                                             aGAAACCCGCCCGGTCATCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGaa                                                      <  7:46237/114‑1 (MQ=255)
                                                                                                                             aGAAACCCGCCCGGTCATCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGaa                                                      >  8:46237/1‑114 (MQ=255)
                                                                                                                                aaCCCGCCCGGTCATCGACGCCCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAggg                <  8:93047/149‑1 (MQ=255)
                                                                                                                                aaCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAggg                <  1:36386/149‑1 (MQ=255)
                                                                                                                                 caatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCa                                     <  3:4181/123‑1 (MQ=255)
                                                                                                                                 caatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCa                                     >  4:4181/5‑127 (MQ=255)
                                                                                                                                 aCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                   <  5:114715/129‑1 (MQ=255)
                                                                                                                                 aCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                   >  6:114715/1‑129 (MQ=255)
                                                                                                                                    tgcccgGTCGTAGATGACCGCGCATTGGCCTGCTCCTGCCCGACCTTGCCAACCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGCCAAAGCCGCCGGCCAGGGCAAt            >  6:158601/2‑149 (MQ=255)
                                                                                                                                    cgcccgGTCATCGACGACCGGTCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAAt            <  5:265244/149‑1 (MQ=255)
                                                                                                                                          gTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATcgc      <  3:80671/149‑1 (MQ=255)
                                                                                                                                               cgacgaCCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATg  >  8:329287/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATG  >  NC_002947/3008172‑3008462

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: ‑