Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 780,004 A→T 7.6% V156E (GTG→GAG)  PP_0670 ← bile acid/Na+ symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947780,0040AT 7.6% 165.5 / 3.4 62V156E (GTG→GAG) PP_0670bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base T (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGCG  >  NC_002947/779859‑780151
                                                                                                                                                  |                                                                                                                                                   
cggcaccggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCt                                                                                                                                                  >  3:24136/1‑149 (MQ=255)
      cggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGgctgc                                                                                                                                            >  2:186303/1‑149 (MQ=255)
       ggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGgctgct                                                                                                                                           <  7:258461/149‑1 (MQ=255)
              aCGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCAccc                                                                                                                                    >  7:295401/1‑149 (MQ=255)
                 cccAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGc                                                                                                                                 <  1:186303/149‑1 (MQ=255)
                    aGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCAt                                                                                                                               >  1:194996/1‑148 (MQ=255)
                     ggCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCa                                                                                                                             <  2:285080/149‑1 (MQ=255)
                      gCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAg                                                                                                                            <  3:303343/149‑1 (MQ=255)
                       cgtcgtTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAgg                                                                                                                           <  3:37869/149‑1 (MQ=255)
                                aGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTac                                                                                                                  >  5:227693/1‑149 (MQ=255)
                                      cgACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACACCACGAAGGCTTCCAGGAACGGGGTGATGTGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCgg                                                                                                             <  1:272028/149‑1 (MQ=255)
                                       gACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATCGCCACCTGGCTGCTGCCCCCCAGCATCAGGGCCAGGTACACCGgca                                                                                                           >  3:157738/1‑149 (MQ=255)
                                        aCTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGATCACCGCGTGGGCCATCGGCATCTCGTTCAACAGCACGAAGGCTTCCTTGAACGTGGTGTTG‑GATATTGCCTCCTGGCTGCTTCCACCCATGCTCAGGGCCAGGTACACCGgcag                                                                                                          <  7:54889/149‑1 (MQ=255)
                                              cgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCAc                                                                                                    <  4:24136/149‑1 (MQ=255)
                                              cgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATCGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCAc                                                                                                    <  4:157738/149‑1 (MQ=255)
                                               accaGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCAc                                                                                                    >  2:280818/4‑148 (MQ=255)
                                                cgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCAc                                                                                                    >  6:62224/1‑147 (MQ=255)
                                                     gcCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTgc                                                                                             >  2:14459/1‑149 (MQ=255)
                                                     gcCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGACCAGGTACACCGGCAGCAGCACCAGTTgc                                                                                             >  5:324632/1‑149 (MQ=255)
                                                        gggCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTgcgc                                                                                           >  8:17264/1‑148 (MQ=255)
                                                                gcgcTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATaa                                                                                  >  6:193658/1‑149 (MQ=255)
                                                                 cgcTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAAc                                                                                 <  5:193658/149‑1 (MQ=255)
                                                                    tGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGTCTGCTGCCACACAGCATCAGGGCCAGGTACACCGGCACCGGCACCAGTTGCTCCAGCAATAACAgg                                                                              >  5:238876/1‑149 (MQ=255)
                                                                      gTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAgggg                                                                            >  6:258907/1‑149 (MQ=255)
                                                                      gTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAgggg                                                                            >  6:300140/1‑149 (MQ=255)
                                                                      gTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAgg                                                                              <  3:135402/147‑1 (MQ=255)
                                                                      gTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAgg                                                                              <  1:14459/147‑1 (MQ=255)
                                                                          gCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGtt                                                                        <  8:295401/149‑1 (MQ=255)
                                                                             ccGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCg                                                                     <  2:194996/149‑1 (MQ=255)
                                                                                 gAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTTCGGa                                                                   >  8:106284/1‑147 (MQ=255)
                                                                                   gggCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAgg                                                               <  6:227693/149‑1 (MQ=255)
                                                                                   gggCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAgg                                                               <  6:91791/149‑1 (MQ=255)
                                                                                        aTCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACACAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAg                                                          >  3:288673/1‑149 (MQ=255)
                                                                                            gCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGt                                                      <  5:300140/149‑1 (MQ=255)
                                                                                            gCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGt                                                      <  4:288673/149‑1 (MQ=255)
                                                                                            gCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGt                                                      >  3:78460/1‑149 (MQ=255)
                                                                                              aaCACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAg                                                                <  1:271795/137‑1 (MQ=255)
                                                                                              aaCACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAg                                                                >  2:271795/1‑137 (MQ=255)
                                                                                              aaCACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGa                                                    <  5:62224/149‑1 (MQ=255)
                                                                                                    aTCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAg                                                          >  1:15042/1‑137 (MQ=255)
                                                                                                    aTCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAg                                                          <  2:15042/137‑1 (MQ=255)
                                                                                                       aaCAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGTTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCttt                                           <  1:4927/149‑1 (MQ=255)
                                                                                                              cGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAAtg                                    <  7:17264/149‑1 (MQ=255)
                                                                                                               gAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCTCCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAAtgt                                   >  4:167849/1‑149 (MQ=255)
                                                                                                               gAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCTCCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAAtgt                                   >  8:184289/1‑149 (MQ=255)
                                                                                                               gAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAAtgt                                   >  1:127616/1‑149 (MQ=255)
                                                                                                                aaGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAAtgtg                                  >  6:702/1‑149 (MQ=255)
                                                                                                                 aGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTgg                                 >  3:155246/1‑149 (MQ=255)
                                                                                                                     ttCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTg                              >  4:252836/1‑148 (MQ=255)
                                                                                                                     ttCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTg                              >  6:196464/1‑148 (MQ=255)
                                                                                                                       ccAGGAACCGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGaa                            <  8:326940/148‑1 (MQ=255)
                                                                                                                             aCGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAgg                                                         <  4:302912/113‑1 (MQ=255)
                                                                                                                             aCGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAgg                                                         >  3:302912/1‑113 (MQ=255)
                                                                                                                                         gagaTTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGAt          >  7:300052/1‑149 (MQ=255)
                                                                                                                                          agaTTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATg         >  3:142595/1‑149 (MQ=255)
                                                                                                                                            aTTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCa       <  5:196464/149‑1 (MQ=255)
                                                                                                                                               gCCTCCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAgg     <  7:184289/148‑1 (MQ=255)
                                                                                                                                               gCCTCCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAgg     <  3:167849/148‑1 (MQ=255)
                                                                                                                                                ccACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGc   <  5:702/149‑1 (MQ=255)
                                                                                                                                                 cACCTGGCTGCTGCCACCCAGGATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGCg  >  3:305208/1‑149 (MQ=255)
                                                                                                                                                 cACCTGGCTGCTGCCACCCAGCTTCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGCg  >  1:235080/1‑149 (MQ=255)
                                                                                                                                                 cACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTggg                                <  8:287123/119‑1 (MQ=255)
                                                                                                                                                  aCCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGGGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGCg  >  5:369438/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATG‑GAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAGCAATAACAGGGGGGTTGCGGACAGGGTCAGGCGTGAGTCGCCTTTGCCAATGTGGGTGAACACCACCACGTAGTCGATGCAGGGCG  >  NC_002947/779859‑780151

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