Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 583,396 T→C 8.6% V421A (GTG→GCG) ‡ selB → selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947583,3960TC 8.6% 154.5 / 4.9 59V421A (GTG→GCG) ‡selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base C (2/3);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

GCAGGTGCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTG  >  NC_002947/583250‑583544
                                                                                                                                                   |                                                                                                                                                    
gCAGGTGCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGgtggtg                                                                                                                                                     <  4:183160/148‑1 (MQ=255)
    gTGATGCGCGATGCAGCGGATCTGGAGCAGGCCCTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTC‑CGGGTGGTGGCCAc                                                                                                                                                >  6:68997/1‑149 (MQ=255)
      gctgcGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCAccc                                                                                                                                              >  1:182189/1‑149 (MQ=255)
       ctgcGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCAccc                                                                                                                                              <  1:158662/148‑1 (MQ=255)
         gcgcGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGCCACCACCCgca                                                                                                                                           <  7:277500/149‑3 (MQ=255)
           gcgATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGCCACCACCCgcaca                                                                                                                                         <  7:4617/149‑5 (MQ=255)
            cgATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAgg                                                                                                                                        <  1:270760/149‑1 (MQ=255)
              aTGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGc                                                                                                                                      <  5:287440/149‑1 (MQ=255)
                  aGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAAc                                                                                                                                  <  6:267662/149‑1 (MQ=255)
                      gATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAActgct                                                                                                                              >  1:85401/1‑149 (MQ=255)
                           ggAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCg                                                                                                                         >  8:33337/1‑149 (MQ=255)
                           ggAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCg                                                                                                                         >  2:97709/1‑149 (MQ=255)
                            gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACTTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  8:60288/147‑1 (MQ=255)
                            gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  2:182189/147‑1 (MQ=255)
                            gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  7:28677/147‑1 (MQ=255)
                            gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  2:85401/147‑1 (MQ=255)
                               cAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGggtggt                                                                                                                                                      >  6:49002/1‑116 (MQ=255)
                               cAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGggtggt                                                                                                                                                      <  5:49002/116‑1 (MQ=255)
                                 ggCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCCCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACg                                                                                                                   <  1:97709/149‑1 (MQ=255)
                                        ccAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATg                                                                                                            >  6:205714/1‑149 (MQ=255)
                                               tGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCCGGGTGGCCACCACCCgcac                                                                                                                                          >  2:252476/1‑110 (MQ=255)
                                               tGCTCGCCAGTGCACCGCGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCCGGGTGGCCACCACCCgcac                                                                                                                                          <  1:252476/113‑4 (MQ=255)
                                                  tCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCt                                                                                                   >  3:136223/1‑148 (MQ=255)
                                                  tCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCt                                                                                                   >  4:118087/1‑148 (MQ=255)
                                                      cAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGACACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAag                                                                                               >  7:180769/1‑148 (MQ=255)
                                                           cACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGcagcag                                                                                         <  1:86489/149‑1 (MQ=255)
                                                            aCCGGGCGGCATTAGCCCGCAGCGACTCGANCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAgg                                                                                        >  1:23694/1‑149 (MQ=255)
                                                            aCCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAgg                                                                                        >  1:18873/1‑149 (MQ=255)
                                                              cGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTc                                                                                                  >  2:53072/1‑137 (MQ=255)
                                                              cGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTc                                                                                                  <  1:53072/137‑1 (MQ=255)
                                                                 ctttgcatAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTg                                                                                   <  7:71732/142‑1 (MQ=255)
                                                                   ggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  4:136223/149‑1 (MQ=255)
                                                                   ggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGGGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  7:133813/149‑1 (MQ=255)
                                                                       ttAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAAt                                                                                                             <  4:83833/117‑1 (MQ=255)
                                                                       ttAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAAt                                                                                                             >  3:83833/1‑117 (MQ=255)
                                                                          gCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGtt                                                                           <  8:109414/148‑1 (MQ=255)
                                                                            ccGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgc                                                                        <  3:118087/149‑1 (MQ=255)
                                                                              gCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgcgc                                                                      <  5:37127/149‑1 (MQ=255)
                                                                                        gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  2:18873/149‑1 (MQ=255)
                                                                                        gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:15017/149‑1 (MQ=255)
                                                                                        gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  2:23694/149‑1 (MQ=255)
                                                                                        gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:205714/149‑1 (MQ=255)
                                                                                          gCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGGTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGc                                                          <  4:162510/149‑1 (MQ=255)
                                                                                               aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:184613/142‑1 (MQ=255)
                                                                                               aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            >  6:184613/1‑142 (MQ=255)
                                                                                               aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAagag                                                     >  7:6181/1‑149 (MQ=255)
                                                                                                      ccGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTc                                                                <  7:3557/131‑1 (MQ=255)
                                                                                                      ccGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTc                                                                >  8:3557/1‑131 (MQ=255)
                                                                                                       cGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  6:142519/113‑1 (MQ=255)
                                                                                                       cGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 >  5:142519/1‑113 (MQ=255)
                                                                                                                gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGa                                                <  5:94514/137‑1 (MQ=255)
                                                                                                                gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGa                                                >  6:94514/1‑137 (MQ=255)
                                                                                                                gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGccc                                    <  7:147377/149‑1 (MQ=255)
                                                                                                                gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGcc                                     >  8:147377/1‑148 (MQ=255)
                                                                                                                     cTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGAcc                               >  4:14078/1‑149 (MQ=255)
                                                                                                                           gTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATc                         <  7:33337/149‑1 (MQ=255)
                                                                                                                            ttACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCAc                                                              >  2:289100/1‑111 (MQ=255)
                                                                                                                            ttACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCAc                                                              <  1:289100/111‑1 (MQ=255)
                                                                                                                              aCCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCt                      >  5:125951/1‑149 (MQ=255)
                                                                                                                                       tgtg‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt             >  8:175125/1‑149 (MQ=255)
                                                                                                                                       cctgccGGGTGGCCACCACCC‑GCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt             >  2:241951/6‑149 (MQ=255)
                                                                                                                                         tg‑CGGGTGGTGGCTACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAATAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCttt           <  5:68997/149‑1 (MQ=255)
                                                                                                                                               gtggtgGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGc      >  5:97284/1‑149 (MQ=255)
                                                                                                                                                tggtggCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGcc     >  6:292183/1‑149 (MQ=255)
                                                                                                                                                   tggCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTg  <  8:6181/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGGTGCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTG  >  NC_002947/583250‑583544

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