Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,264,962 A→T 10.4% D44V (GAC→GTC)  vanB → vanillate O‑demethylase oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,264,9620AT10.4% 218.9 / 15.0 77D44V (GAC→GTC) vanBvanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base T (6/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTGCG  >  NC_002947/4264816‑4265110
                                                                                                                                                  |                                                                                                                                                    
ctgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACgg                                                                                                                                                   >  1:584184/1‑150 (MQ=255)
     gAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTgg                                                                                                                                              <  1:862393/150‑1 (MQ=255)
     gAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTgg                                                                                                                                              <  1:256576/150‑1 (MQ=255)
        cAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgc                                                                                                                                           <  2:432623/150‑1 (MQ=255)
             gACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAg                                                                                                                                      <  2:1589108/150‑1 (MQ=255)
                 tGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTAtt                                                                                                                                  >  1:1027768/1‑150 (MQ=255)
                 tGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTAtt                                                                                                                                  >  2:830057/1‑150 (MQ=255)
                   aTCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTa                                                                                                                                    <  2:603215/146‑1 (MQ=255)
                   aTCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTa                                                                                                                                    >  1:603215/1‑146 (MQ=255)
                       aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtgt                                                                                                                            <  1:1680723/150‑1 (MQ=255)
                       aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtgt                                                                                                                            <  1:1227708/150‑1 (MQ=255)
                       aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtgt                                                                                                                            >  1:1932248/1‑150 (MQ=255)
                       aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtgt                                                                                                                            <  2:1824441/150‑1 (MQ=255)
                            gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGcc                                                                                                                                        <  2:1241464/133‑1 (MQ=255)
                            gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGcc                                                                                                                                        >  1:1241464/1‑133 (MQ=255)
                            gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAac                                                                                                                       <  1:1914005/150‑1 (MQ=255)
                                  gTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                              <  2:64296/137‑1 (MQ=255)
                                  gTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                              >  1:64296/1‑137 (MQ=255)
                                       ccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCCCCCCgagga                                                                                                            >  1:274570/1‑150 (MQ=255)
                                       ccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCgaaga                                                                                                            >  2:889426/1‑150 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgca                                                                                                                                          >  1:971565/1‑112 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgca                                                                                                                                          <  2:971565/113‑2 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAaccac                                                                                                                    >  2:462973/1‑135 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAaccac                                                                                                                    <  1:462973/135‑1 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          >  2:437523/1‑145 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          <  1:437523/145‑1 (MQ=255)
                                              gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCAt                                                                                                     <  1:779706/150‑1 (MQ=255)
                                                   aGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAacc                                                                                                                      <  2:145534/128‑1 (MQ=255)
                                                   aGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAacc                                                                                                                      >  1:145534/1‑128 (MQ=255)
                                                          ggTCTCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGAtt                                                                                         >  2:947186/1‑150 (MQ=255)
                                                           gTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                                                        <  1:1079735/150‑1 (MQ=255)
                                                             aTCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGc                                                                                      <  2:1881670/150‑1 (MQ=255)
                                                                tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                              >  2:491145/1‑107 (MQ=255)
                                                                tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                              <  1:491145/107‑1 (MQ=255)
                                                                 gCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCa                                                                                                      <  1:437281/130‑1 (MQ=255)
                                                                 gCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCa                                                                                                      >  2:437281/1‑130 (MQ=255)
                                                                 gCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                                                        >  1:848105/1‑144 (MQ=255)
                                                                 gCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                                                        <  2:848105/144‑1 (MQ=255)
                                                                   aGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                         <  1:1420096/93‑3 (MQ=38)
                                                                   aGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                         >  2:1420096/1‑91 (MQ=38)
                                                                          aGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAcc                                                                         >  1:88457/1‑150 (MQ=255)
                                                                           gCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCa                                                                                                                                       <  1:1423388/87‑1 (MQ=255)
                                                                           gCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCa                                                                                                                                       >  2:1423388/1‑87 (MQ=255)
                                                                               gCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCt                                                                    >  1:388950/1‑150 (MQ=255)
                                                                                 cgcgGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          <  1:440417/110‑1 (MQ=255)
                                                                                 cgcgGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          >  2:440417/1‑110 (MQ=255)
                                                                                       aGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGc                                                                                                        >  2:188571/1‑106 (MQ=255)
                                                                                       aGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGc                                                                                                        <  1:188571/106‑1 (MQ=255)
                                                                                       aGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                                            <  1:830057/150‑1 (MQ=255)
                                                                                        gATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGcagg                                                                                                                                              >  1:1481578/1‑64 (MQ=38)
                                                                                        gATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAgg                                                                                                                                              <  2:1481578/67‑1 (MQ=38)
                                                                                         aTGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCACCCACCCCGACGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGTCTTCGCGGGGCg                                                          <  2:807081/150‑1 (MQ=255)
                                                                                          tGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCgg                                                         <  2:1159321/150‑1 (MQ=255)
                                                                                           ggCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCa                                                                                                                         >  1:537066/1‑85 (MQ=255)
                                                                                           ggCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCa                                                                                                                         <  2:537066/85‑1 (MQ=255)
                                                                                           ggCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCACTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGt                                                        >  1:883206/1‑150 (MQ=255)
                                                                                              aGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          >  1:1940837/1‑97 (MQ=255)
                                                                                              aGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                          <  2:1940837/97‑1 (MQ=255)
                                                                                                 ctgctgCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGt                                                  >  2:638320/1‑150 (MQ=255)
                                                                                                  tgctgcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCt                                                                    >  2:191464/1‑131 (MQ=255)
                                                                                                  tgctgcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCt                                                                    <  1:191464/131‑1 (MQ=255)
                                                                                                                aGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAg                                   <  1:889426/150‑1 (MQ=255)
                                                                                                                aGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGCAGCCTGCACGAACAg                                   >  2:1347395/1‑150 (MQ=255)
                                                                                                                 gcgcCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAgg                                  <  2:1702950/150‑1 (MQ=255)
                                                                                                                 gcgcCGGCGCCCACACCCACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGAACACCTTCGGAACACCGACACCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAgg                                  <  1:1592909/150‑1 (MQ=255)
                                                                                                                       gcgcCCATCTCGACGCGCACCCGCCCGACGGGCTGGTGCGCCAGTATTCTGTGTTCAACCACCCCGACGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCGCTTAGCCTGCACGAACAGGTGCAAg                            <  2:907903/150‑1 (MQ=255)
                                                                                                                           ccATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCgg                        >  2:1456607/1‑150 (MQ=255)
                                                                                                                                   aCGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGc                >  1:1354432/1‑150 (MQ=255)
                                                                                                                                     gcgCACCAGCCCGTCGGGCAGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTg              >  2:660857/3‑150 (MQ=255)
                                                                                                                                     gcgCACCAGCCCGTCGGGCAGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGCTGATATAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTg              >  1:1475769/3‑150 (MQ=255)
                                                                                                                                        cACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                                                        <  2:1130999/73‑1 (MQ=255)
                                                                                                                                        cACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                                                        >  1:1130999/1‑73 (MQ=255)
                                                                                                                                        cACCAGCCCGTCGGGCAGGTGCGCTAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGTGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCGGGTGTCCGGcggggt           >  1:702570/1‑145 (MQ=255)
                                                                                                                                         aCCAGCCCGTCGGGCAGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTa          >  2:1840639/1‑150 (MQ=255)
                                                                                                                                           cTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAAAGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATc        >  2:332059/1‑150 (MQ=255)
                                                                                                                                            tGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCa       <  1:947186/150‑1 (MQ=255)
                                                                                                                                                 gACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTgcg  >  2:1537578/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTGCG  >  NC_002947/4264816‑4265110

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