Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,052 T→A 5.7% intergenic (‑119/+23) 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,0520TA 5.7% 493.2 / 14.5 141intergenic (‑119/+23)sseB/pepBrhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (66/67);  new base A (2/6);  total (68/73)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 144 total aligned reads displayed.

GTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTT  >  NC_000913/2654904‑2655200
                                                                                                                                                    |                                                                                                                                                    
gTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                                      >  1:17883/1‑149 (MQ=255)
gTCTTCAAGTTCGTTTTTTGTATCGGACAAAAGAATCACAATTCATTCAACGTTCGGGGCGCAAATGCATGATGCGCTAAGCTTAACAGTCCTACAAGGGATTCGCAATTTGATGAATTGGCTTAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                                      >  4:905145/1‑149 (MQ=37)
    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)
                  tGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatat                                                                                                                                    <  2:290092/149‑1 (MQ=255)
                    tttCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTa                                                                                                                                  >  2:1127086/1‑149 (MQ=255)
                      tCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGt                                                                                                                                <  4:333263/149‑1 (MQ=255)
                       cGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGtt                                                                                                                               <  1:794194/149‑1 (MQ=255)
                         gACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTa                                                                                                                             >  2:1399983/1‑149 (MQ=255)
                          aCATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTAc                                                                                                                            >  3:1230783/1‑149 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                                 <  2:568797/139‑1 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                                 >  1:568797/1‑139 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                       >  2:154534/1‑149 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                       <  1:865264/149‑1 (MQ=255)
                                gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  2:1325717/131‑1 (MQ=25)
                                gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                      <  2:1057554/149‑1 (MQ=255)
                                   ccACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAc                                                                                                                   >  2:312381/1‑149 (MQ=255)
                                    cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                                  <  3:1196506/149‑1 (MQ=255)
                                    cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                                  >  1:159175/1‑149 (MQ=255)
                                    cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                                  >  3:65703/1‑149 (MQ=255)
                                     acaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAg                                                                                                                 >  4:1378559/1‑149 (MQ=255)
                                       aaTTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                               <  3:130869/149‑1 (MQ=255)
                                       aaTTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                               >  1:1270747/1‑149 (MQ=255)
                                        attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                              >  1:731673/1‑149 (MQ=255)
                                        attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                              >  2:1023596/1‑149 (MQ=255)
                                            attcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCt                                                                                                          <  1:64884/149‑1 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                               <  2:734171/143‑1 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         >  2:861345/1‑149 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         >  1:1059520/1‑149 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         >  4:583635/1‑149 (MQ=255)
                                                aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCg                                                                                                      >  2:230968/1‑149 (MQ=255)
                                                 aCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                                     >  3:503588/1‑149 (MQ=255)
                                                            gcAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCaa                                                                                          <  4:1230783/149‑1 (MQ=255)
                                                            gcAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCaa                                                                                          <  1:1399983/149‑1 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        >  4:214319/1‑101 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGt                                                                                        >  4:411740/1‑149 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATGTGTTGAATTTGCATAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGt                                                                                        >  2:930587/1‑149 (MQ=255)
                                                                aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACACCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTcc                                                                                      <  2:560683/149‑1 (MQ=255)
                                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                                     <  1:1405227/149‑1 (MQ=255)
                                                                    cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                  <  1:1023596/149‑1 (MQ=255)
                                                                    cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                  <  4:204116/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)
                                                                            cTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGc                                                                          <  4:733961/149‑1 (MQ=255)
                                                                            cTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGc                                                                          >  1:412566/1‑149 (MQ=255)
                                                                               cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGa                                                                       <  3:1424702/149‑1 (MQ=255)
                                                                                  ttATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                    >  1:563263/1‑149 (MQ=255)
                                                                                   tATCAGGCCTACAAGGGATCCGCAATTTGTTGAGTTTGCAGAATTTGTAGGCCGGCTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAAAAATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                   >  2:849998/1‑149 (MQ=38)
                                                                                    aTCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACt                                                                  >  2:690369/1‑149 (MQ=255)
                                                                                       aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgc                                                                                                   <  2:1265120/113‑1 (MQ=255)
                                                                                       aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGtt                                                               >  4:2351/1‑149 (MQ=255)
                                                                                        ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                              >  4:62371/1‑149 (MQ=255)
                                                                                               aaGTGATTCGCAATATCTTGAATTTGCAGAATTTGTAGACCAGATAAGGCGTTTACGCCGCACCCGGCATAGTCGTTTACGTCGTTAACAGATTAGCTATCGTGCGCAGACCAAGTCCCGACGCGCCCGCAGACCACTGTGCAACCGgc                                                       <  1:1469257/149‑1 (MQ=25)
                                                                                               aaGGGCTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGgc                                                       >  4:796089/1‑149 (MQ=255)
                                                                                               aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                  >  1:1167678/1‑122 (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)
                                                                                                        gCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACgg                                              <  2:17883/149‑1 (MQ=255)
                                                                                                           aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                            <  4:1060801/132‑1 (MQ=255)
                                                                                                           aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                            >  3:1060801/1‑132 (MQ=255)
                                                                                                             ttgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGt                                         <  2:244321/149‑1 (MQ=255)
                                                                                                              tgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAATGTTCAACCGGCGCTTTACGGTAAGTc                                        <  4:416055/149‑1 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                    >  4:829022/1‑116 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                    <  3:829022/116‑1 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                   <  4:322241/149‑1 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                   <  3:1332788/149‑1 (MQ=255)
                                                                                                                         cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                             <  3:727377/149‑1 (MQ=255)
                                                                                                                         cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                             >  2:165823/1‑149 (MQ=255)
                                                                                                                          aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGATCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                            >  1:1209650/1‑149 (MQ=255)
                                                                                                                          aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                            >  1:1297473/1‑149 (MQ=255)
                                                                                                                             aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                         <  2:112640/149‑1 (MQ=255)
                                                                                                                             aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                         <  1:304322/149‑1 (MQ=255)
                                                                                                                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGatat                        >  2:1482221/1‑149 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                      <  2:685809/100‑1 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                      >  1:685809/1‑100 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     <  2:536155/149‑1 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcag                    >  4:508793/1‑149 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                      >  3:774529/1‑49 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                      <  4:774529/49‑1 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                 <  3:471464/149‑1 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag                <  4:1309508/149‑1 (MQ=255)
                                                                                                                                        ggatgcggCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTaa                                           <  3:152599/114‑1 (MQ=38)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     <  4:1106316/141‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     >  3:1106316/1‑141 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCt             <  2:331009/149‑1 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                   <  2:721951/92‑1 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                   >  1:721951/1‑92 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc          <  2:1297473/149‑1 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt         <  3:821172/149‑1 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGa       <  2:1270747/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa     >  4:1245062/1‑149 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa     <  2:81939/149‑1 (MQ=255)
                                                                                                                                                    tACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGtt  <  1:1127086/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTT  >  NC_000913/2654904‑2655200

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: