Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,053,520 A→T 6.2% Q17L (CAG→CTG)  PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,053,5200AT 6.2% 300.9 / 8.8 97Q17L (CAG→CTG) PP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base T (3/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGGGC  >  NC_002947/4053371‑4053669
                                                                                                                                                        |                                                                                                                                                       
aTACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCa                                                                                                                                                         <  2:2186984/150‑1 (MQ=255)
         tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCTGGCCGGCgc                                                                                                                                                <  2:1294140/150‑2 (MQ=255)
           gCCATTTCCACGGCAACGCTCTTTTGCACCCAGACAAGGGCTTGCCTCCAAACCGACCTGACAAAAACCGTTTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTgg                                                                                                                                              <  1:2156222/150‑1 (MQ=255)
            ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGc                                                                                                                                                  >  1:544462/1‑145 (MQ=255)
            ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGc                                                                                                                                                  <  2:544462/145‑1 (MQ=255)
                 tCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGa                                                                                                                                        <  1:1259858/150‑1 (MQ=255)
                  ccACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAAATCGACACCCT‑CGTCGT‑ACGCGCCGGCGAGGCCGGCGTGGCCATGa                                                                                                                                        <  2:1376465/150‑1 (MQ=255)
                   cgcGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCATTTCGATCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGc                                                                                                                                      <  2:658760/148‑1 (MQ=255)
                         aTCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAacac                                                                                                                                >  2:369859/1‑150 (MQ=255)
                         aTCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCAC‑AGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAacac                                                                                                                                <  1:776946/150‑1 (MQ=255)
                            ccTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTg                                                                                                                             >  1:2173504/1‑150 (MQ=255)
                              tAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGa                                                                                                                                        >  2:1928595/1‑137 (MQ=255)
                              tAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGa                                                                                                                                        <  1:1928595/137‑1 (MQ=255)
                                gTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCa                                                                                                                         >  2:897468/1‑150 (MQ=255)
                                gTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCa                                                                                                                         >  1:2173590/1‑150 (MQ=255)
                                  cTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAg                                                                                                                       <  1:1876282/150‑1 (MQ=255)
                                      cTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTg                                                                                                                   <  1:1386102/150‑1 (MQ=255)
                                          aGGCAACGGCTGTCCTCCAAGCCGACCTGACAGAAGCCGCCTGAAGAAGGGCGCGGACATAACACTGAAGAACCTCG‑AAATCGGCACCCG‑CGTCGT‑CGGCGCCGGCCAGGCGCGCATGGCCACGAGCGAAAACCTCAGCAAGccgggccg                                                                                                               >  2:668979/1‑143 (MQ=255)
                                          aGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  <  1:2243941/147‑1 (MQ=255)
                                          aGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  >  2:2243941/1‑147 (MQ=255)
                                          aGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGt                                                                                                               <  2:1326728/150‑1 (MQ=255)
                                          aGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGt                                                                                                               >  2:1579591/1‑150 (MQ=255)
                                            gCAACGGCGGTCCTCCAAACCTACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  <  1:663802/145‑1 (MQ=255)
                                            gCAACGGCGGTCCTCCAAACCTACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  >  2:663802/1‑145 (MQ=255)
                                            gCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCttagagcac                                                                                                             >  2:489765/1‑143 (MQ=255)
                                                 ggcggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAc                                                                                                        >  1:1910339/1‑150 (MQ=255)
                                                 ggcggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAc                                                                                                        <  2:2318089/150‑1 (MQ=255)
                                                     gTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTg‑g                                                                                                   <  1:128613/150‑1 (MQ=255)
                                                     gTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTg‑g                                                                                                   <  2:1794187/150‑1 (MQ=255)
                                                      tcctccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑Gt                                                                                                  <  2:409900/150‑1 (MQ=255)
                                                         tccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACt                                                                                               >  1:1296914/1‑150 (MQ=255)
                                                            gaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            >  2:1605707/2‑150 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAATGGCTCGGACATGACACTGAACAACCGCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCTTGGCCATGAGCGAACACCTGAGCAAGCATGGCGTGCCGCACCTG‑GTACTGGa                                                                                            <  1:489765/150‑1 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCTGGCCGGCgcc                                                                                                                                               >  2:1963957/1‑98 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCTGGCCGGCgcc                                                                                                                                               <  1:1963957/100‑3 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            <  1:1579591/150‑1 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            <  2:2173590/150‑1 (MQ=255)
                                                             aaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAg                                                                                           >  2:1860528/1‑150 (MQ=255)
                                                              aCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAg‑c                                                                                         <  2:887882/150‑1 (MQ=255)
                                                                   ccTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAAc                                                                                    >  2:1428619/1‑150 (MQ=255)
                                                                     tGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCg                                                                                  >  2:2184456/1‑150 (MQ=255)
                                                                      gACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGc                                                                                 >  2:97871/1‑150 (MQ=255)
                                                                       aCAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCa                                                                                >  1:1953176/1‑150 (MQ=255)
                                                                         aaaaaCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATc                                                                              <  2:67653/150‑1 (MQ=255)
                                                                           aaaCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTg‑g                                                                                                   <  1:379041/128‑1 (MQ=255)
                                                                           aaaCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTg‑g                                                                                                   >  2:379041/1‑128 (MQ=255)
                                                                             aCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCgcgc                                                                          >  1:1198184/1‑150 (MQ=255)
                                                                               cGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAg                                                                        <  1:1428619/150‑1 (MQ=255)
                                                                                gCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAgg                                                                       >  2:2059916/1‑150 (MQ=255)
                                                                                 ccTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCTGGCCgg                                                                                                                                                   <  1:2081173/75‑1 (MQ=255)
                                                                                 ccTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCTGGCCgg                                                                                                                                                   >  2:2081173/1‑75 (MQ=255)
                                                                                 ccTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCCACCGCATCGCGCAGGc                                                                      >  2:696659/1‑150 (MQ=255)
                                                                                     aaCAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTgg                                                                  <  1:994884/150‑1 (MQ=255)
                                                                                     aaCAAGGGCTCGGACATGACACTGAACAACCTAC‑ACATCCACATTCTGTGTC‑T‑CTGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTgg                                                                  >  1:2156609/1‑150 (MQ=255)
                                                                                      aCAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGc                                                                 <  2:2173504/150‑1 (MQ=255)
                                                                                       cAAGGGCTCGGACATGACACTGAACAACCTCG‑CACTAGACA‑CCG‑CGCCATGGGGCGCCGGCCAGGCCGTGTTGGGCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCCGCAG‑AT‑CGGAAGAGGCAACCTCAGCGCGCAGGCCTGGCg                                                                >  1:405460/1‑150 (MQ=255)
                                                                                          gggCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTAc                                                             <  1:2184456/150‑1 (MQ=255)
                                                                                           ggCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTAcc                                                            <  1:1860528/150‑1 (MQ=255)
                                                                                                 gACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGc                                                      >  1:273112/1‑150 (MQ=255)
                                                                                                    aTGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTgg                                                   <  2:1198184/150‑1 (MQ=255)
                                                                                                    aTGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTgg                                                   <  1:369859/150‑1 (MQ=255)
                                                                                                    aTGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTgg                                                   <  1:841281/150‑1 (MQ=255)
                                                                                                     tGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTggg                                                  <  2:1953176/150‑1 (MQ=255)
                                                                                                      gACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            <  1:885486/108‑1 (MQ=255)
                                                                                                      gACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            >  2:885486/1‑108 (MQ=255)
                                                                                                         acTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCAAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                              >  1:1701189/1‑150 (MQ=255)
                                                                                                           tGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  >  1:1542004/1‑82 (MQ=255)
                                                                                                           tGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTgg                                                                                                                  <  2:1542004/82‑1 (MQ=255)
                                                                                                               atctcTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGCCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGt                                        <  1:2260620/146‑1 (MQ=255)
                                                                                                                aacCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTCTGCTTGCCGCACCTC‑GGTCGGTTG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTc                                       <  1:284542/150‑1 (MQ=255)
                                                                                                                 acCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTcc                                      >  2:972612/1‑149 (MQ=255)
                                                                                                                  ccTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGc                                     >  2:192404/1‑150 (MQ=255)
                                                                                                                   cTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGcc                                    >  2:1057360/1‑150 (MQ=255)
                                                                                                                   cTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGcc                                    >  2:882121/1‑150 (MQ=255)
                                                                                                                    tCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            <  2:297894/94‑1 (MQ=255)
                                                                                                                    tCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGa                                                                                            >  1:297894/1‑94 (MQ=255)
                                                                                                                    tCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGAAAGCTTGGCGTGCCGCACCTG‑GTACTGGAg                                                                                           <  1:1152135/95‑1 (MQ=255)
                                                                                                                    tCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGAAAGCTTGGCGTGCCGCACCTG‑GTACTGGAg                                                                                           >  2:1152135/1‑95 (MQ=255)
                                                                                                                      g‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAc                                 >  2:2316336/1‑150 (MQ=255)
                                                                                                                         aaTCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTa                                                              <  1:1527183/119‑1 (MQ=255)
                                                                                                                         aaTCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTa                                                              >  2:1527183/1‑119 (MQ=255)
                                                                                                                           tCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcc                             >  2:2297415/1‑150 (MQ=255)
                                                                                                                              acacCCT‑CGTCGT‑CGGCGCCGGCGAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcccgg                          >  1:921792/1‑150 (MQ=255)
                                                                                                                               cacCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcccggc                         <  1:897468/150‑1 (MQ=255)
                                                                                                                                  ccT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTg                      >  1:680849/1‑150 (MQ=255)
                                                                                                                                  ccT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTg                      >  2:2027820/1‑150 (MQ=255)
                                                                                                                                           t‑cgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTg                                                                                                              >  2:1130511/1‑56 (MQ=255)
                                                                                                                                           t‑cgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTg                                                                                                              <  1:1130511/56‑1 (MQ=255)
                                                                                                                                           t‑cgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                              >  2:1858930/1‑118 (MQ=255)
                                                                                                                                           t‑cgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                              <  1:1858930/118‑1 (MQ=255)
                                                                                                                                             ccacgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCg             >  1:327770/4‑150 (MQ=255)
                                                                                                                                                  ccGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACAGCCCGGCCTGGCACGACCGCTTcc        <  2:273112/150‑1 (MQ=255)
                                                                                                                                                    ggCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTCGAGTG‑GGGGCGTGAGTTCTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCg      <  2:1814293/150‑1 (MQ=255)
                                                                                                                                                     gCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCgg     <  1:2059916/150‑1 (MQ=255)
                                                                                                                                                       cAGGCCGGGGTGGCGAGGAGCGAACGCCGGAGCAAGCTTGGCGTGCGGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCGGGCGTACCGGCCGCTGGGACGCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTACCgggg   <  2:405460/150‑1 (MQ=255)
                                                                                                                                                       cAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCgggg   <  2:1362137/150‑1 (MQ=255)
                                                                                                                                                       cAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCgggg   >  1:684773/1‑150 (MQ=255)
                                                                                                                                                       cAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCgggg   >  1:1362137/1‑150 (MQ=255)
                                                                                                                                                        aGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGGGc  >  1:529793/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
ATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCG‑AAATCGACACCCT‑CGTCGT‑CGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTG‑GTACTGGAG‑CGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGGGC  >  NC_002947/4053371‑4053669

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