Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,228 T→A 6.4% L153L (CTT→CTA) ‡ hyaD → putative hydrogenase 1 maturation protease HyaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,2280TA 6.4% 250.4 / 5.5 78L153L (CTT→CTA) ‡hyaDputative hydrogenase 1 maturation protease HyaD
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base A (2/3);  total (39/39)
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.94e-01

GGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGAT  >  NC_000913/1036081‑1036377
                                                                                                                                                   |                                                                                                                                                     
ggCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgc                                                                                                                                                     <  1:1334042/150‑1 (MQ=255)
  tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctg                                                                                                                                                   >  2:1771303/2‑150 (MQ=255)
  tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctg                                                                                                                                                   >  1:1208016/2‑150 (MQ=255)
   gcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGGCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctgc                                                                                                                                                  <  2:1340668/150‑1 (MQ=255)
    cTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCg                                                                                                                                                 >  1:89042/1‑150 (MQ=255)
       gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTgg                                                                                                                                              <  1:172471/150‑1 (MQ=255)
       gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgcca                                                                                                                                              <  2:89042/150‑4 (MQ=255)
        cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag                                                                                                                                             <  2:20484/150‑5 (MQ=255)
        cGGATAGCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag                                                                                                                                             <  1:911810/150‑5 (MQ=255)
             atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAtt                                                                                                                                        >  2:402357/1‑150 (MQ=255)
                 gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGc                                                                                                                                    >  2:145806/1‑150 (MQ=255)
                 gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGc                                                                                                                                    <  1:1226242/150‑1 (MQ=255)
                    gACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc                                                                                                                                 >  2:1088697/1‑150 (MQ=255)
                         cTGCCAGCACATCTTGCCCTCGTCGGGCTGCAACCCGCAATGCTCGACGACTCCGGCGGTAGCCTGAGCGATCTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCg                                                                                                                            <  2:365424/150‑1 (MQ=255)
                           gCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc                                                                                                                          <  1:1944543/150‑1 (MQ=255)
                             cAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTa                                                                                                                        >  2:1210815/1‑150 (MQ=255)
                               gCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAt                                                                                                                      <  1:602236/150‑1 (MQ=255)
                                cacaTATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATg                                                                                                                     <  1:242138/150‑1 (MQ=255)
                                  caTATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGaa                                                                                                                   >  1:22668/1‑150 (MQ=255)
                                  caTATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGaa                                                                                                                   >  1:1307560/1‑150 (MQ=255)
                                    tatTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATc                                                                                                                 >  2:246266/1‑150 (MQ=255)
                                    tatTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATc                                                                                                                 >  1:588811/1‑150 (MQ=255)
                                          ccTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTg                                                                                                           >  1:1552541/1‑150 (MQ=255)
                                             cgtcgGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGtct                                                                                                        >  2:391156/1‑150 (MQ=255)
                                                cgGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCaa                                                                                                     >  1:194151/1‑150 (MQ=255)
                                                  gTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAAtt                                                                                                   >  2:53209/1‑150 (MQ=255)
                                                  gTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAAtt                                                                                                   >  2:890034/1‑150 (MQ=255)
                                                    cTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTAt                                                                                                 <  1:1771303/150‑1 (MQ=255)
                                                     tGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACANGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCATNTATg                                                                                                <  2:1437887/150‑1 (MQ=255)
                                                     tGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATg                                                                                                >  2:1884920/1‑150 (MQ=255)
                                                     tGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATg                                                                                                >  2:1891719/1‑150 (MQ=255)
                                                     tGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGAAGCTTGCTGCGTGGGGAATTGTGCCGCAAACGGCTAATGAATCGCGCTGTCTCATTTATg                                                                                                <  2:1430786/150‑1 (MQ=255)
                                                      gCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGTCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGa                                                                                               <  1:399376/150‑1 (MQ=255)
                                                      gCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGa                                                                                               <  1:1031468/150‑1 (MQ=255)
                                                      gCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGa                                                                                               >  1:780571/1‑150 (MQ=255)
                                                        aaCCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGAct                                                                                             >  2:692244/1‑150 (MQ=255)
                                                        aaCCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGAct                                                                                             >  1:1244333/1‑150 (MQ=255)
                                                          cccGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGActgt                                                                                           >  2:1118329/1‑150 (MQ=255)
                                                           ccGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGActgtc                                                                                          <  2:588811/150‑1 (MQ=255)
                                                              cAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGActgtctgt                                                                                       >  2:1128625/1‑150 (MQ=255)
                                                                 tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGa                                                                                    >  1:623426/1‑150 (MQ=255)
                                                                  gCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGAt                                                                                   >  2:1531104/1‑150 (MQ=255)
                                                                   cTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATg                                                                                  >  1:187104/1‑150 (MQ=255)
                                                                    tCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATgg                                                                                 <  2:194151/150‑1 (MQ=255)
                                                                       acgacTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGaaa                                                                              <  1:145806/150‑1 (MQ=255)
                                                                          acTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAAtt                                                                           >  2:908935/1‑150 (MQ=255)
                                                                            tACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTAc                                                                         <  1:1673406/150‑1 (MQ=255)
                                                                             aCGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGTCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACg                                                                        <  1:1012119/150‑1 (MQ=255)
                                                                             aCGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACg                                                                        >  1:899342/1‑150 (MQ=255)
                                                                             aCGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACg                                                                        <  2:22668/150‑1 (MQ=255)
                                                                             aCGGCGGTAGCCCGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACg                                                                        >  2:207920/1‑150 (MQ=255)
                                                                               ggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGaa                                                                      >  2:239401/1‑150 (MQ=255)
                                                                               ggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGaa                                                                      <  1:1071588/150‑1 (MQ=255)
                                                                                gcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAg                                                                     <  2:1591408/150‑1 (MQ=255)
                                                                                 cggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAgg                                                                    <  1:1552488/150‑1 (MQ=255)
                                                                                 cggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAgg                                                                    <  1:1088697/150‑1 (MQ=255)
                                                                                   gTAGCCCGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCg                                                                  <  1:207920/150‑1 (MQ=255)
                                                                                       ccTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCg                                                              <  1:1913283/150‑1 (MQ=255)
                                                                                       ccTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCg                                                              <  1:692244/150‑1 (MQ=255)
                                                                                               aaCTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGcc                                                      <  1:1118329/150‑1 (MQ=255)
                                                                                                 cTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAg                                                    <  1:1128625/150‑1 (MQ=255)
                                                                                                 cTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAg                                                    >  2:511767/1‑150 (MQ=255)
                                                                                                     cACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTAcc                                                <  1:1531104/150‑1 (MQ=255)
                                                                                                      aCGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCg                                               >  1:876471/1‑150 (MQ=255)
                                                                                                            gCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGac                                         <  2:1307560/150‑1 (MQ=255)
                                                                                                              aaCTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGacac                                       <  1:195091/150‑1 (MQ=255)
                                                                                                                 tGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACagg                                    <  1:908935/150‑1 (MQ=255)
                                                                                                                  gCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACagga                                   >  1:1731645/1‑150 (MQ=255)
                                                                                                                           ggAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCa                             <  1:549019/147‑1 (MQ=255)
                                                                                                                           ggAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCa                             >  2:549019/1‑147 (MQ=255)
                                                                                                                           ggAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAggg                          >  1:1137456/1‑150 (MQ=255)
                                                                                                                                aGGCGGCGCTGGCGCAGCTTGCTGCGTGTGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAg                     <  1:1705988/150‑1 (MQ=255)
                                                                                                                                 ggcggcGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGc                    >  1:208571/1‑150 (MQ=255)
                                                                                                                                   cggcgCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCaa                  >  2:1417662/1‑150 (MQ=255)
                                                                                                                                     gcgcTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACg                <  1:1210815/150‑1 (MQ=255)
                                                                                                                                      tccccacgcAGCAAGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGa               >  2:153196/7‑150 (MQ=255)
                                                                                                                                      cgcTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGa               >  1:1389740/1‑150 (MQ=255)
                                                                                                                                       gcTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAg                     >  1:1749411/1‑143 (MQ=255)
                                                                                                                                       gcTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAg                     <  2:1749411/143‑1 (MQ=255)
                                                                                                                                            cgcAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGcc         >  1:309264/1‑150 (MQ=255)
                                                                                                                                               aGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCAtt      <  1:1891719/150‑1 (MQ=255)
                                                                                                                                               aGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAAGCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCAtt      <  1:1884920/150‑1 (MQ=255)
                                                                                                                                                   tgctgcGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGAt  >  2:1587373/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGAT  >  NC_000913/1036081‑1036377

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