Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,045 A→C 5.8% intergenic (‑112/+30) sseB ← / ← pepB rhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,655,0450AC 5.8% 480.7 / 18.0 137intergenic (‑112/+30)sseB/pepBrhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base A (68/61);  new base C (3/5);  total (71/66)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 143 total aligned reads displayed.

CCAGCAGGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCT  >  NC_000913/2654897‑2655193
                                                                                                                                                    |                                                                                                                                                    
ccAGCAGGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATaa                                                                                                                                                      <  4:382128/149‑1 (MQ=255)
 cagcagGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAg                                                                                                                                                     <  2:922354/149‑1 (MQ=255)
 cagcagGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAg                                                                                                                                                     <  3:739584/149‑1 (MQ=255)
   gcagGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGc                                                                                                                                                   >  4:518833/1‑149 (MQ=255)
    cagGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCg                                                                                                                                                  <  3:518833/149‑1 (MQ=255)
     agGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGt                                                                                                                                                 >  1:560497/1‑149 (MQ=255)
     agGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGt                                                                                                                                                 >  2:884733/1‑149 (MQ=255)
       gTCTTCAAGTTCGTTTTTTGTATCGGACAAAAGAATCACAATTCATTCAACGTTCGGGGCGCAAATGCATGATGCGCTAAGCTTAACAGTCCTACAAGGGATTCGCAATTTGATGAATTGGCTTAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                               >  4:905145/1‑149 (MQ=37)
         cTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAc                                                                                                                                             <  4:211006/149‑1 (MQ=255)
           tCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                           <  1:253132/149‑1 (MQ=255)
            cAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                          <  2:405680/149‑1 (MQ=255)
             aaGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         >  1:81939/1‑149 (MQ=255)
                ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTATGCCGGATGCGGCGTAAAcgccttat                                                                                                                                      <  3:302631/149‑5 (MQ=38)
                     tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 <  3:195804/149‑1 (MQ=37)
                     tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                 <  1:894160/149‑1 (MQ=255)
                      ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                <  3:1291817/149‑1 (MQ=255)
                      ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                <  1:607181/149‑2 (MQ=37)
                         tGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatat                                                                                                                             <  2:290092/149‑1 (MQ=255)
                              cGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGtt                                                                                                                        <  1:794194/149‑1 (MQ=255)
                                gACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTa                                                                                                                      >  2:1399983/1‑149 (MQ=255)
                                      aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                          >  1:568797/1‑139 (MQ=255)
                                      aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                          <  2:568797/139‑1 (MQ=255)
                                      aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                <  1:865264/149‑1 (MQ=255)
                                       gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 >  1:1325717/1‑131 (MQ=25)
                                       gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                               >  4:25623/1‑149 (MQ=255)
                                         aCCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTaa                                                                                                             >  2:1469257/1‑149 (MQ=255)
                                          ccACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAc                                                                                                            >  2:312381/1‑149 (MQ=255)
                                           cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                           <  3:1196506/149‑1 (MQ=255)
                                           cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                           >  1:159175/1‑149 (MQ=255)
                                              aaTTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                        <  3:130869/149‑1 (MQ=255)
                                               attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                       >  2:1023596/1‑149 (MQ=255)
                                                ttcattcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTa                                                                                                      <  2:111259/149‑1 (MQ=255)
                                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                        <  2:734171/143‑1 (MQ=255)
                                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                        >  1:734171/1‑143 (MQ=255)
                                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                  >  4:583635/1‑149 (MQ=255)
                                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                  >  1:1059520/1‑149 (MQ=255)
                                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                  >  1:244321/1‑149 (MQ=255)
                                                       aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAgg                                                                                                                                                    >  1:1250350/1‑96 (MQ=255)
                                                       aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAgg                                                                                                                                                    <  2:1250350/96‑1 (MQ=255)
                                                       aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCg                                                                                               >  4:356270/1‑149 (MQ=255)
                                                       aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCg                                                                                               >  2:230968/1‑149 (MQ=255)
                                                        aCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                              >  1:463612/1‑149 (MQ=255)
                                                        aCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                              >  3:503588/1‑149 (MQ=255)
                                                                       aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACACCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTcc                                                                               <  2:560683/149‑1 (MQ=255)
                                                                        tGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                              <  1:1405227/149‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           <  3:574798/149‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           <  4:204116/149‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           <  1:1023596/149‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           >  4:727377/1‑149 (MQ=255)
                                                                            tGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAg                                                                          >  1:536155/1‑149 (MQ=255)
                                                                                 cgcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGccc                                                                     <  2:87715/149‑1 (MQ=255)
                                                                                  gcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCg                                                                    >  3:1125906/1‑149 (MQ=255)
                                                                                      cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGa                                                                <  3:1424702/149‑1 (MQ=255)
                                                                                      cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGa                                                                <  2:850591/149‑1 (MQ=255)
                                                                                       gCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                               >  2:1312341/1‑149 (MQ=255)
                                                                                          tATCAGGCCTACAAGGGATCCGCAATTTGTTGAGTTTGCAGAATTTGTAGGCCGGCTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAAAAATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                            >  2:849998/1‑149 (MQ=38)
                                                                                           aTCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACt                                                           >  2:690369/1‑149 (MQ=255)
                                                                                              aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgc                                                                                            >  1:1265120/1‑113 (MQ=255)
                                                                                              aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGtt                                                        >  4:2351/1‑149 (MQ=255)
                                                                                               ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                       <  2:18283/149‑1 (MQ=255)
                                                                                                      aaGTGATTCGCAATATCTTGAATTTGCAGAATTTGTAGACCAGATAAGGCGTTTACGCCGCACCCGGCATAGTCGTTTACGTCGTTAACAGATTAGCTATCGTGCGCAGACCAAGTCCCGACGCGCCCGCAGACCACTGTGCAACCGgc                                                <  1:1469257/149‑1 (MQ=25)
                                                                                                      aaGGGCTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGgc                                                >  4:796089/1‑149 (MQ=255)
                                                                                                      aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           >  1:1167678/1‑122 (MQ=255)
                                                                                                      aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                           <  2:1167678/122‑1 (MQ=255)
                                                                                                      aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGgc                                                <  4:1209302/149‑1 (MQ=255)
                                                                                                        gggATTCGCAATTTGTTGAATTTGCAGGATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCTGTAGCGCCCGCAGACCACTGTTCAACCGgcgc                                              >  2:678828/1‑149 (MQ=255)
                                                                                                           aTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCttt                                           >  2:1083175/1‑149 (MQ=255)
                                                                                                            ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                          >  3:416055/1‑149 (MQ=255)
                                                                                                             tCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTAc                                         <  1:348031/149‑1 (MQ=255)
                                                                                                                  aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                     >  3:1060801/1‑132 (MQ=255)
                                                                                                                  aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                     <  4:1060801/132‑1 (MQ=255)
                                                                                                                     tgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAATGTTCAACCGGCGCTTTACGGTAAGTc                                 <  4:416055/149‑1 (MQ=255)
                                                                                                                       ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGc                               >  2:304322/1‑149 (MQ=255)
                                                                                                                       ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGc                               >  4:465876/1‑149 (MQ=255)
                                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                             >  4:829022/1‑116 (MQ=255)
                                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                             <  3:829022/116‑1 (MQ=255)
                                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                            <  3:1332788/149‑1 (MQ=255)
                                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                            <  4:322241/149‑1 (MQ=255)
                                                                                                                           aTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGTCGCATCCGGCATATTAGTTTACGCCGTTAGCAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCCCCCGCAGACCACTGTTCAACCGTCGCTTTACCGTAAGTCGCCgtg                           >  4:1261054/1‑147 (MQ=39)
                                                                                                                                cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                      <  3:727377/149‑1 (MQ=255)
                                                                                                                                cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                      >  2:165823/1‑149 (MQ=255)
                                                                                                                                 aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                     >  1:1297473/1‑149 (MQ=255)
                                                                                                                                    aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                  <  2:112640/149‑1 (MQ=255)
                                                                                                                                     tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGatat                 >  2:1482221/1‑149 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                               >  1:685809/1‑100 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                               <  2:685809/100‑1 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca              <  2:536155/149‑1 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcag             >  4:508793/1‑149 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                               <  4:774529/49‑1 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca          >  2:103443/1‑149 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca          <  4:65703/149‑1 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag         <  4:1309508/149‑1 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTCAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag         >  3:1024026/1‑149 (MQ=255)
                                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc       <  2:159175/149‑1 (MQ=255)
                                                                                                                                                gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca              <  4:1106316/141‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCACGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc   <  1:312381/149‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                            >  1:721951/1‑92 (MQ=255)
                                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                            <  2:721951/92‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc   <  2:1297473/149‑1 (MQ=255)
                                                                                                                                                    aGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt  <  3:821172/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCAGGTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCT  >  NC_000913/2654897‑2655193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: