Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 206,064 T→A 7.0% L3857Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947206,0640TA 7.0% 192.6 / 3.1 66L3857Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base A (3/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGG  >  NC_002947/205920‑206211
                                                                                                                                                |                                                                                                                                                   
gCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                 <  7:282890/149‑1 (MQ=39)
  aCGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                               <  4:157650/149‑1 (MQ=39)
    gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                             >  1:327666/1‑149 (MQ=39)
       ggTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGAccc                                                                                                                                          >  5:10444/1‑149 (MQ=39)
        gTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                         >  4:275347/1‑149 (MQ=39)
           aCCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCAGGTCGTGAgccagc                                                                                                                                       <  7:57965/148‑4 (MQ=255)
           aCCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGAGCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                       <  3:212776/148‑1 (MQ=37)
            ccGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGc                                                                                                                                     <  2:327666/149‑1 (MQ=39)
               acgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                               >  5:263686/1‑136 (MQ=40)
               acgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                               <  6:263686/136‑1 (MQ=40)
                                   aaCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATc                                                                                                                >  1:122704/1‑147 (MQ=39)
                                   aaCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCaa                                                                                                              <  4:100145/149‑1 (MQ=39)
                                    aCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                             <  6:10444/149‑1 (MQ=39)
                                     cGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAgg                                                                                                            >  7:182888/1‑149 (MQ=39)
                                      gAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                           >  4:144282/1‑149 (MQ=39)
                                          gCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGc                                                                                                       <  3:275347/149‑1 (MQ=39)
                                          gCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGc                                                                                                       >  7:287212/1‑149 (MQ=39)
                                                   tCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgc                                                                                              >  8:37965/1‑149 (MQ=255)
                                                      aGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACAGATATTCACCGTGCACGTAAACCAGCCGCTGGCTCACGAACAGGTCGAGACACAGAGCAACAACGCCCAGGTCCCCATCACGGCCTGCGAGACCAgcgcg                                                                                           >  8:259014/1‑149 (MQ=14)
                                                        gTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcgcgc                                                                                          >  7:39737/1‑148 (MQ=255)
                                                        gTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGcc                                                                                         <  8:182888/149‑1 (MQ=255)
                                                             cATCACAGCCGACCAGACTTCGGTGGCCGAGACCGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATaca                                                                                    >  2:56518/1‑149 (MQ=255)
                                                                 acaGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCTAGACCAGCGCGCCATACACCCa                                                                                <  5:63084/149‑1 (MQ=255)
                                                                   aGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACg                                                                              >  4:7610/1‑149 (MQ=255)
                                                                         ccAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACTGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTgcgc                                                                        <  1:275111/149‑1 (MQ=37)
                                                                          cAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATacac                                                                                   >  6:89683/1‑137 (MQ=255)
                                                                          cAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATacac                                                                                   <  5:89683/137‑1 (MQ=255)
                                                                           aGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTgcgc                                                                        >  7:178421/1‑147 (MQ=255)
                                                                            gACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                       <  6:301006/83‑1 (MQ=255)
                                                                            gACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                       >  5:301006/1‑83 (MQ=255)
                                                                            gACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAg                                                                     >  8:67078/1‑149 (MQ=255)
                                                                               ttCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                   >  4:18756/1‑148 (MQ=255)
                                                                                  ggTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                <  3:7610/148‑1 (MQ=255)
                                                                                    tGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCAAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCAGCGCAAGGCGAAGAc                                                             >  5:195642/1‑149 (MQ=255)
                                                                                       ccGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTTGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                          <  7:67078/149‑1 (MQ=255)
                                                                                       ccGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                          >  3:316042/1‑149 (MQ=255)
                                                                                        cGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACAccc                                                                                 >  5:105724/1‑125 (MQ=255)
                                                                                        cGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACAccc                                                                                 <  6:105724/125‑1 (MQ=255)
                                                                                               gTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCAGCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacga                                                  <  7:211728/149‑1 (MQ=255)
                                                                                                   aaCCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCt                                              >  1:213057/1‑149 (MQ=40)
                                                                                                      cGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGc                                           >  2:7683/1‑149 (MQ=39)
                                                                                                       gATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGc                                           >  1:90173/1‑148 (MQ=39)
                                                                                                             cACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGc                                           <  3:18756/142‑1 (MQ=37)
                                                                                                               ccGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGTCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGc                                  <  1:56518/149‑1 (MQ=37)
                                                                                                                cGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGcc                                 >  6:30641/1‑149 (MQ=37)
                                                                                                                    cACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTggg                             >  4:19121/1‑149 (MQ=37)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcgcg                                                                                           >  6:251656/1‑86 (MQ=38)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcgcg                                                                                           <  5:251656/86‑1 (MQ=38)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                        >  8:256136/1‑121 (MQ=38)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                        <  7:256136/121‑1 (MQ=38)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTggg                             >  3:107002/1‑148 (MQ=37)
                                                                                                                     aCGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGc                            <  3:19121/149‑1 (MQ=37)
                                                                                                                       gTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCAt                          <  2:90173/149‑1 (MQ=37)
                                                                                                                       gTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCAt                          >  2:256832/1‑149 (MQ=37)
                                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGCTGACGTCTACAACGCCGCTGGCGAGATCAGCCTGGGCATCAc                       >  8:128986/1‑149 (MQ=255)
                                                                                                                             cAGCGGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGTCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGAGATGACGTCTAGAACGACGCTTGCCAGATCAGCCTGGGCATCACGTc                    <  7:259014/149‑1 (MQ=11)
                                                                                                                                 cGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTctgct                <  8:287212/149‑1 (MQ=37)
                                                                                                                                 cGCTGGCTCACGACCTGGTCGTGACCCTGACCAACACCGCCCAGCTCACCATCCAGGCCCGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTctgct                <  8:206962/149‑1 (MQ=11)
                                                                                                                                  tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGa                                       >  3:31970/7‑125 (MQ=255)
                                                                                                                                  tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGa                                       <  4:31970/119‑1 (MQ=255)
                                                                                                                                      ggtCACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGCCCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGa           >  1:280904/3‑149 (MQ=255)
                                                                                                                                        tCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATg         >  7:32290/1‑149 (MQ=35)
                                                                                                                                        tCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATg         >  3:305242/1‑149 (MQ=35)
                                                                                                                                            gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGa     <  8:39737/149‑1 (MQ=34)
                                                                                                                                            gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGa     <  1:256832/149‑1 (MQ=34)
                                                                                                                                               ttGGTCGTGACCCTTAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATg         <  6:57553/141‑1 (MQ=21)
                                                                                                                                               cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGatcc  >  5:57553/1‑146 (MQ=21)
                                                                                                                                                |                                                                                                                                                   
GCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGG  >  NC_002947/205920‑206211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: