Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,165,154 C→T 6.0% Q316Q (CAG→CAA)  PP_3666 ← MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,165,1540CT 6.0% 199.3 / 2.4 66Q316Q (CAG→CAA) PP_3666MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

TCCCGTCGTTCATCGCCAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACC  >  NC_002947/4165006‑4165301
                                                                                                                                                    |                                                                                                                                                   
tCCCGTCGTTCATCGCCAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAAc                                                                                                                                                     <  2:230918/149‑1 (MQ=255)
     tcgtTCATCGCCAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTc                                                                                                                                                <  1:95324/149‑1 (MQ=255)
               ccAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAgccgc                                                                                                                                      >  1:235847/1‑149 (MQ=255)
               ccAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATCGTCTTGCGACCAAACTGGTCGGACAgccgc                                                                                                                                      >  1:146998/1‑149 (MQ=255)
                       gCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCa                                                                                                                               <  3:293955/148‑1 (MQ=255)
                           cGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCa                                                                                                                                  <  2:292773/141‑1 (MQ=255)
                           cGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCa                                                                                                                                  >  1:292773/1‑141 (MQ=255)
                            gAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGa                                                                                                                         <  2:47737/149‑1 (MQ=255)
                                        gATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGata                                                                                                              >  4:263832/1‑148 (MQ=255)
                                             ccAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAGTTTGGTCgca                                                                                                                                             >  7:77782/1‑110 (MQ=255)
                                             ccAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAGTTTGGTCgca                                                                                                                                             <  8:77782/112‑3 (MQ=255)
                                              cAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAg                                                                                                                    >  5:204110/1‑136 (MQ=255)
                                              cAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAg                                                                                                                    <  6:204110/136‑1 (MQ=255)
                                                gAATCACCAGCAGCGGCTGGCCGTCCGGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGc                                                                                                      >  7:87212/1‑148 (MQ=255)
                                                   tCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACTGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGc                                                                                                                   <  7:235882/132‑1 (MQ=255)
                                                   tCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGc                                                                                                                   >  8:235882/1‑132 (MQ=255)
                                                       cagcagCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGTTGGCCAGCGaca                                                                                              >  6:320502/1‑149 (MQ=255)
                                                        agcagcGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAGTTTGGTCgcaa                                                                                                                                            >  3:124397/1‑99 (MQ=255)
                                                        agcagcGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAGTTTGGTCgcaa                                                                                                                                            <  4:124397/102‑4 (MQ=255)
                                                            gcGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCAACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTc                                                                                         >  7:229442/1‑149 (MQ=255)
                                                            gcGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCAACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTc                                                                                         >  7:225126/1‑149 (MQ=255)
                                                                cTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATTTAGGTGGCCAGCGACACCGTCGAgg                                                                                     <  3:263832/149‑1 (MQ=255)
                                                                    ccGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGACGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGAt                                                                                 <  7:257987/149‑1 (MQ=255)
                                                                     cGTCCAGCAGGTGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATg                                                                                <  6:251478/149‑1 (MQ=255)
                                                                     cGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGAt                                                                                 >  8:92623/1‑148 (MQ=255)
                                                                         cagcagGCGAAACAACGGCACGGCCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAAc                                                                             <  5:204179/148‑1 (MQ=255)
                                                                          agcagGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACg                                                                            >  1:157375/1‑148 (MQ=255)
                                                                          agcagGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACg                                                                            >  8:213132/1‑148 (MQ=255)
                                                                               gCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATccc                                                                      >  2:291814/1‑149 (MQ=255)
                                                                                 gAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCg                                                                                        <  3:242178/129‑1 (MQ=255)
                                                                                 gAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCg                                                                                        >  4:242178/1‑129 (MQ=255)
                                                                                  aaaCAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCcgct                                                                   >  2:222308/1‑149 (MQ=255)
                                                                                   aacaacGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGAc                                                                                                                                            >  8:266877/1‑75 (MQ=255)
                                                                                   aacaacGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGAc                                                                                                                                            <  7:266877/75‑1 (MQ=255)
                                                                                         ggCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCAt                                                            >  6:285708/1‑149 (MQ=255)
                                                                                                       cAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTg                                               >  1:210843/1‑148 (MQ=255)
                                                                                                          gAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCAGTTCCCTGGa                                             >  7:258936/1‑147 (MQ=255)
                                                                                                             caacaGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACaggt                                        <  8:87212/149‑1 (MQ=255)
                                                                                                             caacaGCGAGGCCACCACCAGCATGGTCTTGCAACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACaggt                                        <  8:229442/149‑1 (MQ=255)
                                                                                                             caacaGCGAGGCCACCAAGAGCATGGTGTTGCAACCAATCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACaggt                                        <  8:225126/149‑1 (MQ=255)
                                                                                                              aacaGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACagg                                         >  7:286055/1‑147 (MQ=255)
                                                                                                               acaGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACaggtag                                      >  5:258095/1‑149 (MQ=255)
                                                                                                                   cGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTc                                   <  1:58310/148‑1 (MQ=255)
                                                                                                                   cGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTc                                   >  2:58310/1‑148 (MQ=255)
                                                                                                                   cGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCg                                  <  6:258095/149‑1 (MQ=255)
                                                                                                                        ccaccaccaGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATg                             <  7:213132/149‑1 (MQ=255)
                                                                                                                        ccaccaccaGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATg                             >  6:327036/1‑149 (MQ=255)
                                                                                                                         caccaccaGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGt                            <  2:157375/149‑1 (MQ=255)
                                                                                                                            caccaGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAAc                         >  1:266009/1‑149 (MQ=255)
                                                                                                                               cagcagGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGa                      <  8:32068/149‑1 (MQ=255)
                                                                                                                               caGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGa                      <  5:285708/149‑1 (MQ=255)
                                                                                                                               caGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGa                      <  5:327036/149‑1 (MQ=255)
                                                                                                                               caGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGa                      <  7:92623/149‑1 (MQ=255)
                                                                                                                                         ttGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCAt                              >  8:121689/1‑131 (MQ=255)
                                                                                                                                         ttGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCAt                              <  7:121689/131‑1 (MQ=255)
                                                                                                                                         ttGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAgg              >  5:273298/1‑147 (MQ=255)
                                                                                                                                         ttGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTa            >  3:80059/1‑149 (MQ=255)
                                                                                                                                           gCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTa            <  2:210843/147‑1 (MQ=255)
                                                                                                                                            cGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGaa         <  4:120294/149‑1 (MQ=255)
                                                                                                                                             gACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAAc        >  6:2295/1‑149 (MQ=255)
                                                                                                                                             gACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAAc        <  2:235847/149‑1 (MQ=255)
                                                                                                                                                 aaaCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGc                    <  5:274039/133‑1 (MQ=255)
                                                                                                                                                 aaaCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGc                    >  6:274039/1‑133 (MQ=255)
                                                                                                                                                 aaaCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCaa    <  6:273298/149‑1 (MQ=255)
                                                                                                                                                 aaaCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGGCGGCATGTAACTGAGCAACAGGTAGAACGCaa    <  8:258936/149‑1 (MQ=255)
                                                                                                                                                   aCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAAcc  <  4:80059/149‑1 (MQ=255)
                                                                                                                                                   aCTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAAcc  >  8:126252/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCCCGTCGTTCATCGCCAGCATCGCACCGAAAAGGATCTGGATGGCCAGAATCACCAGCAGCGGCTGGCCGTCCAGCAGGCGAAACAACGGCACGGTCAACCCCAGGAACAACAGCGAGGCCACCACCAGCATGGTCTTGCGACCAAACTGGTCGGACAGCCGCCCCATCAGGAAGATCAAGCCGATATAGGTGGCCAGCGACACCGTCGAGGCGATGAACGAATCCCGCTCGCTCATGCCCATTTCACTGGACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACC  >  NC_002947/4165006‑4165301

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