Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,623 T→A 6.1% K183I (AAA→ATA) ‡ fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,6230TA 6.1% 313.0 / 8.5 98K183I (AAA→ATA) ‡fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (49/43);  new base A (3/3);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGCCGCCCG  >  NC_002947/243479‑243767
                                                                                                                                                |                                                                                                                                                
gagGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGGCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGgcg                                                                                                                                             >  2:568114/1‑150 (MQ=255)
    gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                                                                                                                                         <  1:590940/150‑3 (MQ=255)
      gcagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgtg                                                                                                                                       <  2:1143676/150‑1 (MQ=255)
            cagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtt                                                                                                                                 <  1:2149551/150‑1 (MQ=255)
             agGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAttt                                                                                                                                >  1:1729722/1‑150 (MQ=255)
              ggCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtttt                                                                                                                               <  2:756785/150‑1 (MQ=255)
               gCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGTTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTc                                                                                                                              <  1:1496776/150‑1 (MQ=255)
                ccTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCg                                                                                                                             <  1:40610/150‑1 (MQ=255)
                 cTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCgg                                                                                                                            >  2:1337107/1‑150 (MQ=255)
                   gcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTAGGCa                                                                                                                          >  2:1811919/3‑150 (MQ=255)
                      gTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCACCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAgg                                                                                                                       >  2:1007550/1‑150 (MQ=255)
                         ggCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGa                                                                                                                    <  2:1959851/150‑1 (MQ=255)
                            gTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCt                                                                                                                 >  2:542402/1‑150 (MQ=255)
                            gTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                                                                                                                                         >  2:1212346/1‑124 (MQ=255)
                            gTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                                                                                                                                         <  1:1212346/126‑3 (MQ=255)
                               aGTCGCCAACCGCCAGCAGTTCGATATCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGgcgc                                                                                                                                            >  1:194456/1‑120 (MQ=255)
                               aGTCGCCAACCGCCAGCAGTTCGATATCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGgcgc                                                                                                                                            <  2:194456/120‑1 (MQ=255)
                                   gCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGcc                                                                                                          >  2:1295890/1‑150 (MQ=255)
                                     cAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAtt                                                                                                                                                  >  2:1379940/1‑108 (MQ=255)
                                     cAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAtt                                                                                                                                                  <  1:1379940/108‑1 (MQ=255)
                                     cAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCtt                                                                                                        <  2:898724/150‑1 (MQ=255)
                                        ccGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGa                                                                                                     >  2:883334/1‑150 (MQ=255)
                                        ccGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGa                                                                                                     <  1:1337107/150‑1 (MQ=255)
                                         cGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGAt                                                                                                    >  2:1655511/1‑150 (MQ=255)
                                            cagcagTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAAc                                                                                                 <  1:344431/150‑1 (MQ=255)
                                            cagcagTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAAc                                                                                                 <  1:344698/150‑1 (MQ=255)
                                                agTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGaa                                                                                             <  1:2305024/150‑1 (MQ=255)
                                                  ttCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATc                                                                                           <  1:863336/150‑1 (MQ=255)
                                                     gATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAg                                                                                        <  2:967773/150‑1 (MQ=255)
                                                         tcGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTg                                                                                    >  1:431593/1‑150 (MQ=255)
                                                           gAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCtt                                                                                                        >  2:596111/1‑128 (MQ=255)
                                                           gAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGaa                                                                                  <  2:291377/150‑1 (MQ=255)
                                                           gAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAAGTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCtt                                                                                                        <  1:596111/128‑1 (MQ=255)
                                                            aaTACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAAc                                                                                 <  1:68772/150‑1 (MQ=255)
                                                                 cAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAg                                                                                        <  1:341232/138‑1 (MQ=255)
                                                                 cAAGGCGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAg                                                                                        >  2:341232/1‑138 (MQ=255)
                                                                  aaGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTgg                                                                           <  2:480741/150‑1 (MQ=255)
                                                                   aGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTggg                                                                          <  1:883334/150‑1 (MQ=255)
                                                                     ggggTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgt                                                                                                                                        >  2:828581/1‑86 (MQ=255)
                                                                     ggggTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgt                                                                                                                                        <  1:828581/86‑1 (MQ=255)
                                                                       ggTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATACAGGTTGAACCATTGGGGGGt                                                                      >  2:989921/1‑150 (MQ=255)
                                                                        gtgtAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtc                                                                     <  1:1655511/150‑1 (MQ=255)
                                                                           tAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctgg                                                                  >  2:1150632/1‑150 (MQ=255)
                                                                             gggTGCACACAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggtt                                                                <  1:103150/150‑1 (MQ=255)
                                                                               gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCTATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttct                                                              <  2:258158/150‑1 (MQ=255)
                                                                               gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTgggg                                                                         <  1:1146555/139‑1 (MQ=255)
                                                                               gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTgggg                                                                         >  2:1146555/1‑139 (MQ=255)
                                                                               gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttct                                                              >  1:121762/1‑150 (MQ=255)
                                                                               gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttct                                                              <  1:1372915/150‑1 (MQ=255)
                                                                                tGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAAGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTACCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                             <  2:424107/150‑1 (MQ=255)
                                                                                   aaTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTc                                                          >  1:73235/1‑150 (MQ=255)
                                                                                   aaTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTc                                                          >  2:615767/1‑150 (MQ=255)
                                                                                   aaTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTc                                                          >  2:1713938/1‑150 (MQ=255)
                                                                                    aTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCg                                                         >  1:1893268/1‑150 (MQ=255)
                                                                                     tCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGa                                                        >  2:269741/1‑150 (MQ=255)
                                                                                     tCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGa                                                        >  1:546791/1‑150 (MQ=255)
                                                                                     tCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGa                                                        >  1:546697/1‑150 (MQ=255)
                                                                                     tCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGa                                                        >  2:1533825/1‑150 (MQ=255)
                                                                                     tCAGGTCTCCCGCCCCCTCGGCGCCATAGGATTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGa                                                        >  1:1280011/1‑150 (MQ=255)
                                                                                     gcAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGc                                                                                                               <  1:135547/94‑1 (MQ=255)
                                                                                     gcAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGc                                                                                                               >  2:135547/2‑95 (MQ=255)
                                                                                      cAGGTCTCCCGCCCCCTCGGCGCCATAGGATTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGaa                                                       <  2:1280011/150‑1 (MQ=255)
                                                                                         gTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTg                                                                                                                >  2:1472475/1‑90 (MQ=255)
                                                                                         gTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTg                                                                                                                <  1:1472475/90‑1 (MQ=255)
                                                                                          tctcCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATAtt                                                   >  2:152689/1‑150 (MQ=255)
                                                                                            tcCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGgc                                                                                                                                              >  1:1920830/1‑57 (MQ=255)
                                                                                            tcCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGgc                                                                                                                                              <  2:1920830/57‑1 (MQ=255)
                                                                                               cGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                                                                                                                                         <  1:462412/59‑3 (MQ=38)
                                                                                               cGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                                                                                                                                         >  2:462412/1‑57 (MQ=255)
                                                                                                    ccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTg                                                                            <  1:2267263/115‑1 (MQ=255)
                                                                                                    ccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTg                                                                            >  2:2267263/1‑115 (MQ=255)
                                                                                                      tCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGGCTGGTTCTGGTCGAATATTGAGCCGTCCGGc                                       >  1:1693564/1‑150 (MQ=255)
                                                                                                        ggCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCGAGCCTTCGATAACCGAATCCAGGTTGAAc                                                                                 >  1:1852802/1‑106 (MQ=255)
                                                                                                        ggCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAAc                                                                                 <  2:1852802/106‑1 (MQ=255)
                                                                                                        ggCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTgg                                                                           >  1:850899/1‑112 (MQ=255)
                                                                                                        ggCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTgg                                                                           <  2:850899/112‑1 (MQ=255)
                                                                                                           gcCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCACCCAGGTAGAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCt                                  >  2:1399252/1‑150 (MQ=255)
                                                                                                             cATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGc                                <  1:269741/150‑1 (MQ=255)
                                                                                                             cATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGc                                >  2:922023/1‑150 (MQ=255)
                                                                                                             cATAAGATTGCGCTGAAGAAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGc                                >  2:760682/1‑150 (MQ=255)
                                                                                                              aTAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGcc                               <  2:431593/150‑1 (MQ=255)
                                                                                                               tAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGccc                              >  2:1645449/1‑150 (MQ=255)
                                                                                                                aGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCca                             <  2:73235/150‑1 (MQ=255)
                                                                                                                    ttGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATa                         >  1:439166/1‑150 (MQ=255)
                                                                                                                        gcTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGaa                                                                                             >  1:812042/1‑78 (MQ=255)
                                                                                                                        gcTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGaa                                                                                             <  2:812042/78‑1 (MQ=255)
                                                                                                                        gcTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                             >  1:2014449/1‑110 (MQ=255)
                                                                                                                        gcTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                             <  2:2014449/110‑1 (MQ=255)
                                                                                                                           gAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCgcact                  >  1:302360/1‑150 (MQ=255)
                                                                                                                            aaGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGcc                               <  1:1768185/136‑1 (MQ=255)
                                                                                                                            aaGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGcc                               >  2:1768185/1‑136 (MQ=255)
                                                                                                                              ggAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCgcacttgc               >  1:2699/1‑150 (MQ=255)
                                                                                                                                   acaccggCGCCAAATGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACttt          >  2:324696/6‑150 (MQ=255)
                                                                                                                                     ccAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgc        >  1:1564444/1‑150 (MQ=255)
                                                                                                                                      cAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTc                                           <  1:1336063/114‑1 (MQ=255)
                                                                                                                                      cAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTc                                           >  2:1336063/1‑114 (MQ=255)
                                                                                                                                      cAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgcc       <  2:2699/150‑1 (MQ=255)
                                                                                                                                      cAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgcc       <  1:2136035/150‑1 (MQ=255)
                                                                                                                                      cAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgcc       >  1:2119392/1‑150 (MQ=255)
                                                                                                                                           gcCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGTTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGccgcccg  >  2:104546/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGCCGCCCG  >  NC_002947/243479‑243767

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