Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,931,935 G→T 7.1% A529D (GCC→GAC) ‡ PP_3466 ← ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,931,9350GT 7.1% 214.6 / 3.0 69A529D (GCC→GAC) ‡PP_3466ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TCGATGATGCCAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGG  >  NC_002947/3931790‑3932084
                                                                                                                                                 |                                                                                                                                                     
tCGATGATGCCAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGaa                                                                                                                                                   <  1:993253/150‑1 (MQ=255)
    tgatgCCAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCTCGTGCGGCCGGCGAACAGc                                                                                                                                               >  2:1611336/1‑150 (MQ=255)
          cAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTgg                                                                                                                                            >  1:115458/1‑147 (MQ=255)
          cAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTgg                                                                                                                                            <  2:115458/147‑1 (MQ=255)
              gcgcACGTCAGTGCCCCTGTTTGATCACCTCCTGCGCCAGTCCCTACAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATgg                                                                                                                                     <  1:251856/150‑1 (MQ=255)
                gcACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGAGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCATCGACCAGCTGTTCgccggcc                                                                                                                                   <  2:611616/150‑7 (MQ=255)
                       aGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCa                                                                                                                            >  1:565517/1‑150 (MQ=255)
                           cccATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTc                                                                                                                        >  2:21131/1‑150 (MQ=255)
                            ccATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACtt                                                                                                                         <  2:1339094/148‑1 (MQ=255)
                            ccATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACtt                                                                                                                         >  1:1339094/1‑148 (MQ=255)
                                     aTCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTg                                                                                                              <  2:395373/150‑1 (MQ=255)
                                     aTCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTg                                                                                                              <  1:21131/150‑1 (MQ=255)
                                      tCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGCGCGGCCATCGACCAGCTGTTCGCCGGCcgcacgc                                                                                                                             <  1:1123592/134‑5 (MQ=38)
                                      tCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGCGCGGCCATCGACCAGCTGTTCGCCGGCcgcacgc                                                                                                                             >  2:1123592/1‑130 (MQ=37)
                                      tCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGACGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTgg                                                                                                             >  2:1549639/1‑150 (MQ=255)
                                          ctcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATTGCCGCGa                                                                                                                               <  2:881254/128‑1 (MQ=255)
                                          ctcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGa                                                                                                                               >  1:881254/1‑128 (MQ=255)
                                          ctcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTc                                                                                                         >  1:272103/1‑150 (MQ=255)
                                           tcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACt                                                                                                                          <  1:183205/132‑1 (MQ=255)
                                           tcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACt                                                                                                                          >  2:183205/1‑132 (MQ=255)
                                           tcctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCAATCACTTCACGCTCGGTGGCCTCa                                                                                                        <  1:464050/150‑1 (MQ=255)
                                             ctGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATc                                                                                                      <  2:825320/150‑1 (MQ=255)
                                                   cAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCCTCCatcgc                                                                                                <  2:565517/150‑5 (MQ=255)
                                                   cAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGc                                                                                                >  1:960414/1‑150 (MQ=255)
                                                   cAGTAACTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGc                                                                                                >  2:1180244/1‑150 (MQ=255)
                                                     gTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCg                                                                                              >  1:267427/1‑150 (MQ=255)
                                                     cttcgcgCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCg                                                                                              >  2:334666/7‑150 (MQ=255)
                                                      tACCTGTAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGa                                                                                             >  1:1204148/1‑150 (MQ=255)
                                                         cTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGt                                                                                          >  1:267040/1‑150 (MQ=255)
                                                          tGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTg                                                                                         <  1:1046107/150‑1 (MQ=255)
                                                            cAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGc                                                                                       <  2:264616/150‑1 (MQ=255)
                                                             aaCTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGaa                                                                                                                                                   >  2:121232/1‑89 (MQ=255)
                                                             aaCTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGaa                                                                                                                                                   <  1:121232/89‑1 (MQ=255)
                                                               cTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTc                                                                                    <  1:334666/150‑1 (MQ=255)
                                                               cTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCATCGACCAGCTg                                                                                                                                             <  1:1225637/93‑1 (MQ=255)
                                                               cTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCATCGACCAGCTg                                                                                                                                             >  2:1225637/1‑93 (MQ=255)
                                                                  gCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTc                                                                                 <  1:1637324/150‑1 (MQ=255)
                                                                          gCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACCGCCGAAGTGGCCTCGTCCAGCACCa                                                                         <  1:1611336/150‑1 (MQ=255)
                                                                              gtacAGCTGCAGGTCGGCATGGGCCCGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGAt                                                                     <  2:267040/147‑1 (MQ=255)
                                                                                gCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCAc                                                                                                                           >  1:814605/1‑94 (MQ=255)
                                                                                gCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCAc                                                                                                                           <  2:814605/94‑1 (MQ=255)
                                                                                       cAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGc                                                            <  2:960414/150‑1 (MQ=255)
                                                                                        aGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCAc                                                                                                                           <  2:1451898/86‑1 (MQ=255)
                                                                                        aGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCAc                                                                                                                           >  1:1451898/1‑86 (MQ=255)
                                                                                             ggCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGACCGCGCTCACTTCCCGCCCTGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCaa                                                      >  1:115607/1‑150 (MQ=255)
                                                                                                aTCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGcc                                                                                      <  2:284446/115‑1 (MQ=255)
                                                                                                aTCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGcc                                                                                      >  1:284446/1‑115 (MQ=255)
                                                                                                   ggCCAGGGTCGAAGCCCGTTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCAt                                                                                                       <  1:1518243/95‑1 (MQ=255)
                                                                                                   ggCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCAt                                                                                                       >  2:1518243/1‑95 (MQ=255)
                                                                                                    gCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCa                                               <  2:995756/150‑1 (MQ=255)
                                                                                                      cAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACg                                             >  2:646258/1‑150 (MQ=255)
                                                                                                       aGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGa                                            <  2:272103/150‑1 (MQ=255)
                                                                                                             gAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATg                                      >  1:230187/1‑150 (MQ=255)
                                                                                                               aGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACt                                                                                                  >  2:319400/1‑88 (MQ=255)
                                                                                                               aGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACt                                                                                                  <  1:319400/88‑1 (MQ=255)
                                                                                                                 cccGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGa                                  >  2:568772/1‑150 (MQ=255)
                                                                                                                  ccGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGAt                                 >  1:881823/1‑150 (MQ=255)
                                                                                                                         cTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCgttgt                          <  1:1415400/150‑1 (MQ=255)
                                                                                                                          tGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCgttgtt                         >  1:168676/1‑150 (MQ=255)
                                                                                                                          tGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCgttgtt                         >  1:472761/1‑150 (MQ=255)
                                                                                                                           gATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGttt                        >  2:570799/1‑150 (MQ=255)
                                                                                                                            aTCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCt                                                          >  2:193364/1‑115 (MQ=255)
                                                                                                                            aTCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCt                                                          <  1:193364/115‑1 (MQ=255)
                                                                                                                              cAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTg                     <  2:472761/150‑1 (MQ=255)
                                                                                                                                ggATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCAGCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATgg                                     <  2:947990/132‑1 (MQ=255)
                                                                                                                                ggATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATgg                                     >  1:947990/1‑132 (MQ=255)
                                                                                                                                   tGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGcccc                >  1:1371030/1‑150 (MQ=255)
                                                                                                                                        gTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGc                               >  1:1456976/1‑130 (MQ=255)
                                                                                                                                        gTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGc                               <  2:1456976/130‑1 (MQ=255)
                                                                                                                                             gccggcGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCCCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAActg      <  1:1180244/150‑1 (MQ=255)
                                                                                                                                               cggcGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCgg                                                               <  2:604897/91‑1 (MQ=255)
                                                                                                                                               cggcGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCgg                                                               >  1:604897/1‑91 (MQ=255)
                                                                                                                                                 gcGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCTATGCGTTGTTTCTGGCCCCCGGACAACTGCTgg  <  1:1549639/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCGATGATGCCAACGCGCACGTCAGTGCCCATGTTTGATCACCTCCTGCGCCAGTACCTGCAACTGGCCGTCATGCAGTTGCAGCTGCAGGTCGGCATCGGCCAGGGTCGAAGCCCGGTGGCTGATCAGGATGCGCGTGCGGCCGGCGAACAGCTGGTCGATGGCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGG  >  NC_002947/3931790‑3932084

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